Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 65,857 12940 bp→GA [araD][setA] [araD], araA, araB, araC, yabI, thiQ, thiP, thiB, sroA, sgrR, sgrS, sgrT, [setA]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 65857 78796 12940 40 [0] [0] 40 [araD]–[setA] [araD],araA,araB,araC,yabI,thiQ,thiP,thiB,sroA,sgrR,sgrS,sgrT,[setA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 658560 (0.000)40 (0.660)
+GA
38/488 0.3 100% coding (695/696 nt) araD L‑ribulose‑5‑phosphate 4‑epimerase AraD
?NC_000913 78797 = 0 (0.000)coding (1177/1179 nt) setA sugar exporter SetA

TGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAG  >  NC_000913/78580‑79041
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
cgggaataaacgccacggactcttacggtgcaagcgtaaactgcaacggcccgttgtccgtcgccagccagaaggagacttctgtcccttgcggggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgcgtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCAt                                                                                                                                                                                                                     >  1:186018‑M2/218‑251 (MQ=255)
                gcaccttttcggtgcaagcgtaaccttcaacggcccgtgttcctcccccagccataaggagctttctgtcccttgcggggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAgggg                                                                                                                                                                                                     <  1:461142‑M2/50‑1 (MQ=255)
                  actcttgcggtgcaagcgtaacctgcaacggcccgttgtccgtcgccagccagaaggagacttctgtcccttgcggggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGAt                                                                                                                                                                                                   >  1:63384‑M2/200‑251 (MQ=255)
                     cttgcggtgcaagcgtaccctgcaacggcccgttgtccgtcgccagccagaaggagacttctgtcccttgcggggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGAtt                                                                                                                                                                                                  <  1:187393‑M2/53‑1 (MQ=255)
                       tgcggtgcaagcgtaacctgcaacggcccgttgtccgtcgccagccagaaggagacttctgtcccttgcggggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTa                                                                                                                                                                                                 <  1:444393‑M2/54‑1 (MQ=255)
                          ggtgcaagcgtaacctgcaacggcccgttgtccgtcgccagccagaaggagacttctgtcccttgcggggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctggggttttatacagtcattgaTGATGACGTGGACGGTAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCAt                                                                                                                                                                                           >  1:121769‑M2/192‑251 (MQ=255)
                            ggcaagcgtacccgccaacggccgtggtgccgtcgccggccggaaggagccttcggtcccttgcgggttgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                                                                                         <  2:134807‑M2/62‑1 (MQ=255)
                             gcaagcgtaacctgcaccggcccgttgtccgtcgccagccagaaggagacttctgtcccttgcggggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaaa                                                                                                                                                                                        <  1:262891‑M2/63‑1 (MQ=255)
                                     accctgcaacggcccgttgtccgtcgccagccagaaggagacttctgtcccttgcggggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGgtt                                                                                                                                                                                <  1:363728‑M2/71‑1 (MQ=255)
                                         ggcaccggcccgtggtccctcaccagccagaaggagacttctgtcccgtggggggtctcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGttt                                                                                                                                                                            <  2:388838‑M2/75‑1 (MQ=255)
                                            aacggcccgttgtccgtcgccagccagaaggagacttctgtcccttgcggggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacaggcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGc                                                                                                                                                                         <  1:28208‑M2/78‑1 (MQ=255)
                                                                    aagaaggagacttctctcctttgcgggttggcccgccagtgccgggttaagctaaaacctgcctgcgccacgcttaaaatccatcaaaaaaccatgtgtgagtatagcctggtttcgtttgattggcggtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCAGGTCATTTGACCGGGCATGGGTATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAAGgtcgtcgtc                                                                                                                                                 <  2:284311‑M2/102‑1 (MQ=255)
                                                                                          gcgtggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTTCTGCAAGGTAATCCCAATa                                                                                                                           <  1:406384‑M2/124‑1 (MQ=255)
                                                                                            ggggtgtcccgccagtgtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACt                                                                                                                         <  1:314562‑M2/126‑1 (MQ=255)
                                                                                            ggggtgtcccgccagtgtcgggttaagataaaacctgcccgcgccacgctgaaaatccatccaaaaacccggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGGTTAATTCATAAACGCAGGTTGGTTTGCTTTCTAAACGGCAATGGAGATGTCGGTCTTCAAGGTAAGCCCAataat                                                                                                                         >  2:336477‑M2/126‑249 (MQ=255)
                                                                                                 gtcccgccagtgtcgggttaaaataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggttttgtttgattggctgtggttttatacagtcattgaTGATGACGTGGGCGATAGCGGAAAGCCCGGTCATTTTACCGGGCAAGGGGGTTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTTCAAGGGAAGCCCAAAACTGTcc                                                                                                                     >  2:385471‑M2/121‑250 (MQ=255)
                                                                                                      gccagtgtcgggttaagataaaaactgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgggtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGGCGATAGCGGAAAGCCCGGTCAATTGACCGGGCAAGGGGATTAATTCATAAACGCCGGTTGTTTTTCTTCATAAGCGGCAATGGCGTCGACGTGCTGCCAGGGAAGCCCAAaactgtcccaaacg                                                                                                               >  1:319357‑M2/116‑245 (MQ=255)
                                                                                                            gtcgggttaagataaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTtcatc                                                                                                         >  1:481285‑M2/110‑251 (MQ=255)
                                                                                                               gggataagataatacatgcctgcgccacgctgagaatccatcaaaaaaccaggctagagtatagcctggtttcgtttgattagctgtggtattatacagtcattgaTGATGACGTGGACGATAGCGGAACGCCCGGTCCTTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCTTAAGTGGACATGGCGTCGTCGTACTGCAAGGGAAGCCCAATTCAGTCCAGAGCgttcatcacg                                                                                                      >  2:268791‑M2/107‑249 (MQ=255)
                                                                                                                    aagataaaacctgcctgcgccacgctgaaaatccatcaacaaaccaggcttgagtatagcctgatttcgtttgattggctgtggttttatacagtcattaaTGATGACGTGGACGCTAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGCTTAATTCATAAACGCAGGTTGTTTTGCTTCATCAGCGGCCATGGCGTCGTCTTGCTGCAAGGTAAGCCCAATACTCTCCAGACCGTTCATCCTGCCGTg                                                                                                 >  2:36975‑M2/102‑251 (MQ=255)
                                                                                                                         aaaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCAATTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAACCCAATAATGTCCAGACCGTTCATACTGCCGtggcggg                                                                                            >  2:197314‑M2/97‑250 (MQ=255)
                                                                                                                          aaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGgcggcg                                                                                           >  1:292325‑M2/96‑251 (MQ=255)
                                                                                                                          aaacctgcctgcgccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGgcggcg                                                                                           <  2:96760‑M2/156‑1 (MQ=255)
                                                                                                                                gcctgcgacacgctgaacacccatcaaaaaaccaggcttgagtatagcatggtttcgttggattggctggggttttatacgggcattgaTGATGTCGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTGTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAg                                                                                       <  2:320619‑M2/160‑1 (MQ=255)
                                                                                                                                      gcccagctgaaaatccatcaaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGAGAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGAt                                                                               <  2:256560‑M2/168‑1 (MQ=255)
                                                                                                                                      gccacgctgaaaatccatcaaaaaaccaggcttgagtatagcctgggttcggttgattggctgtggttttatacagtaattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGAt                                                                               <  2:204586‑M2/168‑1 (MQ=255)
                                                                                                                                                 aatccatcaaaaaaccaggcttgagtataggctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGAt                                                                    <  2:38156‑M2/179‑1 (MQ=255)
                                                                                                                                                         aaaaaaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTGTtc                                                            >  2:357523‑M2/65‑251 (MQ=255)
                                                                                                                                                          aaaaaccaggcttgagtatagcctggtttcgtttgatttgctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTtct                                                           <  2:57379‑M2/188‑1 (MQ=255)
                                                                                                                                                             aaccaggcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAAGCATCGATGGTAAAGCGATAGGTTTtctctc                                                        >  1:244084‑M2/61‑251 (MQ=255)
                                                                                                                                                               ccaggcttgatgatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTccc                                                      <  2:63384‑M2/193‑1 (MQ=255)
                                                                                                                                                                   gcttgagtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCtt                                                  <  1:101223‑M2/197‑1 (MQ=255)
                                                                                                                                                                        agtatagcctggtttcgtttgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCAcc                                             >  1:425383‑M2/50‑251 (MQ=255)
                                                                                                                                                                                       cgtttggttgtcggttgatttatacagtcattgaTGATGACGTTGACGATAGCGGAATGCCCGGTCAGTTGCACGGGCAAGGGGATTAATTCATACACGCAGGTTGTGTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGTAAGGCATCGATGGTAAAGCTATAGATTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGa                              <  2:129395‑M2/217‑1 (MQ=255)
                                                                                                                                                                                       cggtttattggctggggttttataccgtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCCGGTGGTTTTGCTTCATAAGCGGCAAGGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCAGCATGCAGTGGCGGCGGAAGGCATCGATGGGAAAGCGATAGGTTTTCGCGCCCGCTTTCACCTCTGGCGCTTCCAg                               <  2:371091‑M2/216‑1 (MQ=255)
                                                                                                                                                                                           tgattggctgtggttttatacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATcc                           <  1:147319‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                            atacagtcattgaTGATGACGTGGCCGCTAGCGGAAAGCACGGTCATTTGACAGGGCAAGGGGATTAATTCATAAACGCAGGATGTTTTGCTTCATAAGCGGCACAGGCGTCGTCGTGCTGCACGGTAAGCCCAATCCTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCCTCGATTGTCAAGCGATAGGTTTTATCTCCCGCTTACACCTCTTGCGCTTACAGAGCCACTTCGCAATGGATcccc         >  2:171712‑M2/14‑251 (MQ=255)
                                                                                                                                                                                                            atacagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTTCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATcccc         <  1:251921‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                              acagtcattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGGCGAAATGGATCCCCgg       >  1:219457‑M2/12‑251 (MQ=255)
                                                                                                                                                                                                                   cattgaTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAg  <  1:181153‑M2/245‑1 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
TGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAG  >  NC_000913/78580‑79041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: