Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 66,533 3540 bp→CCCGAGACCCGCACGCGTCGACCTGCA [araD]araB [araD], araA, araB

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 66533 70072 3540 269 [0] [0] 266 [araD]–araB [araD],araA,araB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 665320 (0.000)196 (1.050)
+CCCGAGACCCGCACGCGTCGACCTGCA
96/120 0.0 100% coding (19/696 nt) araD L‑ribulose‑5‑phosphate 4‑epimerase
?NC_000913 70073 = 0 (0.000)intergenic (‑25/‑314) araB/araC L‑ribulokinase/DNA‑binding transcriptional dual regulator

Only 100 of 266 total aligned reads displayed.

GGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGAC  >  NC_000913/69980‑70166
                                                                                             |                                                                                             
tttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcag                                                                                               <  1:6851041‑M2/1‑1 (MQ=255)
 ttgccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcagt                                                                                              >  1:5040955‑M2/93‑93 (MQ=255)
    ccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaataccgggcgcccgagacccgcacgcgtcgccctgcaggttc                                                                                           >  1:3302241‑M2/90‑92 (MQ=255)
    ccccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGtac                                                                                           >  1:3836484‑M2/90‑93 (MQ=255)
     cccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTata                                                                                          >  1:1280773‑M2/89‑93 (MQ=255)
      ccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTAtgc                                                                                         >  1:3716864‑M2/88‑93 (MQ=255)
      ccatgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTAtgc                                                                                         >  1:6219175‑M2/88‑93 (MQ=255)
        atgtgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATggaa                                                                                       >  1:6543320‑M2/86‑93 (MQ=255)
           tgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAgaac                                                                                    >  1:2074940‑M2/83‑93 (MQ=255)
           tgagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGgaa                                                                                    >  1:7085605‑M2/83‑94 (MQ=255)
            gagcgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAgaaaa                                                                                   >  1:2328432‑M2/82‑93 (MQ=255)
              acgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAg                                                                                 <  1:6457838‑M2/15‑1 (MQ=255)
              acgtgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAg                                                                                 <  1:6629514‑M2/15‑1 (MQ=255)
                ctgaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTa                                                                               <  1:3638186‑M2/17‑1 (MQ=255)
                  ggccaggttgtgttgtggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                             <  1:1857090‑M2/19‑1 (MQ=255)
                  gaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATgg                                                                                         >  2:3487371‑M2/76‑82 (MQ=255)
                  aaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                             <  1:2644536‑M2/19‑1 (MQ=255)
                  aaccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                             <  1:5444379‑M2/19‑1 (MQ=255)
                    ccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGtagagt                                                                           >  1:1752989‑M2/74‑93 (MQ=255)
                    ccaggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTagaga                                                                           <  1:4796333‑M2/21‑1 (MQ=255)
                      aggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGaa                                                                                     <  2:3254270‑M2/11‑1 (MQ=255)
                      aggttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGtagagaga                                                                         >  1:6957490‑M2/72‑93 (MQ=255)
                        gttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAAc                                                                                   >  2:1943922‑M2/70‑82 (MQ=255)
                        gttgtgctttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGCGAAAc                                                                                   >  2:91440‑M2/70‑82 (MQ=255)
                         ttgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGc                                                                      >  1:4239882‑M2/69‑94 (MQ=255)
                          tgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAg                                                                                 <  2:4584014‑M2/15‑1 (MQ=255)
                          tgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAg                                                                                 <  2:1129380‑M2/15‑1 (MQ=255)
                          tgtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTgct                                                                     >  1:1093820‑M2/68‑93 (MQ=255)
                           gtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGt                                                                                <  2:3873600‑M2/16‑1 (MQ=255)
                           gtgttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGGATGGAGAAACAGt                                                                                >  2:2466978‑M2/67‑82 (MQ=255)
                             gttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTag                                                                              >  2:62724‑M2/65‑82 (MQ=255)
                             gttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTag                                                                              <  2:1311828‑M2/18‑1 (MQ=255)
                              ttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTaga                                                                             >  2:417934‑M2/64‑82 (MQ=255)
                              ttttggcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGCATGGAGAAACAGTAGAGAGCTGCGAtac                                                                 >  1:6544943‑M2/64‑93 (MQ=255)
                               ttcggcagcgccaggttggcgcctaatacccggcgcccgagacccgcacgcgtcgccctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaa                                                                <  1:3370620‑M2/32‑1 (MQ=255)
                                  agcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAg                                                             <  1:3577338‑M2/35‑1 (MQ=255)
                                   gcagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGGATGGAGAAACAgtagagagct                                                                        >  2:6879923‑M2/59‑80 (MQ=255)
                                    gagcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCg                                                           <  1:3055308‑M2/37‑1 (MQ=255)
                                     agcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGcga                                                          >  1:7332354‑M2/57‑93 (MQ=255)
                                      gcgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCgtt                                                         >  1:1228830‑M2/56‑93 (MQ=255)
                                       cgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGCAACAGTAGAGAGATGCGa                                                                    >  2:2941326‑M2/55‑82 (MQ=255)
                                       cgccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGa                                                                    <  2:662119‑M2/28‑1 (MQ=255)
                                         ccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATa                                                                  >  2:6106759‑M2/53‑82 (MQ=255)
                                         ccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATa                                                                  >  2:2772459‑M2/53‑82 (MQ=255)
                                         ccaggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAACAGCGTCagc                                                      >  1:2147059‑M2/53‑93 (MQ=255)
                                          caggttggcttctaatccccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaa                                                                 <  2:1215053‑M2/31‑1 (MQ=255)
                                          caggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaa                                                                 >  2:5814289‑M2/52‑82 (MQ=255)
                                          caggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCagga                                                     >  1:315328‑M2/52‑93 (MQ=255)
                                           aggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaa                                                                <  2:2631261‑M2/32‑1 (MQ=255)
                                            ggttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATaaaa                                                               <  2:7272523‑M2/33‑1 (MQ=255)
                                              ttggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAg                                                             <  2:7399391‑M2/35‑1 (MQ=255)
                                              atggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGa                                                 <  1:3470044‑M2/47‑1 (MQ=255)
                                                tgcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATc                                               <  1:735611‑M2/49‑1 (MQ=255)
                                                ggcttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCg                                                           >  2:5791341‑M2/46‑82 (MQ=255)
                                                   ttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCa                                                        >  2:3007496‑M2/43‑82 (MQ=255)
                                                   ttctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCcgg                                            >  1:7229826‑M2/43‑93 (MQ=255)
                                                   gtctaatccccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGc                                            <  1:6264012‑M2/52‑1 (MQ=255)
                                                    tctaatacccggcgcccgagacccgcaggcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAg                                                       <  2:4805302‑M2/41‑1 (MQ=255)
                                                     ctaatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAgg                                                      <  2:5661603‑M2/42‑1 (MQ=255)
                                                      taatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGt                                                     >  2:2564003‑M2/40‑82 (MQ=255)
                                                      taatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGt                                                     >  2:1803391‑M2/40‑82 (MQ=255)
                                                      caatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTaa                                         <  1:4976494‑M2/55‑1 (MQ=255)
                                                       aatccccggcgcacgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                    <  2:2097‑M2/44‑1 (MQ=255)
                                                       aatacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTa                                                    <  2:4232234‑M2/44‑1 (MQ=255)
                                                        atacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATc                                       <  1:4680339‑M2/57‑1 (MQ=255)
                                                        atacccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATc                                       >  1:2811169‑M2/38‑94 (MQ=255)
                                                          acccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGa                                                 <  2:4334133‑M2/47‑1 (MQ=255)
                                                          acccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATcta                                     >  1:2287860‑M2/36‑93 (MQ=255)
                                                           cccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGAt                                                <  2:5292064‑M2/48‑1 (MQ=255)
                                                           cccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGAt                                                >  2:6983234‑M2/35‑82 (MQ=255)
                                                            ccggcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTtag                                   >  1:968593‑M2/34‑93 (MQ=255)
                                                               gcgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGc                                            >  2:7235083‑M2/31‑82 (MQ=255)
                                                                cgcccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGAt                               <  1:2848485‑M2/65‑1 (MQ=255)
                                                                   ccgagacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaa                            <  1:5697653‑M2/68‑1 (MQ=255)
                                                                      tgacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAAt                         <  1:711905‑M2/71‑1 (MQ=255)
                                                                      cgacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAAt                         <  1:5229166‑M2/71‑1 (MQ=255)
                                                                      agacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGAtaaaaaa                         >  1:6345787‑M2/24‑93 (MQ=255)
                                                                       gacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTa                                    <  2:46526‑M2/60‑1 (MQ=255)
                                                                       gacccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAatc                        >  1:6163947‑M2/23‑93 (MQ=255)
                                                                         cccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATg                                  >  2:3357136‑M2/21‑82 (MQ=255)
                                                                          gcgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTa                     <  1:3942913‑M2/75‑1 (MQ=255)
                                                                          ccgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGctt                     >  1:5329578‑M2/20‑93 (MQ=255)
                                                                           cgcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGa                                >  2:1616217‑M2/19‑82 (MQ=255)
                                                                            gcacgcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGAt                               >  2:6442061‑M2/18‑82 (MQ=255)
                                                                                gcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGcac               >  1:92480‑M2/14‑93 (MQ=255)
                                                                                gcgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGcaa               >  1:4057413‑M2/14‑93 (MQ=255)
                                                                                 tgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATa              <  1:2067219‑M2/82‑1 (MQ=255)
                                                                                 cgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaaaa                          <  2:6549228‑M2/70‑1 (MQ=255)
                                                                                 cgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaaaa                          >  2:3701537‑M2/13‑82 (MQ=255)
                                                                                 cgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaaaa                          >  2:6034067‑M2/13‑82 (MQ=255)
                                                                                 cgtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCatt              >  1:1378391‑M2/13‑93 (MQ=255)
                                                                                  gtcgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAg             <  1:6784377‑M2/83‑1 (MQ=255)
                                                                                    ggacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCa           <  1:6038386‑M2/85‑1 (MQ=255)
                                                                                    cgacctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGc                       >  2:2349120‑M2/10‑82 (MQ=255)
                                                                                      tcctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCaaa         <  1:1003‑M2/87‑1 (MQ=255)
                                                                                       tctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAg        <  1:3225475‑M2/88‑1 (MQ=255)
                                                                                        ctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATg                   >  2:5984734‑M2/6‑82 (MQ=255)
                                                                                        ctgcaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATg                   >  2:4995103‑M2/6‑82 (MQ=255)
                                                                                           taGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGtgtg    <  1:4454330‑M2/92‑1 (MQ=255)
                                                                                           gaGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGtgtg    <  1:2002677‑M2/92‑1 (MQ=255)
                                                                                             |                                                                                             
GGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGAC  >  NC_000913/69980‑70166

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: