Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,661,486 C→G 6.9% D16H (GAC→CAC)  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,4860CG 6.9% 423.8 / 15.0 127D16H (GAC→CAC) iscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (57/61);  new base G (6/3);  total (63/64)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

AGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATT  >  NC_000913/2661340‑2661635
                                                                                                                                                   |                                                                                                                                                     
aGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCAc                                                                                                                                                    >  2:1233496/1‑150 (MQ=255)
aGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCAc                                                                                                                                                    >  2:531666/1‑150 (MQ=255)
 gcgcCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCAcc                                                                                                                                                   >  1:1292005/1‑150 (MQ=255)
       aCCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCTTAAAGTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTg                                                                                                                                             >  2:882191/1‑150 (MQ=255)
         cAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGt                                                                                                                                           <  2:1409644/150‑1 (MQ=255)
         cAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGt                                                                                                                                           <  2:421105/150‑1 (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)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCATCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                              >  1:711684/1‑150 (MQ=255)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTAACAAAGGTTCCCTCCATCATCATAAACTGCATCATTTTCTCTGAAACACGAGAGTCCACCCGCGTGGTAGCGGagaactcc                                                                                                                              >  1:1605251/1‑144 (MQ=255)
                        aTTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCgag                                                                                                                            <  1:524741/150‑1 (MQ=255)
                             gggCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCTTCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGAGGAGTAGTCGAGATaaa                                                                                                                       >  2:911629/1‑150 (MQ=255)
                              ggCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGAATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                                                                                                      >  2:1229779/1‑150 (MQ=255)
                                cGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCg                                                                                                                    >  2:784111/1‑150 (MQ=255)
                                   tatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTa                                                                                                                 <  2:367329/150‑1 (MQ=255)
                                   tatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTa                                                                                                                 <  1:202959/150‑1 (MQ=255)
                                   tatCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTa                                                                                                                 <  1:1521103/150‑1 (MQ=255)
                                      cTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGTTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAtt                                                                                                              <  1:556775/150‑1 (MQ=255)
                                      cTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAtt                                                                                                              >  2:671283/1‑150 (MQ=255)
                                       tACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAAttt                                                                                                             <  2:887917/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)
                                                       tGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                             >  2:948210/1‑150 (MQ=255)
                                                       tGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                             >  2:559832/1‑150 (MQ=255)
                                                       tGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                             <  1:1635608/150‑1 (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)
                                                                  aCGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCACGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCAc                                                                                  <  1:1320544/150‑1 (MQ=255)
                                                                   cGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACt                                                                                 <  2:484267/150‑1 (MQ=255)
                                                                    ggTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTc                                                                                <  1:1688259/150‑1 (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)
                                                                              ggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGa                                                                      >  2:747912/1‑150 (MQ=255)
                                                                                 ggCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCCGGTCCACCGGCGTGGTTGGGCAGGAGTCGAGATAAATAGGTAATTTCATTGCTCTATTAAATACGTACATCACTCAATTCAAGGAAta                                                                   >  1:1557435/1‑149 (MQ=255)
                                                                                     gggTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAACGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGc                                                               >  1:1325472/1‑150 (MQ=255)
                                                                                         tACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTAcc                                                           <  2:145326/150‑1 (MQ=255)
                                                                                          aCCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCg                                                          <  2:396143/150‑1 (MQ=255)
                                                                                          aCCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCg                                                          <  1:779650/150‑1 (MQ=255)
                                                                                          aCCAAAGGTTCCGTCCATCGTCATAAACTGCATAATTATCTCGGCACCACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTAATTCCTCTATAAACTCCGTCCATCACTCAATGCAACGAATCAGGCTTCCg                                                          >  1:287352/1‑150 (MQ=255)
                                                                                           ccgaaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTAATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGATACCgg                                                         >  2:877167/4‑150 (MQ=255)
                                                                                           ccAAAGGTTCCGTCCATCGTCATAAACTGCATCTTTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCgg                                                         >  2:570564/1‑150 (MQ=255)
                                                                                             aaaGGTTCCGTCCATCGTCATATACTGGATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATATATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCt                                                       <  2:49400/150‑1 (MQ=255)
                                                                                             aaaGGTTCCGTCCATCGTCATAATCTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCt                                                       <  1:1077830/150‑1 (MQ=255)
                                                                                              aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTg                                                      <  1:1335271/150‑1 (MQ=255)
                                                                                                ggTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGa                                                    <  1:671283/150‑1 (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)
                                                                                                               nntNAACTGNATCATTTTCTCTGCAACACGAGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGttt                                     <  2:1132640/148‑1 (MQ=255)
                                                                                                                  aaaCTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACg                                  <  2:1212106/150‑1 (MQ=255)
                                                                                                                     cTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                               <  1:913889/150‑1 (MQ=255)
                                                                                                                     cTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                               <  1:1101996/150‑1 (MQ=255)
                                                                                                                     atGCAGCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                               <  1:1456331/149‑1 (MQ=255)
                                                                                                                       gCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTa                             <  1:531666/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)
                                                                                                                                 ctcGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCCCTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTc                   >  1:572865/1‑150 (MQ=255)
                                                                                                                                  tcGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTcc                  <  1:882191/150‑1 (MQ=255)
                                                                                                                                   ccgcaaccacGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTcc                  >  1:1048643/8‑150 (MQ=255)
                                                                                                                                    ggCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCgg                >  2:535046/1‑150 (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:1298775/150‑1 (MQ=255)
                                                                                                                                           acGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTTCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGAt          <  2:1565013/150‑1 (MQ=255)
                                                                                                                                               gggTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACTACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTg      <  2:818367/150‑1 (MQ=255)
                                                                                                                                               gggTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTg      >  2:682127/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)
                                                                                                                                                   |                                                                                                                                                     
AGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATT  >  NC_000913/2661340‑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: