Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,661,487 C→G 5.2% V15V (GTG→GTC)  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,4870CG 5.2% 418.3 / 8.1 137V15V (GTG→GTC) iscScysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (59/69);  new base G (4/3);  total (64/73)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 139 total aligned reads displayed.

AGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTC  >  NC_000913/2661340‑2661636
                                                                                                                                                    |                                                                                                                                                     
aGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCAc                                                                                                                                                     >  1:1887840/1‑150 (MQ=255)
 gcgcCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCAcc                                                                                                                                                    <  1:620102/150‑1 (MQ=255)
       aCCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTg                                                                                                                                              <  1:1526695/150‑1 (MQ=255)
        ccAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTgg                                                                                                                                             >  1:119316/1‑150 (MQ=255)
        ccAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTgg                                                                                                                                             <  1:1516949/150‑1 (MQ=255)
           taTCCGCAATCTGATTACGGTAGATATCTCCCGCTTCTGCAGCCTGCCAGCCGAAACGGTGAGATCGGGAGGCCGGGTTACCACAGGTTCAGGCCATCGTCATAATCTGCATCATTTTCTCGGCAA‑CACGCGGGGCCCCCGTCGTGGTTg                                                                                                                                          <  1:903744/149‑1 (MQ=255)
           gATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTg                                                                                                                                          >  1:1394919/1‑150 (MQ=255)
           gATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCAACCGGCGGGGTTg                                                                                                                                          >  2:1234556/1‑150 (MQ=255)
             tCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGTgttgcc                                                                                                                                         <  1:1743277/150‑6 (MQ=255)
              ccgCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCATGGCAGCCGAAACGGTGAGAACCGGAGGCCGGGTTACCAAACGTTCCGTTCATCGTCATCAACTGCATTATTTTCTCGGCAA‑CACGGGGGTCCACCGGCGTGGTTGCgg                                                                                                                                       <  2:250580/148‑1 (MQ=255)
              cGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                       <  2:722770/150‑1 (MQ=255)
              cGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                       <  2:1843839/150‑1 (MQ=255)
              cGGCAATCTGATCACGGGCGATAGCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGGGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                       <  1:1203597/150‑1 (MQ=255)
                gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                                     <  1:778722/150‑1 (MQ=255)
                gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                                     <  1:1118945/150‑1 (MQ=255)
                gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGT‑GTTGCCGag                                                                                                                                     <  2:446213/150‑1 (MQ=38)
                 cAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGagt                                                                                                                                    >  2:645162/1‑150 (MQ=255)
                 cAATCTGATTACGGGCGATATCTACCGCTTCTTAAGCCTGCCAGCCGAAACTGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGagt                                                                                                                                    <  2:1923958/150‑1 (MQ=255)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGTAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGT‑GTTGCCGagaaaatg                                                                                                                               <  1:1386673/150‑7 (MQ=37)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCGACCGGCGGGGTTGCGGAGGAGTCg                                                                                                                               >  2:1981826/1‑150 (MQ=255)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                               <  2:1689455/150‑1 (MQ=255)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                               >  2:759532/1‑150 (MQ=255)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                               >  2:1610264/1‑150 (MQ=255)
                       gATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCCGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCga                                                                                                                              <  1:1893438/150‑1 (MQ=255)
                            cGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGCAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATaa                                                                                                                         >  2:1454340/1‑150 (MQ=255)
                              ggCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                                                                                                       >  1:456428/1‑150 (MQ=255)
                               gCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATc                                                                                                                      >  2:660753/1‑150 (MQ=255)
                                      cTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGGGGTTGCGGAGTAGTCGAGATAACTCGGTAAtt                                                                                                               >  1:520981/1‑150 (MQ=255)
                                       tACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATTAATCGGTAAttt                                                                                                              >  1:1759665/1‑150 (MQ=255)
                                         ccGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                                                                                                            <  2:116924/150‑1 (MQ=255)
                                          cGCTGCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAt                                                                                                           <  2:538697/150‑1 (MQ=255)
                                               cttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                                      <  1:528206/150‑1 (MQ=255)
                                                   ctccTGCCAGCCGAAACGGTGAGAACGGGAGGCGGGGTTATCAAAGGTTCCTTCCATCGTCATTATCTGCATCATTTTCTCGGCAA‑CACGCGTGTCTACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCtata                                                                                                  <  2:1074085/148‑1 (MQ=255)
                                                     ccTGCCAGCCGAAACGTTGAGAACGGGAGGCCGGGTTACAAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTACACCTGCGTGGTTGAGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATaaa                                                                                                <  1:372983/150‑1 (MQ=255)
                                                       tGCCAGCCGAAACTGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                              <  2:1633967/150‑1 (MQ=255)
                                                          cAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAACCTCCg                                                                                           >  2:1987342/1‑150 (MQ=255)
                                                            gCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTa                                                                                         <  2:1072696/150‑1 (MQ=255)
                                                            gCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTa                                                                                         <  1:1675608/150‑1 (MQ=255)
                                                            gCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTATCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCTAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTa                                                                                         >  1:822745/1‑150 (MQ=255)
                                                             ccGAAACGGTGAGAACGGGTGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTATCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTAc                                                                                        <  2:304707/150‑1 (MQ=255)
                                                             ccGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTAc                                                                                        <  1:759532/150‑1 (MQ=255)
                                                             ccGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTAc                                                                                        <  1:663683/150‑1 (MQ=255)
                                                               gAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCCCCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACAt                                                                                      >  2:1022633/1‑150 (MQ=255)
                                                                aaccGGTGAGAACGGCAGGCCGGGTTACCAAAGGTTCCGTCCATCGCCATAAACTGCATCATTTTCTCTTCAA‑CACGCTTGTACACCTGCGTGGTTGCGGAGTAGTCGAGATAAATCGTTAATTTCATTGCTCTATAAACTCCGTACATc                                                                                     <  1:81982/147‑1 (MQ=255)
                                                                    ggTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCCCCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTc                                                                                 >  1:877485/1‑150 (MQ=255)
                                                                       gagaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGTGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAAt                                                                              <  1:1340265/150‑1 (MQ=255)
                                                                       gagaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAAt                                                                              <  2:979677/150‑1 (MQ=255)
                                                                        agaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATg                                                                             <  2:896345/150‑1 (MQ=255)
                                                                         gaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTC‑TAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCa                                                                           >  2:680466/1‑150 (MQ=255)
                                                                            cGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAg                                                                         <  2:423730/150‑1 (MQ=255)
                                                                            cGGGAGGACGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGAGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAg                                                                         >  1:78399/1‑150 (MQ=255)
                                                                             gggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAgg                                                                        <  2:1299380/150‑1 (MQ=255)
                                                                              ggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGa                                                                       <  1:788149/150‑1 (MQ=255)
                                                                              ggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGa                                                                       >  1:539667/1‑150 (MQ=255)
                                                                                 ggCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATc                                                                    >  1:1511343/1‑150 (MQ=255)
                                                                                   ccGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAg                                                                  >  2:1869976/1‑150 (MQ=255)
                                                                                    cGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTAGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAgg                                                                 <  1:1474944/150‑1 (MQ=255)
                                                                                          aCCAAAGGGTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCg                                                           >  1:1421147/1‑150 (MQ=255)
                                                                                             aaaGGTTCCTTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCt                                                        >  2:259333/1‑150 (MQ=255)
                                                                                              aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTg                                                       <  2:888762/150‑1 (MQ=255)
                                                                                                    ccGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTCTAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTa                                                 <  2:1317679/150‑1 (MQ=255)
                                                                                                      gTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACg                                               <  2:1124682/150‑1 (MQ=255)
                                                                                                        ccATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGAc                                             >  2:684370/1‑150 (MQ=255)
                                                                                                         cATCGTCATAAACTGCATCATTTTCTCGGCAA‑CCCGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGGACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGAcc                                            >  1:970418/1‑150 (MQ=255)
                                                                                                         cATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCTCTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGAcc                                            >  1:121230/1‑150 (MQ=255)
                                                                                                          aTCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACTGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCg                                           <  1:1822630/150‑1 (MQ=255)
                                                                                                          aTCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCGAGGAATCAGGCTACCGGCTGGATGTACGACCg                                           <  1:467643/150‑1 (MQ=255)
                                                                                                          aTCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTTCCGGCTGGATGTACGACCg                                           >  2:1534158/1‑150 (MQ=255)
                                                                                                          aTCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCg                                           <  1:245455/150‑1 (MQ=255)
                                                                                                           tcatCATAAACTGCATCATTATCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGACGTCATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgt                                          <  2:185607/147‑1 (MQ=255)
                                                                                                            cGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgtg                                         >  1:1738452/1‑150 (MQ=255)
                                                                                                             gTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgtgt                                        >  1:1243888/1‑150 (MQ=255)
                                                                                                              tCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGtt                                       <  1:429251/150‑1 (MQ=255)
                                                                                                               cATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAACCTCCGTACATCACTCAATGCAAGGAATCAGGATACCGGCTGGATGTACGACCGTGttt                                      >  2:147148/1‑150 (MQ=255)
                                                                                                                aTACACTGCAGCATAATCTCGGCAA‑CACGCGGGGCCACCGGCGTGGTGGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTa                                     <  2:368278/150‑1 (MQ=255)
                                                                                                                aTAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTTCTCTATAAACTCCGTACATCACTAAATGCAAGGAATCAGGCTACCGGCTGGATGTAAGACCGTGtata                                     >  1:732475/1‑150 (MQ=255)
                                                                                                                 tAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTAc                                    >  2:856111/1‑150 (MQ=255)
                                                                                                                   aaCTTCATCATTTTCTCGGAAA‑CACTCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACgg                                  <  1:1351019/150‑1 (MQ=255)
                                                                                                                    aCTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGa                                 <  1:835452/150‑1 (MQ=255)
                                                                                                                     cTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                                <  2:968521/150‑1 (MQ=255)
                                                                                                                     atGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                                <  2:1117068/149‑1 (MQ=255)
                                                                                                                        catcatTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAt                             <  2:154866/150‑1 (MQ=255)
                                                                                                                        catcatTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAt                             <  1:31374/150‑1 (MQ=255)
                                                                                                                          tcatTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTTGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAttt                           <  2:568529/150‑1 (MQ=255)
                                                                                                                           catTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTa                          >  2:641580/1‑150 (MQ=255)
                                                                                                                              tttCTCGGCAA‑CACGCGTGTCCACCGGCGTGGTTGCGGAGTTGTCGAGATAAATCTGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGCCCGTGTTTACGGAGTAATTAGCa                       <  2:119316/150‑1 (MQ=255)
                                                                                                                               ttCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGTTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCAc                      >  2:145255/1‑150 (MQ=255)
                                                                                                                                 ctcGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTc                    >  1:1638541/1‑150 (MQ=255)
                                                                                                                                  tcGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCCAGGAATCAGGCTCCCGGCTGGATGTACGCCCGTGTTTACGGAGGATTTAGCACTcc                   >  1:985023/1‑150 (MQ=255)
                                                                                                                                  tcGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTCCGACCGTGTTTACGGAGTATTTAGCACTcc                   >  1:1396073/1‑150 (MQ=255)
                                                                                                                                    ggCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCgg                 <  2:1162725/150‑1 (MQ=255)
                                                                                                                                    cgcaaccacGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCg                  >  2:759047/7‑150 (MQ=255)
                                                                                                                                     gCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATAAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGc                >  2:1729965/1‑150 (MQ=255)
                                                                                                                                      cAA‑CACGCGGGTCCACCGGCGTGGTTGAGGAGTAGTCGAGATAAATCTGTAATTTCATTGCTCTATAAACTCCGTTCATCAGTCAATGCAAGGAATCATGCTACCGGCTGGATGTACGACCGTGTTTACGTTGTATTTAGCACTCCGGcc               <  2:1246212/150‑1 (MQ=255)
                                                                                                                                       aaccacGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCAATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTTGCACTCCGGcc               >  1:1266850/4‑150 (MQ=255)
                                                                                                                                              cgGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCTATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCt        >  1:379428/1‑150 (MQ=255)
                                                                                                                                                ggTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGa      <  2:640285/150‑1 (MQ=255)
                                                                                                                                                  tCCACCGGCGTGGTTGCGGAGTAGGCGAGATAAATCGGTTATTTCATTGCTCTTTAAACTCCGTACATCACTCATTGCAAGGAATCTGGCTAGCGTCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAAt    <  1:980653/150‑1 (MQ=255)
                                                                                                                                                    cACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAAttg  >  2:443032/1‑149 (MQ=255)
                                                                                                                                                    cACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTATATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTc  <  2:1228842/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAA‑CACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTC  >  NC_000913/2661340‑2661636

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