Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,661,491 G→C 7.9% P14R (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,4910GC 7.9% 382.6 / 17.6 125P14R (CCG→CGG) iscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (54/61);  new base C (7/3);  total (61/64)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 125 total aligned reads displayed.

ACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATT  >  NC_000913/2661347‑2661635
                                                                                                                                                 |                                                                                                                                                
aCCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCTTAAAGTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTg                                                                                                                                             >  2:882191/1‑150 (MQ=255)
   aGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGATACGGTGAGAACGGGAGGCCGGGTTACAAAAGGTTCCCTCCATCTTCATAAACTTCATCATTTTCTCGGCAACTCGCGGGTCCGACGGCGTgggt                                                                                                                                          >  2:1172123/1‑148 (MQ=255)
    gATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTg                                                                                                                                         >  2:1466253/1‑150 (MQ=255)
     aTCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGTgttgc                                                                                                                                         <  2:1400046/150‑5 (MQ=255)
       cGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCTACACGCGGGTCCACCGGCGTGGTAGCgg                                                                                                                                      >  2:219345/1‑150 (MQ=255)
       cGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                      <  2:92657/150‑1 (MQ=255)
       cGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                      <  1:1533748/150‑1 (MQ=255)
         gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                                    <  1:702193/150‑1 (MQ=255)
         gCAATCTGATTACGGGCGAAATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                                    <  1:1619978/150‑1 (MQ=255)
          cAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAACCTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGTGTTGCCGaga                                                                                                                                   <  1:1020908/150‑2 (MQ=37)
          cAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGagt                                                                                                                                   >  1:1565013/1‑150 (MQ=255)
            aTCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGTGTTGCCGagaaa                                                                                                                                 <  2:217699/150‑4 (MQ=37)
               tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCTGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                              >  2:1293636/1‑150 (MQ=255)
               tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                              <  2:1621682/150‑1 (MQ=255)
               tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTAACAAAGGTTCCCTCCATCATCATAAACTGCATCATTTTCTCTGAAACACGAGAGTCCACCCGCGTGGTAGCGGagaactcc                                                                                                                              >  1:1605251/1‑144 (MQ=255)
                 aTTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                            <  1:524741/150‑1 (MQ=255)
                 aTTACGGGCGATAGCTACCGCTTCTTCAGCCTCCCAGCCGAAACGGTGATAACGGGAGGCCGGGTTACCAACGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                            <  2:1196288/150‑1 (MQ=255)
                      gggCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCTTCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGAGGAGTAGTCGAGATaaa                                                                                                                       >  2:911629/1‑150 (MQ=255)
                         cGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCg                                                                                                                    >  2:784111/1‑150 (MQ=255)
                            tatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTa                                                                                                                 <  2:367329/150‑1 (MQ=255)
                            tatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTa                                                                                                                 <  1:1521103/150‑1 (MQ=255)
                            tatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTa                                                                                                                 <  1:202959/150‑1 (MQ=255)
                             atCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTaa                                                                                                                >  1:12982/1‑150 (MQ=255)
                               cTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGTTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAtt                                                                                                              <  1:556775/150‑1 (MQ=255)
                                  ccGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                                                                                                           <  1:1005548/150‑1 (MQ=255)
                                   cGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAt                                                                                                          >  1:728759/1‑150 (MQ=255)
                                    gCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAtt                                                                                                         >  1:183068/1‑150 (MQ=255)
                                      ttcttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGc                                                                                                       >  2:1235933/1‑150 (MQ=255)
                                       tcttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGct                                                                                                      >  2:1397222/1‑150 (MQ=255)
                                        cttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                                     >  2:779650/1‑150 (MQ=255)
                                         ttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctct                                                                                                    >  1:1687246/1‑150 (MQ=255)
                                         ttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctct                                                                                                    >  1:574816/1‑150 (MQ=255)
                                           cAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCCTCGTCATAAACTGCAGCAGTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCCTTGCTCtat                                                                                                  >  1:1210244/1‑150 (MQ=255)
                                            aGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCtata                                                                                                 <  2:586952/150‑1 (MQ=255)
                                                tGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                             >  2:948210/1‑150 (MQ=255)
                                                tGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                             >  2:559832/1‑150 (MQ=255)
                                                   cAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCg                                                                                          >  2:1307214/1‑150 (MQ=255)
                                                   cAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCg                                                                                          <  1:1141730/150‑1 (MQ=255)
                                                    aGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGt                                                                                         <  1:1460552/150‑1 (MQ=255)
                                                     cccGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTa                                                                                        <  1:1648534/149‑1 (MQ=255)
                                                      ccGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTAc                                                                                       <  2:286635/150‑1 (MQ=255)
                                                           aCGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCAc                                                                                  <  2:1415660/150‑1 (MQ=255)
                                                           aCGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCACGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCAc                                                                                  <  1:1320544/150‑1 (MQ=255)
                                                            cGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACt                                                                                 <  2:484267/150‑1 (MQ=255)
                                                             ggTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGGTGCCTAGTAGTCGAGATAAATCGGTAATTTCATTTCTCTATAAACTCCGGTCATCACta                                                                                >  2:869879/1‑149 (MQ=255)
                                                               tGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATATACTCCGTACATCACTCaa                                                                              >  2:73974/1‑150 (MQ=255)
                                                                gagaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAAt                                                                             >  1:398848/1‑150 (MQ=255)
                                                                  taaCGGGAGTCAGGGTTACCAAAGGGTCCGTCCTTCGTCATAAACTGCTTCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGc                                                                           <  1:124623/149‑1 (MQ=255)
                                                                  gaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCCCCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGc                                                                           >  1:31464/1‑150 (MQ=255)
                                                                     cGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAg                                                                        <  1:638316/150‑1 (MQ=255)
                                                                             cGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCTTCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAgg                                                                <  1:7115/150‑1 (MQ=255)
                                                                              gggTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCCATAAACTCCGTACATCACTCAATGCAAGGAATCAGGc                                                               >  1:796461/1‑150 (MQ=255)
                                                                                  tACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTAcc                                                           <  2:145326/150‑1 (MQ=255)
                                                                                   aCCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCg                                                          <  1:779650/150‑1 (MQ=255)
                                                                                   aCCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCg                                                          <  2:396143/150‑1 (MQ=255)
                                                                                    ccgaaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTAATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGATACCgg                                                         >  2:877167/4‑150 (MQ=255)
                                                                                    ccAAAGGTTCCGTCCATCGTCATAAACTGCATCTTTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCgg                                                         >  2:570564/1‑150 (MQ=255)
                                                                                    ccAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCgg                                                         >  1:1598977/1‑150 (MQ=255)
                                                                                     cAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATGACTCAATGCAAGGAATCAGGCTACCGGc                                                        >  1:1416657/1‑150 (MQ=255)
                                                                                      aaaGGTTCCGTCCATCGTCATATACTGGATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATATATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCt                                                       <  2:49400/150‑1 (MQ=255)
                                                                                      aaaGGTTCCGTCCATCGTCATAATCTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCt                                                       <  1:1077830/150‑1 (MQ=255)
                                                                                      aaaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCt                                                       >  2:776134/1‑150 (MQ=255)
                                                                                       aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTg                                                      <  1:1335271/150‑1 (MQ=255)
                                                                                         ggTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGa                                                    <  1:671283/150‑1 (MQ=255)
                                                                                          gTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGAt                                                   >  2:1291548/1‑150 (MQ=255)
                                                                                           ttCCGTCCATGGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATg                                                  <  1:273859/150‑1 (MQ=255)
                                                                                           ttCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATg                                                  <  2:1292005/150‑1 (MQ=255)
                                                                                             ccGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTa                                                <  2:728759/150‑1 (MQ=255)
                                                                                             ccGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTa                                                <  1:770483/150‑1 (MQ=255)
                                                                                               gTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACg                                              <  2:398848/150‑1 (MQ=255)
                                                                                                 ccATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGAc                                            >  1:1662125/1‑150 (MQ=255)
                                                                                                   aTCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCg                                          >  1:1450239/1‑150 (MQ=255)
                                                                                                    tCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgt                                         >  2:472000/1‑150 (MQ=255)
                                                                                                     cGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgtg                                        >  1:1641111/1‑150 (MQ=255)
                                                                                                           aaaCTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACg                                  <  2:1212106/150‑1 (MQ=255)
                                                                                                              cTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                               <  1:1101996/150‑1 (MQ=255)
                                                                                                              cTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                               <  1:913889/150‑1 (MQ=255)
                                                                                                                gCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTa                             <  1:531666/150‑1 (MQ=255)
                                                                                                                   tcatTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAttt                          <  2:523738/150‑1 (MQ=255)
                                                                                                                     atTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTAAg                        >  1:1281161/1‑150 (MQ=255)
                                                                                                                      ttttCTCGGCAACACGCGGGGCCACCGGCGTGGTTGCGGAGTAGTCGAGATTAATCGGTAATTTCATTGCTCTCTATACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGc                       <  1:92533/150‑1 (MQ=255)
                                                                                                                       ctactcCGCAACCACGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTACATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGTGTATTTAGc                       >  2:766386/4‑150 (MQ=255)
                                                                                                                        ttCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCAc                     <  1:1046802/150‑1 (MQ=255)
                                                                                                                         actcCGCAACCACGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCAc                     >  2:410904/2‑150 (MQ=255)
                                                                                                                          ctcGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCCCTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTc                   >  1:572865/1‑150 (MQ=255)
                                                                                                                          ctcCGCAACCACGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCCGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACt                    >  2:710530/1‑150 (MQ=255)
                                                                                                                            ccgcaaccacGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTcc                  >  1:1048643/8‑150 (MQ=255)
                                                                                                                             ggCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCgg                <  1:1720765/150‑1 (MQ=255)
                                                                                                                              gcaaccacGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCACGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCgg                >  1:1488322/6‑150 (MQ=255)
                                                                                                                              gcaaccacGCCGGTGGACCCGCGTGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCATGGAATCAGGCTACCGGCTGGATGTACGCCCGTGTTTACGGAGTATTTAGCACTCCgtg               >  1:289312/6‑148 (MQ=255)
                                                                                                                               cAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATGGGTAATTTCATTGCTCTATAAAGTCCGTACATCCCTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGCCCGTGTTTACGGAGTATTTAGCACTCCGGcc              >  2:872747/1‑150 (MQ=255)
                                                                                                                               cAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGcc              >  2:1405489/1‑150 (MQ=255)
                                                                                                                                aaCCCGGGGGTCCACCGGGGTTGTTGCGGTGTTGTGGAGCTAAATAGGTAATTTCAATGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGGTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCAGGCCt             <  2:592289/150‑1 (MQ=255)
                                                                                                                                aaCACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCt             >  2:909176/1‑150 (MQ=255)
                                                                                                                                    acGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTTCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGAt          <  2:1565013/150‑1 (MQ=255)
                                                                                                                                        gggTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTg      >  2:682127/1‑150 (MQ=255)
                                                                                                                                         ggTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCCGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGa     >  2:1041150/1‑150 (MQ=255)
                                                                                                                                          gTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGaa    <  1:1397222/150‑1 (MQ=255)
                                                                                                                                           tCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAAt   <  2:1075408/150‑1 (MQ=255)
                                                                                                                                            ccACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAAtt  >  2:1190647/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATT  >  NC_000913/2661347‑2661635

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: