Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,750,086 A→C 5.3% intergenic (+26/+29) yfjD → / ← grpE putative inner membrane protein/nucleotide exchange factor GrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,750,0860AC 5.3% 242.3 / 5.0 75intergenic (+26/+29)yfjD/grpEputative inner membrane protein/nucleotide exchange factor GrpE
Reads supporting (aligned to +/- strand):  ref base A (39/32);  new base C (2/2);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

AGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCAT  >  NC_000913/2749943‑2750234
                                                                                                                                               |                                                                                                                                                    
aGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGTTAAAAGTTTTTGCTGTGAAAGCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGTCATTCGCATGCAGGGCCg                                                                                                                                                <  2:626380/150‑1 (MQ=255)
aGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCGGGCCg                                                                                                                                                <  2:472922/150‑1 (MQ=255)
aGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCGGGCCg                                                                                                                                                <  2:465774/150‑1 (MQ=255)
   cACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGa                                                                                                                                             >  1:1584884/1‑150 (MQ=255)
   cACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGa                                                                                                                                             >  2:1620696/1‑150 (MQ=255)
      ccGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAatt                                                                                                                                          <  1:1536370/150‑1 (MQ=255)
       cgcgTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAatta                                                                                                                                         <  1:1620696/150‑1 (MQ=255)
       cgcgTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAatta                                                                                                                                         <  2:181474/150‑1 (MQ=255)
         cgtgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAattatt                                                                                                                                       <  2:1537530/150‑1 (MQ=255)
         cgtgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAattatt                                                                                                                                       <  2:537653/150‑1 (MQ=255)
          gtgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAAttatta                                                                                                                                      <  1:1613871/150‑1 (MQ=255)
            gcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACg                                                                                                                                    <  2:840317/150‑1 (MQ=255)
              gtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGaa                                                                                                                                  >  2:653164/1‑150 (MQ=255)
               tATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGaaa                                                                                                                                 >  2:221122/1‑150 (MQ=255)
                aTTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAg                                                                                                                                >  1:1277272/1‑150 (MQ=255)
                aTTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAg                                                                                                                                >  2:112242/1‑150 (MQ=255)
                      gAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGaaa                                                                                                                          <  1:1713541/150‑1 (MQ=255)
                          aCGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTaa                                                                                                                      >  1:1664113/1‑150 (MQ=255)
                            gatatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGc                                                                                                                    <  1:554093/150‑1 (MQ=255)
                                 cgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTg                                                                                                               <  1:221122/150‑1 (MQ=255)
                                     atTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCttt                                                                                                           <  2:794525/150‑1 (MQ=255)
                                        ctcGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGTCCGTGAATTATTACGAAAGCAAAAGCTTAATTTCTGCTTTCgt                                                                                                        <  2:308672/150‑2 (MQ=255)
                                           gACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTAc                                                                                                     <  2:1034626/150‑1 (MQ=255)
                                                   ggACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCa                                                                                             >  1:1830862/1‑150 (MQ=255)
                                                   ggACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCa                                                                                             >  1:69935/1‑150 (MQ=255)
                                                   cgACAATAAGATTAAGCAGATAAAAGTTATACCTGTGAAACCGCTGCGCCAGAGTGTCGCTGAGTAAACAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCa                                                                                             >  2:1363215/2‑150 (MQ=255)
                                                      cAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATcg                                                                                          >  2:1905483/1‑150 (MQ=255)
                                                         tatGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATcgccg                                                                                       >  1:1955253/1‑150 (MQ=255)
                                                           tGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCa                                                                                     >  2:27282/1‑150 (MQ=255)
                                                           tGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCa                                                                                     >  1:202294/1‑150 (MQ=255)
                                                           tGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCa                                                                                     <  2:838378/150‑1 (MQ=255)
                                                           tGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCa                                                                                     >  2:1941725/1‑150 (MQ=255)
                                                             aTTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGGGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCTGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACg                                                                                   <  2:1335934/150‑1 (MQ=255)
                                                             aTTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACg                                                                                   <  1:222596/150‑1 (MQ=255)
                                                               tAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGaa                                                                                 >  1:590693/1‑150 (MQ=255)
                                                                aaGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAAt                                                                                <  1:469807/150‑1 (MQ=255)
                                                                 aGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATc                                                                               >  2:1292289/1‑150 (MQ=255)
                                                                  gCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCg                                                                              >  1:1084238/1‑150 (MQ=255)
                                                                    aGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTa                                                                            >  1:255911/1‑150 (MQ=255)
                                                                      gTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACg                                                                          <  2:245715/150‑1 (MQ=255)
                                                                        aaaaGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGAc                                                                        <  2:1824954/150‑1 (MQ=255)
                                                                                      tgAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATaa                                                          >  2:283054/1‑150 (MQ=255)
                                                                                       gAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAAc                                                         <  1:283054/150‑1 (MQ=255)
                                                                                       gAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAAc                                                         >  2:359793/1‑150 (MQ=255)
                                                                                       gAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAAc                                                         >  2:366845/1‑150 (MQ=255)
                                                                                        aaaCCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAAcc                                                        >  1:1024431/1‑150 (MQ=255)
                                                                                              cTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTg                                                  >  2:1150653/1‑150 (MQ=255)
                                                                                               tGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGc                                                 <  2:1955253/150‑1 (MQ=255)
                                                                                               tGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGc                                                 <  2:1584884/150‑1 (MQ=255)
                                                                                               ggcgcGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGc                                                 <  2:912631/149‑1 (MQ=255)
                                                                                                gcgcGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCa                                                <  1:610133/150‑1 (MQ=255)
                                                                                                   cgAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATaa                                             <  1:653164/150‑1 (MQ=255)
                                                                                                      gagTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGc                                          <  1:1905483/150‑1 (MQ=255)
                                                                                                          gtgtCGGAGTAACGAAAACTGCACGGCATTATCATGAAGGGACGTGAATTAATACCAAAGCAAAAATTAAAAATTATCTTTCGCTAAAGTAACCATCGCCGCACGAATCGTACGACCAATCAGCGAATAACCATTCTCCATAAAGCACAg                                      >  2:1769458/1‑150 (MQ=255)
                                                                                                           tgGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGt                                     >  1:386188/1‑150 (MQ=255)
                                                                                                            ggCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAATGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTa                                    >  1:1652298/1‑150 (MQ=255)
                                                                                                            ggCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTa                                    >  2:1196771/1‑150 (MQ=255)
                                                                                                                  gTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGccc                                        <  2:1788006/140‑1 (MQ=255)
                                                                                                                  gTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGccc                                        >  1:1788006/1‑140 (MQ=255)
                                                                                                                    aaCGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCt                            >  1:864635/1‑150 (MQ=255)
                                                                                                                     aCGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTg                           >  1:1202832/1‑150 (MQ=255)
                                                                                                                          caCGGCCCTGCATGCGAATGCCGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCa                      >  1:1030852/2‑150 (MQ=255)
                                                                                                                            cGGCCCTGCATGCGAATGCCGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCCCACGAATCGTACGACCATTCAGCGTATAACCCATCTACATAATGCCCAGTACCTTACCTGCCGCAAc                    >  2:878240/1‑150 (MQ=255)
                                                                                                                            cGGCCCGGCATTCGCATGCAGGGCCGTGAATTGTTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAAc                    >  1:1255538/1‑150 (MQ=255)
                                                                                                                            cGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAAc                    >  1:243052/1‑150 (MQ=255)
                                                                                                                            cGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAAc                    >  1:283088/1‑150 (MQ=255)
                                                                                                                             ggCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACg                   <  1:1819849/150‑1 (MQ=255)
                                                                                                                             ggCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACg                   <  2:768836/150‑1 (MQ=255)
                                                                                                                                ccGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGtca                >  1:494056/1‑150 (MQ=255)
                                                                                                                                   gCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGtcatca             <  1:112242/150‑1 (MQ=255)
                                                                                                                                       tCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGAtt         <  2:1084238/150‑1 (MQ=255)
                                                                                                                                        cGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCACACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTc        >  1:161558/1‑150 (MQ=255)
                                                                                                                                         gCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTcc       >  2:341096/1‑150 (MQ=255)
                                                                                                                                             gCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTccacca   >  1:1963271/1‑150 (MQ=255)
                                                                                                                                              cAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCAt  >  2:1595131/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCAT  >  NC_000913/2749943‑2750234

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: