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

TAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTC  >  NC_000913/2476662‑2476946
                                                                                                                                                 |                                                                                                                                              
tAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGc                                                                                                                                        >  2:165239/1‑151 (MQ=255)
         ccAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTTTTCTCTTAGCCAGATAATAAACGCTTCATTTTTTGGCCCCTGCCGGCCCGTTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGgc                                                                                                                                <  2:638069/150‑1 (MQ=255)
            ggACCTCGTGGTTTTTTACATTATAAAGTGTTTTTTCCCTGTTCTCTAAGCCAAATAAAAACCGCTTAATTTTTTGCCCTCTGCCGCCCCGGTAAGGGGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGgcgc                                                                                                                              <  2:1056937/149‑1 (MQ=255)
              aaCTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGcgacg                                                                                                                          >  1:6876/1‑151 (MQ=255)
                       ttttATTCTTTATGAAGTGTTTTTTCCCTTTCCTCTTGGCCAAAAAAAAAACTCTTCATTTTTGGGCCCCTCCCTCCCCGGTAATGTGGTGATGAAAAAGTGCTGAGGCCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTc                                                                                                                     <  1:794396/147‑1 (MQ=255)
                       ttttATTCATTATGAGGTGTTTTTTCCCTGTTCTTTTAGCCGAATAATAAACGCTTCAATTTTTGGCCCCTGCCCCCCCGGTAATGGGGTGATGTAATAGGGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTca                                                                                                                    <  2:702612/148‑1 (MQ=255)
                       ttttATTAATAAGGAGTTGGTTTTTACTTGTTCTCTAAGCCAAAAAATAAACTCTACATTTTTTGGCCTTTTCCTCCCGGTAAATGTGGTAAGGTAATATGGGTGTTGGCTTTCCACCGCCATATCGCCAAACGGCGcgacga                                                                                                                         <  1:563394/143‑1 (MQ=255)
                           aTTCATTATGAAGTGTTTTTTCCCTTTTCTCTTACCCAAAAAAAAAACGCTTCATTTTTTGCCCCCTGCCTGCCCGGTAATGTGGTGATATAAAAATGCTGATGCCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACt                                                                                                               <  1:959846/149‑1 (MQ=255)
                             tCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGc                                                                                                            >  1:395677/1‑150 (MQ=255)
                              cATTATGAAGTGGTTTTTCCCTGTTCTCTTAGCCATATAATAAACGTTTCAATTTTTGCCCCCTGCCTCCCGGTTATGGGGGGAATGTAAGATGGCTGATGCCATTCCGCCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCa                                                                                                           <  1:85442/150‑1 (MQ=255)
                                ttATGAGGGGTTTTTTCCCTTTTCTCTGAGCGAAATAAAAACCTCTTAATTTTTTGCCCACGGCCCCCCCGGTAATGTGGTGTTGAAAAATTGCTTATGGCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCaa                                                                                                          <  1:1011124/149‑1 (MQ=255)
                                ttATGAATTGTTTTTTCCCGTTTCTTTTAGCAAAATAAAAAACGCTTTATTTTTGGGCCCCTGCCTGCCCGGTAATGGGGTGATGTAATAATGCTGAGGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAg                                                                                                         <  1:670303/150‑1 (MQ=255)
                                ttAGGAATTGTTTTTTACCTGTTCTCTTAGCAAAATAATAAACCCTTCATTTTTTGCCCCCTCCCTCCCCGGTAATGGGGTGATTTAATAATGCTGATGGCTTTCCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCa                                                                                                           <  1:1281373/148‑1 (MQ=255)
                                 tATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACCGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGcc                                                                                                       >  2:649841/1‑151 (MQ=255)
                                     aaGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCtttt                                                                                                   >  2:265392/1‑151 (MQ=255)
                                       gtgtTTTTTCCCTTTTCTTTAACCGAAAAAAAAAACTCTTCATTTTTTGCCCCCTGCCCCCCCGTTAATGGGGTATTTAAAAAATGCTGATGCCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCtttt                                                                                                   <  1:831104/149‑1 (MQ=255)
                                             tttACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTAc                                                                                             >  2:451668/1‑149 (MQ=255)
                                              ttACCTTTTCCTTTAGCAAGAAAAAAAACGATTAAATTTTTGGCCACTGCCGCCCGGATATTGTGTTGATGTATTAATGCGAGGGGCATTGAACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCa                                                                                           <  2:852565/150‑1 (MQ=255)
                                                aCCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGc                                                                                        >  2:779741/1‑151 (MQ=255)
                                                 ccTTTTCTCTTACCCAAAAAAAAACCGCTTCATTTTTTGCCCCCTCCCTCCCCGGTGATGTGGTGATGTAATAATGCTGATGGCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCg                                                                                       <  2:405893/151‑1 (MQ=255)
                                                 ccTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCg                                                                                       >  2:1006793/1‑151 (MQ=255)
                                                   tGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGtt                                                                                     >  1:223767/1‑151 (MQ=255)
                                                       ctctTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTttct                                                                                 >  1:289975/1‑151 (MQ=255)
                                                         ctTAGCCATAAAAAAAACGTTTTATTTTTGGCCCCCTGCCTGCCCGGAAATGTGGTGATGTAAAAATTCGGAGGGCTTTTCCCCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTttc                                                                                  <  1:274228/148‑1 (MQ=255)
                                                         ctTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTttcttc                                                                               >  1:786477/1‑151 (MQ=255)
                                                                aaataataAACGTTTAATTTTTTGCCCACTCCCTGCCCAGAAATGTGTTGATATAATAATGCTAATGCCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCAgcgc                                                                        <  1:272592/149‑1 (MQ=255)
                                                                 gATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGcc                                                                       >  1:285691/1‑151 (MQ=255)
                                                                  ataataAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGcc                                                                       >  1:637197/1‑150 (MQ=255)
                                                                    aataaACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCAcc                                                                    >  1:1066642/1‑151 (MQ=255)
                                                                    aataaACGCTCCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCAcc                                                                    >  1:731404/1‑151 (MQ=255)
                                                                           gCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAatat                                                             >  2:582649/1‑151 (MQ=255)
                                                                              tCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGAt                                                          >  2:1169764/1‑151 (MQ=255)
                                                                                aaTTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTc                                                        >  2:718514/1‑151 (MQ=255)
                                                                                  tttttGGCCACTGCCTCCCGGGTAATGTGGTGATGTAATAATGCTGATGGCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCa                                                       <  2:632442/150‑1 (MQ=255)
                                                                                     ttGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATcgc                                                   >  2:326952/1‑151 (MQ=255)
                                                                                     ttGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATcgc                                                   >  2:1242527/1‑151 (MQ=255)
                                                                                         cccctgcccgcCCGGTAATGTGGTTATGAAATAATGCTGATGGCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCa                                               <  1:126253/148‑1 (MQ=255)
                                                                                         cccccccccccccGGAGATGTGGTGATATAAAAATCCTGATGGCTTTTCCCCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATcg                                                    <  1:32646/136‑1 (MQ=255)
                                                                                         ccACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGc                                                >  1:1229274/1‑150 (MQ=255)
                                                                                               cctgccCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAAc                                          >  1:767509/1‑150 (MQ=255)
                                                                                               cctgccCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAAc                                          >  1:602923/1‑150 (MQ=255)
                                                                                                ctgccCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAAc                                          >  2:521147/1‑149 (MQ=255)
                                                                                                   cccGGTAATGTGGTGATTTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGc                                       <  2:354409/149‑1 (MQ=255)
                                                                                                   cccGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCt                                      >  2:789644/1‑150 (MQ=255)
                                                                                                    ccGGTATTGGGGTGTTGTAATAATGTTGATGGCTTTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTaa                                    <  2:334862/151‑1 (MQ=255)
                                                                                                           aaTGTGGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATc                                 >  1:68013/1‑149 (MQ=255)
                                                                                                             tgtgGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAg                               >  1:29043/1‑149 (MQ=255)
                                                                                                             tgtgGTGATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGa                              >  2:1201792/1‑150 (MQ=255)
                                                                                                             tgtgGTAATGTAATAATGCTGATGGCATTCCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAg                               <  2:1207137/149‑1 (MQ=255)
                                                                                                              gtggtgATGAAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGa                              <  2:887644/149‑1 (MQ=255)
                                                                                                              gtggtaATGTAAAAATGCTAATGACATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGaa                             <  1:309444/150‑1 (MQ=255)
                                                                                                                ggggATGAAATAATGCGGAGGGCATGCCTCCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCt                                      <  1:716376/136‑1 (MQ=255)
                                                                                                                 gtgATGTAATAATGCTGATGGCATTTCACCGCCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGAt                         >  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)
                                                                                                                                                 |                                                                                                                                              
TAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATCGCCACCCCAATATGATTCATCGCGGCAATAACAGCTAAATCAGAACGATCAAAGCCAATTCCAGAAGATGTC  >  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: