Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,661,482 G→C 7.9% P17R (CCG→CGG)  iscS ← cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,661,4820GC 7.9% 378.0 / 18.6 124P17R (CCG→CGG) iscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (55/58);  new base C (7/3);  total (63/61)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

Only 100 of 129 total aligned reads displayed.

GCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTG  >  NC_000913/2661333‑2661631
                                                                                                                                                      |                                                                                                                                                     
 cgGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGACTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGTTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCgg                                                                                                                                                      >  2:273859/1‑150 (MQ=255)
   gATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGt                                                                                                                                                    >  2:1101996/1‑150 (MQ=255)
   gATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGt                                                                                                                                                    >  1:367329/1‑150 (MQ=255)
   gATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCATCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCgtgt                                                                                                                                                    >  2:1456331/1‑150 (MQ=255)
       aGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCAc                                                                                                                                                >  2:1233496/1‑150 (MQ=255)
        gcgcCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCAcc                                                                                                                                               >  1:1292005/1‑150 (MQ=255)
              aCCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCTTAAAGTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTg                                                                                                                                         >  2:882191/1‑150 (MQ=255)
                cAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGt                                                                                                                                       <  2:1409644/150‑1 (MQ=255)
                cAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGt                                                                                                                                       <  2:421105/150‑1 (MQ=255)
                 aGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGATACGGTGAGAACGGGAGGCCGGGTTACAAAAGGTTCCCTCCATCTTCATAAACTTCATCATTTTCTCGGCAA‑CTCGCGGGTCCGACGGCGTgggt                                                                                                                                      >  2:1172123/1‑148 (MQ=255)
                  gATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTg                                                                                                                                     >  2:1466253/1‑150 (MQ=255)
                   aTCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGTgttgc                                                                                                                                     <  2:1400046/150‑5 (MQ=255)
                     cGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                  <  2:92657/150‑1 (MQ=255)
                     cGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                  <  1:1533748/150‑1 (MQ=255)
                       gCAATCTGATTACGGGCGAAATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                                <  1:1619978/150‑1 (MQ=255)
                        cAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAACCTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGT‑GTTGCCGaga                                                                                                                               <  1:1020908/150‑2 (MQ=37)
                        cAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGagt                                                                                                                               >  1:1565013/1‑150 (MQ=255)
                          aTCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGT‑GTTGCCGagaaa                                                                                                                             <  2:217699/150‑4 (MQ=37)
                             tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCTGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                          >  2:1293636/1‑150 (MQ=255)
                             tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                          <  2:1621682/150‑1 (MQ=255)
                             tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCATCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                          >  1:711684/1‑150 (MQ=255)
                             tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTAACAAAGGTTCCCTCCATCATCATAAACTGCATCATTTTCTCTGAAA‑CACGAGAGTCCACCCGCGTGGTAGCGGagaactcc                                                                                                                          >  1:1605251/1‑144 (MQ=255)
                                    gggCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCTTCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGAGGAGTAGTCGAGATaaa                                                                                                                   >  2:911629/1‑150 (MQ=255)
                                     ggCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGAATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                                                                                                  >  2:1229779/1‑150 (MQ=255)
                                       cGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCg                                                                                                                >  2:784111/1‑150 (MQ=255)
                                          tatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTa                                                                                                             <  2:367329/150‑1 (MQ=255)
                                          tatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTa                                                                                                             <  1:202959/150‑1 (MQ=255)
                                             cTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGTTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAtt                                                                                                          <  1:556775/150‑1 (MQ=255)
                                             cTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAtt                                                                                                          >  2:671283/1‑150 (MQ=255)
                                              tACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                                                                                                         <  2:887917/150‑1 (MQ=255)
                                                ccGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                                                                                                       <  1:1005548/150‑1 (MQ=255)
                                                 cGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAt                                                                                                      >  1:728759/1‑150 (MQ=255)
                                                  gCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAtt                                                                                                     >  1:183068/1‑150 (MQ=255)
                                                    ttcttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGc                                                                                                   >  2:1235933/1‑150 (MQ=255)
                                                     tcttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGct                                                                                                  >  2:1397222/1‑150 (MQ=255)
                                                      cttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                                 >  2:779650/1‑150 (MQ=255)
                                                       ttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctct                                                                                                >  1:1687246/1‑150 (MQ=255)
                                                       ttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctct                                                                                                >  1:574816/1‑150 (MQ=255)
                                                              tGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                         >  2:948210/1‑150 (MQ=255)
                                                              tGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                         >  2:559832/1‑150 (MQ=255)
                                                              tGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                         <  1:1635608/150‑1 (MQ=255)
                                                                 cAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCg                                                                                      <  1:1141730/150‑1 (MQ=255)
                                                                 cAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCg                                                                                      <  1:531395/150‑1 (MQ=255)
                                                                  aGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGt                                                                                     <  1:1460552/150‑1 (MQ=255)
                                                                    ccGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTAc                                                                                   <  2:286635/150‑1 (MQ=255)
                                                                         aCGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCACGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCAc                                                                              <  1:1320544/150‑1 (MQ=255)
                                                                          cGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACt                                                                             <  2:484267/150‑1 (MQ=255)
                                                                           ggTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTc                                                                            <  1:1688259/150‑1 (MQ=255)
                                                                           ggTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGGTGCCTAGTAGTCGAGATAAATCGGTAATTTCATTTCTCTATAAACTCCGGTCATCACta                                                                            >  2:869879/1‑149 (MQ=255)
                                                                             tGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATATACTCCGTACATCACTCaa                                                                          >  2:73974/1‑150 (MQ=255)
                                                                              gagaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAAt                                                                         >  1:398848/1‑150 (MQ=255)
                                                                                taaCGGGAGTCAGGGTTACCAAAGGGTCCGTCCTTCGTCATAAACTGCTTCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGc                                                                       <  1:124623/149‑1 (MQ=255)
                                                                                gaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCCCCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGc                                                                       >  1:31464/1‑150 (MQ=255)
                                                                                   cGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAg                                                                    <  1:638316/150‑1 (MQ=255)
                                                                                     ggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGa                                                                  >  2:747912/1‑150 (MQ=255)
                                                                                        ggCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCCGGTCCACCGGCGTGGTTGGGCAGGAGTCGAGATAAATAGGTAATTTCATTGCTCTATTAAATACGTACATCACTCAATTCAAGGAAta                                                               >  1:1557435/1‑149 (MQ=255)
                                                                                           cGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCTTCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAgg                                                            <  1:7115/150‑1 (MQ=255)
                                                                                            gggTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGc                                                           <  1:609826/150‑1 (MQ=255)
                                                                                            gggTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCCATAAACTCCGTACATCACTCAATGCAAGGAATCAGGc                                                           >  1:796461/1‑150 (MQ=255)
                                                                                            gggTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAACGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGc                                                           >  1:1325472/1‑150 (MQ=255)
                                                                                                tACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTAcc                                                       <  2:145326/150‑1 (MQ=255)
                                                                                                 aCCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCg                                                      <  1:779650/150‑1 (MQ=255)
                                                                                                 aCCAAAGGTTCCGTCCATCGTCATAAACTGCATAATTATCTCGGCAC‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTAATTCCTCTATAAACTCCGTCCATCACTCAATGCAACGAATCAGGCTTCCg                                                      >  1:287352/1‑150 (MQ=255)
                                                                                                  ccgaaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTAATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGATACCgg                                                     >  2:877167/4‑150 (MQ=255)
                                                                                                  ccAAAGGTTCCGTCCATCGTCATAAACTGCATCTTTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCgg                                                     >  2:570564/1‑150 (MQ=255)
                                                                                                    aaaGGTTCCGTCCATCGTCATATACTGGATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATATATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCt                                                   <  2:49400/150‑1 (MQ=255)
                                                                                                    aaaGGTTCCGTCCATCGTCATAATCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCt                                                   <  1:1077830/150‑1 (MQ=255)
                                                                                                     aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTg                                                  <  1:1335271/150‑1 (MQ=255)
                                                                                                       ggTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGa                                                <  1:671283/150‑1 (MQ=255)
                                                                                                        gTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGAt                                               >  2:1610655/1‑150 (MQ=255)
                                                                                                         ttCCGTCCATGGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATg                                              <  1:273859/150‑1 (MQ=255)
                                                                                                         ttCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATg                                              <  2:1292005/150‑1 (MQ=255)
                                                                                                             gTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACg                                          <  2:398848/150‑1 (MQ=255)
                                                                                                               ccATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGAc                                        >  1:1662125/1‑150 (MQ=255)
                                                                                                                cATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATAACTCAATGCAAGGAATCAGGCTACCGGCTGCATGTACGccc                                       >  2:317613/1‑150 (MQ=255)
                                                                                                                 aTCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCg                                      >  1:1450239/1‑150 (MQ=255)
                                                                                                                  tCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgt                                     >  2:472000/1‑150 (MQ=255)
                                                                                                                   cGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgtg                                    >  1:1641111/1‑150 (MQ=255)
                                                                                                                      nntNAACTGNATCATTTTCTCTGCAA‑CACGAGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGttt                                 <  2:1132640/148‑1 (MQ=255)
                                                                                                                         aaaCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACg                              <  2:1212106/150‑1 (MQ=255)
                                                                                                                            cTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                           <  1:913889/150‑1 (MQ=255)
                                                                                                                            cTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                           <  1:1101996/150‑1 (MQ=255)
                                                                                                                            atGCAGCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                           <  1:1456331/149‑1 (MQ=255)
                                                                                                                              gCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTa                         <  1:531666/150‑1 (MQ=255)
                                                                                                                                 tcatTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAttt                      <  2:523738/150‑1 (MQ=255)
                                                                                                                                   atTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTAAg                    >  1:1281161/1‑150 (MQ=255)
                                                                                                                                    ttttCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGc                   <  1:1233496/150‑1 (MQ=255)
                                                                                                                                     ctactcCGCAACCACGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTACATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGTGTATTTAGc                   >  2:766386/4‑150 (MQ=255)
                                                                                                                                      ttCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCAc                 <  1:1046802/150‑1 (MQ=255)
                                                                                                                                       actcCGCAACCACGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCAc                 >  2:410904/2‑150 (MQ=255)
                                                                                                                                        ctcCGCAACCACGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCCGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACt                >  2:710530/1‑150 (MQ=255)
                                                                                                                                         tcGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTcc              <  1:882191/150‑1 (MQ=255)
                                                                                                                                          ccgcaaccacGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTcc              >  1:1048643/8‑150 (MQ=255)
                                                                                                                                           ggCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCgg            >  2:535046/1‑150 (MQ=255)
                                                                                                                                            gcaaccacGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCACGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCgg            >  1:1488322/6‑150 (MQ=255)
                                                                                                                                             cAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGcc          >  2:1405489/1‑150 (MQ=255)
                                                                                                                                              aa‑CCCGGGGGTCCACCGGGGTTGTTGCGGTGTTGTGGAGCTAAATAGGTAATTTCAATGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGGTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCAGGCCt         <  2:592289/150‑1 (MQ=255)
                                                                                                                                              aa‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCt         <  2:1298775/150‑1 (MQ=255)
                                                                                                                                                      gggTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACTACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTg  <  2:818367/150‑1 (MQ=255)
                                                                                                                                                      gggTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTg  >  2:682127/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCGGATCAGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTG  >  NC_000913/2661333‑2661631

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