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→GTCiscScysteine 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.

AGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTC  >  NC_000913/2661340‑2661636
                                                                                                                                                    |                                                                                                                                                     
aGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCAc                                                                                                                                                     >  1:1887840/1‑150 (MQ=255)
 gcgcCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCAcc                                                                                                                                                    <  1:620102/150‑1 (MQ=255)
       aCCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTg                                                                                                                                              <  1:1526695/150‑1 (MQ=255)
        ccAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTgg                                                                                                                                             >  1:119316/1‑150 (MQ=255)
        ccAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTgg                                                                                                                                             <  1:1516949/150‑1 (MQ=255)
           taTCCGCAATCTGATTACGGTAGATATCTCCCGCTTCTGCAGCCTGCCAGCCGAAACGGTGAGATCGGGAGGCCGGGTTACCACAGGTTCAGGCCATCGTCATAATCTGCATCATTTTCTCGGCAACACGCGGGGCCCCCGTCGTGGTTg                                                                                                                                          <  1:903744/149‑1 (MQ=255)
           gATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTg                                                                                                                                          >  1:1394919/1‑150 (MQ=255)
           gATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCAACCGGCGGGGTTg                                                                                                                                          >  2:1234556/1‑150 (MQ=255)
             tCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACCACGCCGGTGGACCCGCGTgttgcc                                                                                                                                         <  1:1743277/150‑6 (MQ=255)
              ccgCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCATGGCAGCCGAAACGGTGAGAACCGGAGGCCGGGTTACCAAACGTTCCGTTCATCGTCATCAACTGCATTATTTTCTCGGCAACACGGGGGTCCACCGGCGTGGTTGCgg                                                                                                                                       <  2:250580/148‑1 (MQ=255)
              cGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                       <  2:722770/150‑1 (MQ=255)
              cGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                       <  2:1843839/150‑1 (MQ=255)
              cGGCAATCTGATCACGGGCGATAGCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGGGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCgg                                                                                                                                       <  1:1203597/150‑1 (MQ=255)
                gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                                     <  1:778722/150‑1 (MQ=255)
                gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGag                                                                                                                                     <  1:1118945/150‑1 (MQ=255)
                gCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGTGTTGCCGag                                                                                                                                     <  2:446213/150‑1 (MQ=38)
                 cAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGagt                                                                                                                                    >  2:645162/1‑150 (MQ=255)
                 cAATCTGATTACGGGCGATATCTACCGCTTCTTAAGCCTGCCAGCCGAAACTGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGagt                                                                                                                                    <  2:1923958/150‑1 (MQ=255)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGTAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCCGCAACCACGCCGGTGGACCCGCGTGTTGCCGagaaaatg                                                                                                                               <  1:1386673/150‑7 (MQ=37)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCGACCGGCGGGGTTGCGGAGGAGTCg                                                                                                                               >  2:1981826/1‑150 (MQ=255)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                               <  2:1689455/150‑1 (MQ=255)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                               >  2:759532/1‑150 (MQ=255)
                      tGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCg                                                                                                                               >  2:1610264/1‑150 (MQ=255)
                       gATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCCGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCga                                                                                                                              <  1:1893438/150‑1 (MQ=255)
                            cGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGCAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATaa                                                                                                                         >  2:1454340/1‑150 (MQ=255)
                              ggCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAAt                                                                                                                       >  1:456428/1‑150 (MQ=255)
                               gCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATc                                                                                                                      >  2:660753/1‑150 (MQ=255)
                                      cTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGGGGTTGCGGAGTAGTCGAGATAACTCGGTAAtt                                                                                                               >  1:520981/1‑150 (MQ=255)
                                       tACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATTAATCGGTAAttt                                                                                                              >  1:1759665/1‑150 (MQ=255)
                                         ccGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCa                                                                                                            <  2:116924/150‑1 (MQ=255)
                                          cGCTGCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCAt                                                                                                           <  2:538697/150‑1 (MQ=255)
                                               cttcAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGctc                                                                                                      <  1:528206/150‑1 (MQ=255)
                                                   ctccTGCCAGCCGAAACGGTGAGAACGGGAGGCGGGGTTATCAAAGGTTCCTTCCATCGTCATTATCTGCATCATTTTCTCGGCAACACGCGTGTCTACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCtata                                                                                                  <  2:1074085/148‑1 (MQ=255)
                                                     ccTGCCAGCCGAAACGTTGAGAACGGGAGGCCGGGTTACAAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTACACCTGCGTGGTTGAGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATaaa                                                                                                <  1:372983/150‑1 (MQ=255)
                                                       tGCCAGCCGAAACTGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACt                                                                                              <  2:1633967/150‑1 (MQ=255)
                                                          cAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAACCTCCg                                                                                           >  2:1987342/1‑150 (MQ=255)
                                                            gCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTa                                                                                         <  2:1072696/150‑1 (MQ=255)
                                                            gCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTa                                                                                         <  1:1675608/150‑1 (MQ=255)
                                                            gCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTATCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCTAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTa                                                                                         >  1:822745/1‑150 (MQ=255)
                                                             ccGAAACGGTGAGAACGGGTGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTATCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTAc                                                                                        <  2:304707/150‑1 (MQ=255)
                                                             ccGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTAc                                                                                        <  1:759532/150‑1 (MQ=255)
                                                             ccGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTAc                                                                                        <  1:663683/150‑1 (MQ=255)
                                                               gAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCCCCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACAt                                                                                      >  2:1022633/1‑150 (MQ=255)
                                                                aaccGGTGAGAACGGCAGGCCGGGTTACCAAAGGTTCCGTCCATCGCCATAAACTGCATCATTTTCTCTTCAACACGCTTGTACACCTGCGTGGTTGCGGAGTAGTCGAGATAAATCGTTAATTTCATTGCTCTATAAACTCCGTACATc                                                                                     <  1:81982/147‑1 (MQ=255)
                                                                    ggTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCCCCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTc                                                                                 >  1:877485/1‑150 (MQ=255)
                                                                       gagaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGTGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAAt                                                                              <  1:1340265/150‑1 (MQ=255)
                                                                       gagaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAAt                                                                              <  2:979677/150‑1 (MQ=255)
                                                                        agaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATg                                                                             <  2:896345/150‑1 (MQ=255)
                                                                         gaACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCTAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCa                                                                           >  2:680466/1‑150 (MQ=255)
                                                                            cGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAg                                                                         <  2:423730/150‑1 (MQ=255)
                                                                            cGGGAGGACGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGAGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAg                                                                         >  1:78399/1‑150 (MQ=255)
                                                                             gggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAgg                                                                        <  2:1299380/150‑1 (MQ=255)
                                                                              ggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGa                                                                       <  1:788149/150‑1 (MQ=255)
                                                                              ggAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGa                                                                       >  1:539667/1‑150 (MQ=255)
                                                                                 ggCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATc                                                                    >  1:1511343/1‑150 (MQ=255)
                                                                                   ccGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAg                                                                  >  2:1869976/1‑150 (MQ=255)
                                                                                    cGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTAGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAgg                                                                 <  1:1474944/150‑1 (MQ=255)
                                                                                          aCCAAAGGGTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCg                                                           >  1:1421147/1‑150 (MQ=255)
                                                                                             aaaGGTTCCTTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCt                                                        >  2:259333/1‑150 (MQ=255)
                                                                                              aaGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTg                                                       <  2:888762/150‑1 (MQ=255)
                                                                                                    ccGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTCTAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTa                                                 <  2:1317679/150‑1 (MQ=255)
                                                                                                      gTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACg                                               <  2:1124682/150‑1 (MQ=255)
                                                                                                        ccATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGAc                                             >  2:684370/1‑150 (MQ=255)
                                                                                                         cATCGTCATAAACTGCATCATTTTCTCGGCAACCCGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGGACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGAcc                                            >  1:970418/1‑150 (MQ=255)
                                                                                                         cATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCTCTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGAcc                                            >  1:121230/1‑150 (MQ=255)
                                                                                                          aTCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACTGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCg                                           <  1:1822630/150‑1 (MQ=255)
                                                                                                          aTCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCGAGGAATCAGGCTACCGGCTGGATGTACGACCg                                           <  1:467643/150‑1 (MQ=255)
                                                                                                          aTCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTTCCGGCTGGATGTACGACCg                                           >  2:1534158/1‑150 (MQ=255)
                                                                                                          aTCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCg                                           <  1:245455/150‑1 (MQ=255)
                                                                                                           tcatCATAAACTGCATCATTATCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGACGTCATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgt                                          <  2:185607/147‑1 (MQ=255)
                                                                                                            cGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgtg                                         >  1:1738452/1‑150 (MQ=255)
                                                                                                             gTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCgtgt                                        >  1:1243888/1‑150 (MQ=255)
                                                                                                              tCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGtt                                       <  1:429251/150‑1 (MQ=255)
                                                                                                               cATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAACCTCCGTACATCACTCAATGCAAGGAATCAGGATACCGGCTGGATGTACGACCGTGttt                                      >  2:147148/1‑150 (MQ=255)
                                                                                                                aTACACTGCAGCATAATCTCGGCAACACGCGGGGCCACCGGCGTGGTGGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTa                                     <  2:368278/150‑1 (MQ=255)
                                                                                                                aTAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTTCTCTATAAACTCCGTACATCACTAAATGCAAGGAATCAGGCTACCGGCTGGATGTAAGACCGTGtata                                     >  1:732475/1‑150 (MQ=255)
                                                                                                                 tAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTAc                                    >  2:856111/1‑150 (MQ=255)
                                                                                                                   aaCTTCATCATTTTCTCGGAAACACTCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACgg                                  <  1:1351019/150‑1 (MQ=255)
                                                                                                                    aCTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGa                                 <  1:835452/150‑1 (MQ=255)
                                                                                                                     cTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                                <  2:968521/150‑1 (MQ=255)
                                                                                                                     atGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAg                                <  2:1117068/149‑1 (MQ=255)
                                                                                                                        catcatTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAt                             <  2:154866/150‑1 (MQ=255)
                                                                                                                        catcatTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAt                             <  1:31374/150‑1 (MQ=255)
                                                                                                                          tcatTTTCTCGGCAACACGCGGGTCCACCGGCGTTGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTAttt                           <  2:568529/150‑1 (MQ=255)
                                                                                                                           catTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTa                          >  2:641580/1‑150 (MQ=255)
                                                                                                                              tttCTCGGCAACACGCGTGTCCACCGGCGTGGTTGCGGAGTTGTCGAGATAAATCTGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGCCCGTGTTTACGGAGTAATTAGCa                       <  2:119316/150‑1 (MQ=255)
                                                                                                                               ttCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGTTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCAc                      >  2:145255/1‑150 (MQ=255)
                                                                                                                                 ctcGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTc                    >  1:1638541/1‑150 (MQ=255)
                                                                                                                                  tcGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCCAGGAATCAGGCTCCCGGCTGGATGTACGCCCGTGTTTACGGAGGATTTAGCACTcc                   >  1:985023/1‑150 (MQ=255)
                                                                                                                                  tcGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTCCGACCGTGTTTACGGAGTATTTAGCACTcc                   >  1:1396073/1‑150 (MQ=255)
                                                                                                                                    ggCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCgg                 <  2:1162725/150‑1 (MQ=255)
                                                                                                                                    cgcaaccacGCCGGTGGACCCGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCg                  >  2:759047/7‑150 (MQ=255)
                                                                                                                                     gCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATAAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGc                >  2:1729965/1‑150 (MQ=255)
                                                                                                                                      cAACACGCGGGTCCACCGGCGTGGTTGAGGAGTAGTCGAGATAAATCTGTAATTTCATTGCTCTATAAACTCCGTTCATCAGTCAATGCAAGGAATCATGCTACCGGCTGGATGTACGACCGTGTTTACGTTGTATTTAGCACTCCGGcc               <  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)
                                                                                                                                                    |                                                                                                                                                     
AGCGCCGACCAGATCGGCAATCTGATTACGGGCGATATCTACCGCTTCTTCAGCCTGCCAGCCGAAACGGTGAGAACGGGAGGCCGGGTTACCAAAGGTTCCGTCCATCGTCATAAACTGCATCATTTTCTCGGCAACACGCGGGTCCACCGGCGTGGTTGCGGAGTAGTCGAGATAAATCGGTAATTTCATTGCTCTATAAACTCCGTACATCACTCAATGCAAGGAATCAGGCTACCGGCTGGATGTACGACCGTGTTTACGGAGTATTTAGCACTCCGGCCTGATTCTGAATTC  >  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: