Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,476,804 T→C T276A (ACG→GCG)  dsdC ← dsd operon activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,476,8040TC98.5% 245.4 / ‑3.8 66T276A (ACG→GCG) dsdCdsd operon activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (33/32);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

TAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTC  >  NC_000913/2476662‑2476946
                                                                                                                                                 |                                                                                                                                              
tAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGc                                                                                                                                        >  2:165239/1‑151 (MQ=255)
         ccAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTTTTCTCTTAGCCAGATAATAAACGCTTCATTTTTTGGCCCCTGCCGGCCCGT‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGgc                                                                                                                                <  2:638069/150‑1 (MQ=255)
            ggACCTCGTGGTTTTTTACATTATAAAGTGTTTTTTCCCTGTTCTCTAAGCCAAATAAAAACCGCTTAATTTTTTGCCCTCTGCCGCCCCGG‑T‑AAGGGGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGgcgc                                                                                                                              <  2:1056937/149‑1 (MQ=255)
              aaCTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGcgacg                                                                                                                          >  1:6876/1‑151 (MQ=255)
                       ttttATTCTTTATGAAGTGTTTTTTCCCTTTCCTCTTGGCCAAAAAAAAAACTCTTCATTTTTGGGCCCCTCCCTCCCCGG‑T‑AATGTGGTG‑ATGAAAAAGTGCTGAGGCCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTc                                                                                                                     <  1:794396/147‑1 (MQ=255)
                       ttttATTCATTATGAGGTGTTTTTTCCCTGTTCTTTTAGCCGAATAATAAACGCTTCAATTTTTGGCCCCTGCCCCCCCGG‑T‑AATGGGGTG‑ATGTAATAGGGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTca                                                                                                                    <  2:702612/148‑1 (MQ=255)
                       ttttATTAATAAGGAGTTGGTTTTTACTTGTTCTCTAAGCCAAAAAATAAACTCTACATTTTTTGGCCTTTTCCT‑CCCGG‑TAAATGTGGTA‑AGGTAATATGGGTGTTGGCTTTCCACCGCCATATCGCCAAACGGCGcgacga                                                                                                                         <  1:563394/143‑1 (MQ=255)
                           aTTCATTATGAAGTGTTTTTTCCCTTTTCTCTTACCCAAAAAAAAAACGCTTCATTTTTTGCCCCCTGCCTGCCCGG‑T‑AATGTGGTG‑ATATAAAAATGCTGATGCCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACt                                                                                                               <  1:959846/149‑1 (MQ=255)
                             tCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGc                                                                                                            >  1:395677/1‑150 (MQ=255)
                              cATTATGAAGTGGTTTTTCCCTGTTCTCTTAGCCATATAATAAACGTTTCAATTTTTGCCCCCTGCCT‑CCCGG‑T‑TATGGGGGGAATGTAAGATGGCTGATGCCATTCCGCCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCa                                                                                                           <  1:85442/150‑1 (MQ=255)
                                ttATGAGGGGTTTTTTCCCTTTTCTCTGAGCGAAATAAAAACCTCTTAATTTTTTGCCCACGGCCCCCCCGG‑T‑AATGTGGTG‑TTGAAAAATTGCTTATGGCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCaa                                                                                                          <  1:1011124/149‑1 (MQ=255)
                                ttATGAATTGTTTTTTCCCGTTTCTTTTAGCAAAATAAAAAACGCTTTATTTTTGGGCCCCTGCCTGCCCGG‑T‑AATGGGGTG‑ATGTAATAATGCTGAGGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAg                                                                                                         <  1:670303/150‑1 (MQ=255)
                                ttAGGAATTGTTTTTTACCTGTTCTCTTAGCAAAATAATAAACCCTTCATTTTTTGCCCCCTCCCTCCCCGG‑T‑AATGGGGTG‑ATTTAATAATGCTGATGGCTTTCCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCa                                                                                                           <  1:1281373/148‑1 (MQ=255)
                                 tATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACCGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGcc                                                                                                       >  2:649841/1‑151 (MQ=255)
                                     aaGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCtttt                                                                                                   >  2:265392/1‑151 (MQ=255)
                                       gtgtTTTTTCCCTTTTCTTTAACCGAAAAAAAAAACTCTTCATTTTTTGCCCCCTGCCCCCCCGT‑T‑AATGGGGTA‑TTTAAAAAATGCTGATGCCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCtttt                                                                                                   <  1:831104/149‑1 (MQ=255)
                                             tttACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTAc                                                                                             >  2:451668/1‑149 (MQ=255)
                                              ttACCTTTTCCTTTAGCAAGAAAAAAAACGATTAAATTTTTGGCCACTGCC‑GCCCGGAT‑ATTGTGTTG‑ATGTATTAATGCGAGGGGCATTGAACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCa                                                                                           <  2:852565/150‑1 (MQ=255)
                                                aCCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGc                                                                                        >  2:779741/1‑151 (MQ=255)
                                                 ccTTTTCTCTTACCCAAAAAAAAACCGCTTCATTTTTTGCCCCCTCCCTCCCCGG‑T‑GATGTGGTG‑ATGTAATAATGCTGATGGCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCg                                                                                       <  2:405893/151‑1 (MQ=255)
                                                 ccTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCg                                                                                       >  2:1006793/1‑151 (MQ=255)
                                                   tGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGtt                                                                                     >  1:223767/1‑151 (MQ=255)
                                                       ctctTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTttct                                                                                 >  1:289975/1‑151 (MQ=255)
                                                         ctTAGCCATAAAAAAAACGTTTTATTTTTGGCCCCCTGCCTGCCCGG‑A‑AATGTGGTG‑ATGTAAAAATTCGGAGGGCTTTTCCCCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTttc                                                                                  <  1:274228/148‑1 (MQ=255)
                                                         ctTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTttcttc                                                                               >  1:786477/1‑151 (MQ=255)
                                                                aaataataAACGTTTAATTTTTTGCCCACTCCCTGCCCAG‑A‑AATGTGTTG‑ATATAATAATGCTAATGCCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCAgcgc                                                                        <  1:272592/149‑1 (MQ=255)
                                                                 gATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGcc                                                                       >  1:285691/1‑151 (MQ=255)
                                                                  ataataAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGcc                                                                       >  1:637197/1‑150 (MQ=255)
                                                                    aataaACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCAcc                                                                    >  1:1066642/1‑151 (MQ=255)
                                                                    aataaACGCTCCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCAcc                                                                    >  1:731404/1‑151 (MQ=255)
                                                                           gCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAatat                                                             >  2:582649/1‑151 (MQ=255)
                                                                              tCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGAt                                                          >  2:1169764/1‑151 (MQ=255)
                                                                                aaTTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTc                                                        >  2:718514/1‑151 (MQ=255)
                                                                                  tttttGGCCACTGCCTCCCGGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCa                                                       <  2:632442/150‑1 (MQ=255)
                                                                                     ttGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATcgc                                                   >  2:326952/1‑151 (MQ=255)
                                                                                     ttGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATcgc                                                   >  2:1242527/1‑151 (MQ=255)
                                                                                         cccctgcccgcCCGG‑T‑AATGTGGTT‑ATGAAATAATGCTGATGGCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCa                                               <  1:126253/148‑1 (MQ=255)
                                                                                         cccccccccccccGG‑A‑GATGTGGTG‑ATATAAAAATCCTGATGGCTTTTCCCCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATcg                                                    <  1:32646/136‑1 (MQ=255)
                                                                                         ccACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGc                                                >  1:1229274/1‑150 (MQ=255)
                                                                                               cctgccCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAAc                                          >  1:767509/1‑150 (MQ=255)
                                                                                               cctgccCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAAc                                          >  1:602923/1‑150 (MQ=255)
                                                                                                ctgccCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAAc                                          >  2:521147/1‑149 (MQ=255)
                                                                                                   cccGG‑T‑AATGTGGTG‑ATTTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGc                                       <  2:354409/149‑1 (MQ=255)
                                                                                                   cccGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCt                                      >  2:789644/1‑150 (MQ=255)
                                                                                                    ccGG‑T‑ATTGGGGTG‑TTGTAATAATGTTGATGGCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTaa                                    <  2:334862/151‑1 (MQ=255)
                                                                                                           aaTGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATc                                 >  1:68013/1‑149 (MQ=255)
                                                                                                             tgtgGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAg                               >  1:29043/1‑149 (MQ=255)
                                                                                                             tgtgGTG‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGa                              >  2:1201792/1‑150 (MQ=255)
                                                                                                             tgtgGTA‑ATGTAATAATGCTGATGGCATTCCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAg                               <  2:1207137/149‑1 (MQ=255)
                                                                                                              gtggtg‑ATGAAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGa                              <  2:887644/149‑1 (MQ=255)
                                                                                                              gtggta‑ATGTAAAAATGCTAATGACATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGaa                             <  1:309444/150‑1 (MQ=255)
                                                                                                                gggg‑ATGAAATAATGCGGAGGGCATGCCTCCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCt                                      <  1:716376/136‑1 (MQ=255)
                                                                                                                 gtg‑ATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGAt                         >  1:236755/1‑151 (MQ=255)
                                                                                                                        taataaTGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGc                   >  1:547958/1‑151 (MQ=255)
                                                                                                                         aataatGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCaaa                     <  2:754834/148‑1 (MQ=255)
                                                                                                                         aaaaatGCTGAGGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGcc                  <  2:1209401/148‑1 (MQ=255)
                                                                                                                          ataatGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGcc                  >  2:1068636/1‑150 (MQ=255)
                                                                                                                            aatGCTGAGGCCTTTCCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCaa                <  2:826695/150‑1 (MQ=255)
                                                                                                                             acgcagATGGAATTCCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGa                          <  2:900128/134‑1 (MQ=255)
                                                                                                                               gCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAAtt              <  2:155441/149‑1 (MQ=255)
                                                                                                                                cTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTcc            >  1:292556/1‑150 (MQ=255)
                                                                                                                                 tGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTcc            >  2:3401/1‑149 (MQ=255)
                                                                                                                                  gATGACTTTCCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCaga         <  1:651216/151‑1 (MQ=255)
                                                                                                                                          ttCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAATCCAATTCCAGAAGATGTc  >  2:1079668/1‑150 (MQ=255)
                                                                                                                                          ttCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTc  <  2:637232/150‑1 (MQ=255)
                                                                                                                                          tgctccGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGt   <  1:682808/145‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGG‑T‑AATGTGGTG‑ATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTC  >  NC_000913/2476662‑2476946

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