Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,534,891 C→A 8.7% intergenic (+22/+74) nhoA → / ← yddE arylamine N‑acetyltransferase/PF02567 family protein YddE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,534,8910CA 8.7% 255.4 / 11.0 80intergenic (+22/+74)nhoA/yddEarylamine N‑acetyltransferase/PF02567 family protein YddE
Reads supporting (aligned to +/- strand):  ref base C (27/45);  new base A (4/3);  total (31/49)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

ATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTA  >  NC_000913/1534746‑1535040
                                                                                                                                                 |                                                                                                                                                     
aTTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATgcggcg                                                                                                                                                   <  2:1007306/150‑1 (MQ=255)
aTTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATgcggcg                                                                                                                                                   >  2:648022/1‑150 (MQ=255)
   atatGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTaa                                                                                                                                                >  1:94032/1‑150 (MQ=255)
      tGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGGTGCCGGATGCGGCGTAAACg                                                                                                                                             <  2:1024020/150‑1 (MQ=255)
        cTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGcc                                                                                                                                           <  2:180714/150‑1 (MQ=255)
         tgtgATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                          <  2:1621801/150‑1 (MQ=255)
          gtgATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGTAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                         >  1:631931/1‑150 (MQ=255)
           tgATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                        >  1:11581/1‑150 (MQ=255)
            gATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                       >  2:1596511/1‑150 (MQ=255)
              tGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATgcggc                                                                                                                                                    <  1:1220737/135‑1 (MQ=255)
              tGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATgcggc                                                                                                                                                    >  2:1220737/1‑135 (MQ=255)
                cAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                   <  2:1189904/150‑1 (MQ=255)
                 aaGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTTGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGTCGGGAAAATAATTTATGTCAGGTTGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                  <  2:1174359/150‑1 (MQ=255)
                 aaGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                  >  2:1600153/1‑150 (MQ=255)
                 aaGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                  <  2:1469289/150‑1 (MQ=255)
                  aGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCa                                                                                                                                 <  1:500662/150‑1 (MQ=255)
                  aGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCa                                                                                                                                 <  2:94032/150‑1 (MQ=255)
                   gAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATAAGGCGTTTACGCCGCATCCGGCaa                                                                                                                                <  1:701180/150‑2 (MQ=255)
                      cAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACa                                                                                                                             <  2:996755/150‑1 (MQ=255)
                       aGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACAt                                                                                                                            >  1:525645/1‑150 (MQ=255)
                           tGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGc                                                                                                                        >  2:1556632/1‑150 (MQ=255)
                            ggTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGcc                                                                                                                       <  1:52150/150‑1 (MQ=255)
                            ggTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGcc                                                                                                                       >  2:52150/1‑150 (MQ=255)
                             gTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGccc                                                                                                                      >  1:984188/1‑150 (MQ=255)
                              tCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCg                                                                                                                     >  2:1625518/1‑150 (MQ=255)
                               cTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCgg                                                                                                                    <  1:1556632/150‑1 (MQ=255)
                               cTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCgg                                                                                                                    <  1:531801/150‑1 (MQ=255)
                                tGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGt                                                                                                                   <  2:1285801/150‑1 (MQ=255)
                                  ggCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTg                                                                                                                 <  1:1480019/150‑1 (MQ=255)
                                     gTGGATGATGCGAAACATAGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAg                                                                                                              <  2:619151/150‑1 (MQ=255)
                                       ggATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGAc                                                                                                            <  1:1380980/150‑1 (MQ=255)
                                         atgatgCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCt                                                                                                          <  1:1372570/150‑1 (MQ=255)
                                         atgatgCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCt                                                                                                          <  1:648661/150‑1 (MQ=255)
                                               cGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAgg                                                                                                    <  1:1275693/150‑1 (MQ=255)
                                                gAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGt                                                                                                                   <  2:1100091/134‑1 (MQ=255)
                                                gAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGt                                                                                                                   >  1:1100091/1‑134 (MQ=255)
                                                 aaaTATGGCTTTACGGTGGATGAGTTAGCGCTGGTGATGTCGGCGTTTGATACGCACCAGGAGGTGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCg                                                                                                  <  2:1091170/150‑1 (MQ=39)
                                                  aaCATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGGTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGa                                                                                                 <  1:1622848/150‑1 (MQ=255)
                                                   aCATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGaa                                                                                                >  2:868243/1‑150 (MQ=255)
                                                          tttACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATc                                                                                         <  1:1348687/150‑1 (MQ=255)
                                                           ttACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCa                                                                                        <  2:1328901/150‑1 (MQ=255)
                                                           ttACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCa                                                                                        <  2:361798/150‑1 (MQ=255)
                                                              ccGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGAc                                                                                     >  1:1026341/1‑150 (MQ=255)
                                                                gTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACaa                                                                                   >  2:53500/1‑150 (MQ=255)
                                                                 tGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAAc                                                                                  <  1:678346/150‑1 (MQ=255)
                                                                          tAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCaga                                                                         >  2:592463/1‑150 (MQ=255)
                                                                           aGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCagag                                                                        >  2:424748/1‑150 (MQ=255)
                                                                              gcTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTc                                                                     <  1:112629/150‑1 (MQ=255)
                                                                                  gTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAAtt                                                                 <  1:736057/150‑1 (MQ=255)
                                                                                     aTGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGcc                                                              <  2:468176/150‑1 (MQ=255)
                                                                                      tGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGccc                                                             <  1:1600153/150‑1 (MQ=255)
                                                                                       ggcggcGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCa                                                            >  1:1453194/1‑150 (MQ=255)
                                                                                        gcggcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCAt                                                           >  2:372668/1‑150 (MQ=255)
                                                                                         cggcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCAtt                                                          <  1:53500/150‑1 (MQ=255)
                                                                                         cggcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCAtt                                                          <  1:13712/150‑1 (MQ=255)
                                                                                             gTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGc                                                      <  1:514984/150‑1 (MQ=255)
                                                                                              tttGATTGGCACCCGGAGGCGGGAAAATAATTTATGGCTGGTTGCCGGATGCGGCGTAAAAGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCGCATTCGGCa                                                     <  1:1609182/150‑1 (MQ=255)
                                                                                              tttGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCa                                                     <  1:1200664/150‑1 (MQ=255)
                                                                                                tGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGTATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATg                                                   <  2:1108521/150‑1 (MQ=255)
                                                                                                tGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATg                                                   >  1:637749/1‑150 (MQ=255)
                                                                                                tGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATg                                                   <  1:424748/150‑1 (MQ=255)
                                                                                                  aTACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGaa                                                 <  1:372668/150‑1 (MQ=255)
                                                                                                      gCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaa                                             >  1:640711/1‑150 (MQ=255)
                                                                                                       cACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACaaa                                            <  2:952339/150‑1 (MQ=255)
                                                                                                         cccGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAAt                                          >  1:1545021/1‑150 (MQ=255)
                                                                                                          ccGGAGGCGGGAAAATAATTTATGTCAAGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATc                                         >  2:57387/1‑150 (MQ=255)
                                                                                                              aGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCg                                     <  1:868243/150‑1 (MQ=255)
                                                                                                                   ggAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCtgtg                                <  1:1596511/150‑1 (MQ=255)
                                                                                                                     taaaTAATTTATTGCATGTTGGCGGATGGGGCGGAAACGCCTTATCCGGCATAAATTAGCCCGGGTGTAGTCCTGACCTGCGAAGCGCATCAGACAACGGTTTGATAGGTCAATTTCCCATTCGGCATGAAAAAAAATCACCGCTGTGcc                              <  2:11581/149‑1 (MQ=255)
                                                                                                                     aaaaTAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAATTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGcc                              <  2:1546011/150‑1 (MQ=255)
                                                                                                                            tttATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATg                       >  2:289176/1‑150 (MQ=255)
                                                                                                                            tttATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCAGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATg                       <  1:1426905/150‑1 (MQ=255)
                                                                                                                               aTGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTg                    <  2:1287950/150‑1 (MQ=255)
                                                                                                                               aTGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCAGACATTAGCCCGGTTGTAGCCCTGACAGGCGCAGCGCATCCGACAACGGGGTCAGAGTTCAATTGCCCAGTCGGCATGATACAGAATCCCCGCTGTGCCTGAAGTgggg                    >  2:470639/1‑148 (MQ=255)
                                                                                                                                tGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGa                   >  1:1465578/1‑150 (MQ=255)
                                                                                                                                 gTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGAc                  >  1:1293819/1‑150 (MQ=255)
                                                                                                                                    aGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACttt               <  2:759606/150‑1 (MQ=255)
                                                                                                                                     ggTTGCCGGATGCGGCGTAAACGCCTTATCCGGCACACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACtttt              >  1:1418933/1‑150 (MQ=255)
                                                                                                                                        tGCCGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTcc           >  1:954800/1‑150 (MQ=255)
                                                                                                                                        tGCCGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTcc           >  1:36464/1‑150 (MQ=255)
                                                                                                                                         gCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCg          <  1:962043/150‑1 (MQ=255)
                                                                                                                                         gCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCg          <  1:1187326/150‑1 (MQ=255)
                                                                                                                                         gCCGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACTGCTGTGCCTGAAATGGTGACTTTTTCCg          >  1:683563/1‑150 (MQ=255)
                                                                                                                                           cGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCTCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCggt        >  2:1472546/1‑150 (MQ=255)
                                                                                                                                                 cggcgTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGGTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTa  >  1:496533/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
ATTATATGCTGTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTA  >  NC_000913/1534746‑1535040

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