Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,948,179 T→A 6.2% intergenic (‑28/+1) yebC ← / ← nudB putative transcriptional regulator YebC/dihydroneopterin triphosphate diphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,948,1790TA 6.2% 498.2 / 11.2 155intergenic (‑28/+1)yebC/nudBputative transcriptional regulator YebC/dihydroneopterin triphosphate diphosphatase
Reads supporting (aligned to +/- strand):  ref base T (73/72);  new base A (7/3);  total (80/75)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

CCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACC  >  NC_000913/1948033‑1948328
                                                                                                                                                  |                                                                                                                                                     
 cGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAg                                                                                                                                                   >  1:1282/1‑150 (MQ=255)
  gCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAgg                                                                                                                                                  <  1:169882/150‑1 (MQ=255)
   cccAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGc                                                                                                                                                 >  1:2107636/1‑150 (MQ=255)
    ccAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCa                                                                                                                                                <  2:1851475/150‑1 (MQ=255)
      aGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGc                                                                                                                                              >  1:697636/1‑150 (MQ=255)
       gCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCg                                                                                                                                             <  2:228873/150‑1 (MQ=255)
        cTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGt                                                                                                                                            >  2:708926/1‑150 (MQ=255)
         ttAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGtt                                                                                                                                           <  1:1996692/150‑1 (MQ=255)
          tAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTa                                                                                                                                          <  1:510492/150‑1 (MQ=255)
           aGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTaa                                                                                                                                         >  2:570589/1‑150 (MQ=255)
             ccGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAtt                                                                                                                                       <  2:1621695/150‑1 (MQ=255)
              cgcgGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTa                                                                                                                                      >  1:1034080/1‑150 (MQ=255)
                   ttACCAGCTCACGAAAGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGGCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAAc                                                                                                                                 <  2:1026552/150‑1 (MQ=255)
                       cAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGtt                                                                                                                             <  2:2458775/150‑1 (MQ=255)
                       cAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGtt                                                                                                                             <  2:336609/150‑1 (MQ=255)
                        aGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTc                                                                                                                            <  1:241039/150‑1 (MQ=255)
                        aGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTc                                                                                                                            <  1:345111/150‑1 (MQ=255)
                          cTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCaa                                                                                                                          <  2:2110787/150‑1 (MQ=255)
                           tCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAAt                                                                                                                         >  2:420310/1‑150 (MQ=255)
                             aCGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCg                                                                                                                       >  1:2401474/1‑150 (MQ=255)
                                   gATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGcc                                                                                                                 >  1:519474/1‑150 (MQ=255)
                                    aTTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCg                                                                                                                <  2:744466/150‑1 (MQ=255)
                                       ttAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGtt                                                                                                             <  2:275913/150‑1 (MQ=255)
                                        tAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTg                                                                                                            >  1:1533278/1‑150 (MQ=255)
                                        tAGTGAAGANTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTg                                                                                                            >  1:1524313/1‑150 (MQ=255)
                                              aGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAg                                                                                                      >  1:1101626/1‑150 (MQ=255)
                                              aGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAg                                                                                                      <  1:1434806/150‑1 (MQ=255)
                                              aGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTACAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAg                                                                                                      >  2:270467/1‑150 (MQ=255)
                                               gATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAgg                                                                                                     >  1:388185/1‑150 (MQ=255)
                                                 ttttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGAc                                                                                                   <  2:1804297/150‑1 (MQ=255)
                                                  tttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACt                                                                                                  >  2:2239596/1‑150 (MQ=255)
                                                   ttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACtt                                                                                                 <  2:764345/150‑1 (MQ=255)
                                                    tACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTa                                                                                                >  2:1391255/1‑150 (MQ=255)
                                                    tACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTa                                                                                                <  1:430789/150‑1 (MQ=255)
                                                    tACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTa                                                                                                <  2:1282/150‑1 (MQ=255)
                                                     aCCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAg                                                                                               >  2:1225396/1‑150 (MQ=255)
                                                      ccGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGt                                                                                              >  2:1150745/1‑150 (MQ=255)
                                                               gCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGt                                                                                              <  2:1041337/141‑1 (MQ=255)
                                                               acATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGcgccgc                                                                                     >  1:683822/2‑150 (MQ=255)
                                                                cATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGcgccgcc                                                                                    <  2:922712/150‑1 (MQ=255)
                                                                cATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGcgccgcc                                                                                    <  2:1647898/150‑1 (MQ=255)
                                                                 aTCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGcgccgccg                                                                                   >  1:1135443/1‑150 (MQ=255)
                                                                  tCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGcgccgccgc                                                                                  >  2:1743381/1‑150 (MQ=255)
                                                                    cTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcag                                                                                <  2:2445729/150‑1 (MQ=255)
                                                                     tGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTTCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcagc                                                                               >  1:1908919/1‑150 (MQ=255)
                                                                     tGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcagc                                                                               <  1:1109960/150‑1 (MQ=255)
                                                                     tGCGCAGCTTTACGATGTCTGGTGGTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGGTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcagc                                                                               <  1:1288014/150‑1 (MQ=255)
                                                                         cAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAgcg                                                                                         >  2:2186743/1‑136 (MQ=255)
                                                                         cAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAgcg                                                                                         <  1:2186743/136‑1 (MQ=255)
                                                                          aGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGgcg                                                                          <  2:1907246/150‑1 (MQ=255)
                                                                          aGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGgcg                                                                          >  2:1581476/1‑150 (MQ=255)
                                                                           gCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGgcgc                                                                         <  1:2325696/150‑1 (MQ=255)
                                                                           gCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGgcgc                                                                         >  1:1867412/1‑150 (MQ=255)
                                                                            cTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCa                                                                        >  1:1038314/1‑150 (MQ=255)
                                                                                  gATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGc                                                                  <  2:119505/150‑1 (MQ=255)
                                                                                  gATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGc                                                                  >  1:2233264/1‑150 (MQ=255)
                                                                                       cTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGTTTGCACCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAACCCACtt                                                             >  2:236600/1‑150 (MQ=255)
                                                                                       cTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACtt                                                             >  2:1661091/1‑150 (MQ=255)
                                                                                                gCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAg                                                    >  1:2244743/1‑150 (MQ=255)
                                                                                                 cccATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                   >  1:2294541/1‑150 (MQ=255)
                                                                                                  ccATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGAt                                                  >  1:976584/1‑150 (MQ=255)
                                                                                                   cATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATg                                                 <  1:969494/150‑1 (MQ=255)
                                                                                                    aTTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGt                                                >  2:1876223/1‑150 (MQ=255)
                                                                                                      ttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTc                                              <  1:541203/150‑1 (MQ=255)
                                                                                                       tACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTGAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCa                                             >  2:533967/1‑150 (MQ=255)
                                                                                                        aCTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTATTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAg                                            <  2:1538039/150‑1 (MQ=255)
                                                                                                        aCTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAg                                            <  2:574929/150‑1 (MQ=255)
                                                                                                         cTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGt                                           >  2:2339059/1‑150 (MQ=255)
                                                                                                           aTGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGa                                         >  2:2068564/1‑150 (MQ=255)
                                                                                                           aTGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGa                                         >  1:284482/1‑150 (MQ=255)
                                                                                                            tGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGaa                                        <  1:708926/150‑1 (MQ=255)
                                                                                                              aCCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGaaaa                                      >  2:1849768/1‑150 (MQ=255)
                                                                                                              aCCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGaaaa                                      >  2:544667/1‑150 (MQ=255)
                                                                                                               ccTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAAc                                     <  2:1524313/150‑1 (MQ=255)
                                                                                                                   ccATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTTCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATc                                 <  2:1908919/150‑1 (MQ=255)
                                                                                                                    cATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCt                                <  2:297814/150‑1 (MQ=255)
                                                                                                                    cATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCt                                >  2:699231/1‑150 (MQ=255)
                                                                                                                    cATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCt                                <  1:28224/150‑1 (MQ=255)
                                                                                                                     aTAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTg                               >  1:1640709/1‑150 (MQ=255)
                                                                                                                      tAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGc                              >  2:858198/1‑150 (MQ=255)
                                                                                                                      tAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGc                              <  2:683822/150‑1 (MQ=255)
                                                                                                                       aaaaaTATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCTAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGcc                             <  1:828913/150‑1 (MQ=255)
                                                                                                                       aaaaaTATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGcc                             >  2:1977844/1‑150 (MQ=255)
                                                                                                                        aaaaTATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCg                            >  1:790739/1‑150 (MQ=255)
                                                                                                                             atCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGt                       >  2:2090968/1‑150 (MQ=255)
                                                                                                                              tctcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTg                      >  2:148261/1‑150 (MQ=255)
                                                                                                                               ctcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGc                     <  1:113243/150‑1 (MQ=255)
                                                                                                                               ctcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGc                     >  1:567507/1‑150 (MQ=255)
                                                                                                                                tcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCg                    >  1:2064373/1‑150 (MQ=255)
                                                                                                                                  cAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGa                  <  2:815003/150‑1 (MQ=255)
                                                                                                                                   aGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGaa                 <  1:848668/150‑1 (MQ=255)
                                                                                                                                    gATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAg                <  2:860195/150‑1 (MQ=255)
                                                                                                                                       aacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGc                                                       >  1:1407110/5‑108 (MQ=255)
                                                                                                                                       aacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGc                                                       <  2:1407110/104‑1 (MQ=255)
                                                                                                                                       aacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAgcgc             >  1:2322268/5‑150 (MQ=255)
                                                                                                                                        acgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCa            >  1:168362/4‑150 (MQ=255)
                                                                                                                                        acgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCa            >  1:118559/4‑150 (MQ=255)
                                                                                                                                          cctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAg          >  1:496685/1‑150 (MQ=255)
                                                                                                                                          cctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAg          <  2:1348261/150‑1 (MQ=255)
                                                                                                                                                tgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCACGCCACTTGTAAGCCAGATGTTCAGTGAAAAAGATCTGCCGCTCGTGCGGAAGCGCAAGACAGaa    >  2:835827/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACC  >  NC_000913/1948033‑1948328

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