Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,899 A→C 8.7% intergenic (+98/+41) rplY → / ← yejK 50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,282,8990AC 8.7% 235.9 / 9.9 87intergenic (+98/+41)rplY/yejK50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base A (43/35);  new base C (3/5);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAAC  >  NC_000913/2282752‑2283048
                                                                                                                                                   |                                                                                                                                                     
ggcACGCGTCCAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCGGGCGTGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAgg                                                                                                                                                     <  1:1039671/148‑1 (MQ=255)
gTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAgg                                                                                                                                                     <  1:352665/150‑1 (MQ=255)
 tCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGc                                                                                                                                                    <  1:1593637/150‑1 (MQ=255)
  cACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCg                                                                                                                                                   <  1:1470755/150‑1 (MQ=255)
    cccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGtt                                                                                                                                                 >  2:21272/1‑150 (MQ=255)
     ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTc                                                                                                                                                <  2:307140/150‑1 (MQ=255)
      cGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGtgaa                                                                                                                                               >  2:810642/1‑147 (MQ=255)
      cGTACAAACCGAAGCTGCAGCACATAGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTACGGCGGGGTTTTTTGTATCTGCAGAATATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCa                                                                                                                                               >  2:332372/1‑150 (MQ=255)
          cAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgcc                                                                                                                                           <  2:1311041/150‑1 (MQ=255)
          cAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCATGATTACGGTATTGCTATTTTTTGCAGGCAAGATAAGGCGTTCAcgcc                                                                                                                                           >  1:211428/1‑150 (MQ=255)
               cGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATc                                                                                                                                      >  2:1633115/1‑150 (MQ=255)
                gAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCt                                                                                                                                     <  2:6188/150‑1 (MQ=255)
                  aGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTgg                                                                                                                                   <  1:1492582/150‑1 (MQ=255)
                  aGCTGCACCACTTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCCCTCCCGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTgg                                                                                                                                   <  1:498591/150‑1 (MQ=255)
                   gCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGc                                                                                                                                  <  1:1605419/150‑1 (MQ=255)
                      gcagcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTTTTGCTATTTTTTGCACGCCAGATGCGGCGTGAACGCCTTATCTGGcctg                                                                                                                               <  1:1131571/150‑4 (MQ=255)
                      gcagcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcctg                                                                                                                               <  1:1326959/150‑4 (MQ=255)
                      gcagcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcctg                                                                                                                               <  1:810642/150‑4 (MQ=255)
                      gcagcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAGGCGTTCACGCCGCATCTGGCATTc                                                                                                                              <  1:758902/150‑1 (MQ=255)
                           acaTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                                                                                          >  2:876773/1‑150 (MQ=255)
                            caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                                                                                                         <  1:739205/150‑1 (MQ=255)
                            caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                                                                                                         <  2:428348/150‑1 (MQ=255)
                            caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAac                                                                                                                         <  2:613667/150‑1 (MQ=255)
                             aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGCTTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAaca                                                                                                                        >  2:871927/1‑150 (MQ=255)
                               cGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAAcaca                                                                                                                      <  1:889285/150‑1 (MQ=255)
                                gACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACaa                                                                                                                     >  1:954488/1‑150 (MQ=255)
                                     cgttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAATTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCg                                                                                                                <  1:252217/150‑1 (MQ=255)
                                       ttcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGAt                                                                                                              <  2:276013/150‑1 (MQ=255)
                                                    gCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAg                                                                                                 >  2:802979/1‑150 (MQ=255)
                                                      tGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTTTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGAc                                                                                               >  2:1662516/1‑150 (MQ=255)
                                                      tGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGAc                                                                                               >  1:394762/1‑150 (MQ=255)
                                                       gAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                                                                                              >  2:719715/1‑150 (MQ=255)
                                                       gAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                                                                                              >  2:497291/1‑150 (MQ=255)
                                                         aTAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                                                                                            >  1:112236/1‑150 (MQ=255)
                                                           cagatacaaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCACCAGACGTCCg                                                                                          >  1:63189/8‑150 (MQ=255)
                                                            aGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCgg                                                                                         >  2:1669966/1‑150 (MQ=255)
                                                            aGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCgg                                                                                         >  2:796792/1‑150 (MQ=255)
                                                            aGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCgg                                                                                         <  1:1382879/150‑1 (MQ=255)
                                                                gTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCt                                                                                     >  1:527790/1‑150 (MQ=255)
                                                                gTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCt                                                                                     >  2:830734/1‑150 (MQ=255)
                                                                 caaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg                                                                                    >  2:245115/2‑150 (MQ=255)
                                                                 caaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGACCGGCGCTg                                                                                    >  2:870548/2‑150 (MQ=255)
                                                                  aaaaaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc                                                                                   >  1:1418486/1‑150 (MQ=255)
                                                                  aaaaaaCCCCGCTCCGGCGGGGTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGTGTTCACGCCGCATCTGGCAGTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCCGCa                                                                                  >  1:1616931/1‑150 (MQ=255)
                                                                  aaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc                                                                                   >  2:1718957/1‑150 (MQ=255)
                                                                  aaaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATAACGGTATAGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGc                                                                                   >  2:1278350/1‑150 (MQ=255)
                                                                   aaaaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCa                                                                                  >  2:1671321/1‑150 (MQ=255)
                                                                   aaaaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCa                                                                                  >  1:1674897/1‑150 (MQ=255)
                                                                   aaaaaCCCCGCCGGAGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCCCGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCa                                                                                  >  1:1566209/1‑150 (MQ=255)
                                                                       aCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTg                                                                              >  1:1161456/1‑150 (MQ=255)
                                                                         cccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGt                                                                            >  2:1134958/1‑150 (MQ=255)
                                                                         cccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTCTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGt                                                                            <  1:724133/150‑1 (MQ=255)
                                                                             cTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTcgcg                                                                        >  2:1698872/1‑150 (MQ=255)
                                                                                   cggGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGGGCAAAtt                                                                  >  1:997411/1‑150 (MQ=255)
                                                                                       gTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGCTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGGCCGGCGCTGCAATTGGTCGCGCAAATTcggc                                                              <  1:1523280/150‑1 (MQ=255)
                                                                                       gTTTTTTGTATCTGCAGATTATGCCTGAGTACGGTATTGCTATTTTTTGCAGGCCAGATACGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcggc                                                              <  1:476059/150‑1 (MQ=255)
                                                                                          ttttGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCggcggc                                                           <  1:1734112/150‑1 (MQ=255)
                                                                                             tGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTg                                                        >  2:1464246/1‑150 (MQ=255)
                                                                                             tGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTg                                                        <  2:1161456/150‑1 (MQ=255)
                                                                                              gTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGATCGCGCAAATTCGGCGGCGTGc                                                       >  2:504004/1‑150 (MQ=255)
                                                                                                aTCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCt                                                     <  2:1095367/150‑1 (MQ=255)
                                                                                                 tCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCtt                                                    >  1:971709/1‑150 (MQ=255)
                                                                                                   tGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTa                                                  <  2:112236/150‑1 (MQ=255)
                                                                                                   tGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTa                                                  >  2:1291284/1‑150 (MQ=255)
                                                                                                   tGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAAGACGGCGGCGTGccgtga                                                  >  2:305962/1‑146 (MQ=255)
                                                                                                       gATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATgg                                              <  2:1478453/150‑1 (MQ=255)
                                                                                                        aTTATGCCTGATTACTGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGt                                             <  1:1677519/150‑1 (MQ=255)
                                                                                                            tGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAg                                         >  1:1652772/1‑150 (MQ=255)
                                                                                                             gCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAgg                                        <  1:876773/150‑1 (MQ=255)
                                                                                                             gCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAgg                                        <  2:958789/150‑1 (MQ=255)
                                                                                                             gCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAgg                                        <  2:1533259/150‑1 (MQ=255)
                                                                                                              ccTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGt                                       >  2:1394891/1‑150 (MQ=255)
                                                                                                                 gATTACGGTATTGCTATTTTTAGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTAtt                                    <  1:332372/150‑2 (MQ=255)
                                                                                                                         tATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTg                            <  2:944683/150‑1 (MQ=255)
                                                                                                                         tATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTg                            <  2:954488/150‑1 (MQ=255)
                                                                                                                              cTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCAGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATcc                       >  2:1107733/1‑150 (MQ=255)
                                                                                                                               tATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATccc                      >  1:968832/1‑150 (MQ=255)
                                                                                                                                aTTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCa                     >  1:1426589/1‑150 (MQ=255)
                                                                                                                                   ttttGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCagaa                  >  1:349821/1‑150 (MQ=255)
                                                                                                                                     ttGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCagaaga                <  2:1169548/150‑1 (MQ=255)
                                                                                                                                        cAGGCCAGATATGGCGTTCACGCCGCATCTGGCTTTCAGAACACAATACCGTTTAATTGCCGCGAGAGGTCCGGCGGTGCAATTGGTCGCGCATATTCGGCGGCGTGCCTTTTATGGTCAAGGTATGGGTTGCTGGATCCCAGAAGATGc             <  2:431569/150‑1 (MQ=255)
                                                                                                                                        cAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGc             <  2:1187515/150‑1 (MQ=255)
                                                                                                                                          ggCAAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAACCCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTTGTCGCGCTAATTCGGAGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGt           <  2:863946/150‑1 (MQ=255)
                                                                                                                                           gCCAGATGCGGCGTGAACGCCTTATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGtt          >  2:1475946/1‑150 (MQ=255)
                                                                                                                                           gCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGtt          <  2:1136561/150‑1 (MQ=255)
                                                                                                                                            ccAGATGCGGCGTGAACGCCTTATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTc         >  1:1018058/1‑150 (MQ=255)
                                                                                                                                             cAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGACGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCa        >  2:1307943/1‑150 (MQ=255)
                                                                                                                                                  aaGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTaa   >  2:357619/1‑150 (MQ=255)
                                                                                                                                                   aGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAAc  >  1:233571/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAAC  >  NC_000913/2282752‑2283048

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