Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,534,905 T→C 10.1% intergenic (+36/+60) 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,9050TC10.1% 213.3 / 13.2 79intergenic (+36/+60)nhoA/yddEarylamine N‑acetyltransferase/PF02567 family protein YddE
Reads supporting (aligned to +/- strand):  ref base T (26/45);  new base C (5/3);  total (31/48)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCAC  >  NC_000913/1534756‑1535044
                                                                                                                                                     |                                                                                                                                           
gtgATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGTAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                                             >  1:631931/1‑150 (MQ=255)
 tgATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                            >  1:11581/1‑150 (MQ=255)
  gATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                           >  2:1596511/1‑150 (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                                                                                                                           >  2:52150/1‑150 (MQ=255)
                  ggTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGcc                                                                                                                           <  1:52150/150‑1 (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                                                                                                                       >  1:1100091/1‑134 (MQ=255)
                                      gAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGt                                                                                                                       <  2:1100091/134‑1 (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:361798/150‑1 (MQ=255)
                                                 ttACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCa                                                                                            <  2:1328901/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:13712/150‑1 (MQ=255)
                                                                               cggcgTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCAtt                                                              <  1:53500/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:424748/150‑1 (MQ=255)
                                                                                      tGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATg                                                       >  1:637749/1‑150 (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:36464/1‑150 (MQ=255)
                                                                                                                              tGCCGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTcc               >  1:954800/1‑150 (MQ=255)
                                                                                                                               gCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCg              <  1:1187326/150‑1 (MQ=255)
                                                                                                                               gCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCg              <  1:962043/150‑1 (MQ=255)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACTGCTGTGCCTGAAATGGTGACTTTTTCCg              >  1:683563/1‑150 (MQ=255)
                                                                                                                                 cGGATAAGGCGTTTACGCCGCATCCGGCATACATTAGCTCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCggt            >  2:1472546/1‑150 (MQ=255)
                                                                                                                                  ggataaggcgTTTACGCCGCATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCggtt           >  2:574333/7‑150 (MQ=255)
                                                                                                                                       cggcgTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGGTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTa      >  1:496533/1‑150 (MQ=255)
                                                                                                                                         gcgTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATc    <  2:637749/150‑1 (MQ=255)
                                                                                                                                           gTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCAc  <  2:1338592/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GTGATGCAAGAACAGTTTGGTCTGGGCGTGGATGATGCGAAACATGGCTTTACCGTGGATGAGTTAGCGCTGGTGATGGCGGCGTTTGATACGCACCCGGAGGCGGGAAAATAATTTATGTCAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCAC  >  NC_000913/1534756‑1535044

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: