Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,291,110 A→T 5.1% K766I (AAA→ATA)  yraJ → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,291,1100AT 5.1% 1356.8 / 45.6 373K766I (AAA→ATA) yraJputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (165/188);  new base T (16/3);  total (182/191)
Fisher's exact test for biased strand distribution p-value = 1.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 376 total aligned reads displayed.

TCGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAAT  >  NC_000913/3290962‑3291256
                                                                                                                                                    |                                                                                                                                                  
tcGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTaaa                                                                                                                                                   >  1:83165/1‑150 (MQ=255)
   gTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCg                                                                                                                                                >  1:3753151/1‑150 (MQ=255)
   gTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCg                                                                                                                                                <  1:2214907/150‑1 (MQ=255)
    tAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGt                                                                                                                                               >  1:909841/1‑150 (MQ=255)
      gATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCACCGGTAAACCGTTg                                                                                                                                             >  2:1109030/1‑150 (MQ=255)
        tactacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGcc                                                                                                                                           >  1:2382693/1‑150 (MQ=255)
        tactacGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATAGCAACGGTAAACCGTTGcc                                                                                                                                           >  2:3055048/1‑150 (MQ=255)
            acGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCttt                                                                                                                                       <  1:3436128/150‑1 (MQ=255)
             cGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTg                                                                                                                                      <  1:1692811/150‑1 (MQ=255)
                     cccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCtt                                                                                                                              >  1:1121917/1‑150 (MQ=255)
                     cccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCtt                                                                                                                              >  2:1502002/1‑150 (MQ=255)
                      ccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTg                                                                                                                             <  2:1844747/150‑1 (MQ=255)
                      ccGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTg                                                                                                                             >  1:625648/1‑150 (MQ=255)
                       cGATAACGTCGATCTTGAACAAACAACCCAGCTTGTGGTGCCCAACAGAGGTGCAATGGGAGCGGCGCGTTTCAACGCCAAGATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGc                                                                                                                            <  2:2936959/150‑1 (MQ=255)
                        gATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGcc                                                                                                                           >  1:3697461/1‑150 (MQ=255)
                         aTAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCa                                                                                                                          <  2:4241590/150‑1 (MQ=255)
                         aTAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCa                                                                                                                          <  1:1471369/150‑1 (MQ=255)
                             cGTCGATCTTGCCCAAACCCCCCAGGGTGTGGTGCCCACCAGAGCTGCAATGGTAGCGGCGCGTTTCACCGCCAATATCGGTTATCGCGTAATTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCaa                                                                                                                      <  2:2519425/150‑1 (MQ=255)
                              gTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAAc                                                                                                                     >  1:1240302/1‑150 (MQ=255)
                                cGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACga                                                                                                                   <  2:4152848/150‑1 (MQ=255)
                                 gATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgat                                                                                                                  <  1:4284446/150‑1 (MQ=255)
                                  aTCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgatg                                                                                                                 >  1:2821961/1‑150 (MQ=255)
                                   tCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgatga                                                                                                                >  1:3635828/1‑150 (MQ=255)
                                    cTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgatgat                                                                                                               >  1:1932797/1‑150 (MQ=255)
                                    cTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACgatgat                                                                                                               >  1:590945/1‑150 (MQ=255)
                                     ttGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATa                                                                                                              >  2:1685484/1‑150 (MQ=255)
                                      tGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATAc                                                                                                             >  2:874258/1‑150 (MQ=255)
                                       gAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACg                                                                                                            >  2:1977814/1‑150 (MQ=255)
                                       gAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACg                                                                                                            <  1:2548944/150‑1 (MQ=255)
                                         acaaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACggg                                                                                                          >  2:3582376/1‑150 (MQ=255)
                                          caaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgggg                                                                                                         <  2:1737422/150‑1 (MQ=255)
                                          caaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgggg                                                                                                         >  2:1482530/1‑150 (MQ=255)
                                          caaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTGTGGCGCTCTTGCCAGCAACGATGATACgggg                                                                                                         >  2:3317725/1‑150 (MQ=255)
                                          caaacaaCACAGTGTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACgggg                                                                                                         <  2:4190169/150‑1 (MQ=255)
                                           aaacaaCACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCCGCGATCGCAACGGTAAACCGTTGCCCTGTGGCGCTCTTGCCAGCACCGCTGATACGGGGc                                                                                                        >  2:3465214/1‑150 (MQ=255)
                                              caacaCAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGcaac                                                                                                     <  1:1595427/150‑1 (MQ=255)
                                              caacaCAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGcaac                                                                                                     <  1:1788378/150‑1 (MQ=255)
                                               aacaCAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGcaaca                                                                                                    >  2:1316720/1‑150 (MQ=255)
                                                 cacaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACaaa                                                                                                  <  1:546538/150‑1 (MQ=255)
                                                  acaGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAg                                                                                                 >  1:3083673/1‑150 (MQ=255)
                                                     gTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAt                                                                                              <  2:2829571/150‑1 (MQ=255)
                                                     gTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAt                                                                                              >  1:2019517/1‑150 (MQ=255)
                                                     gTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAt                                                                                              >  1:2197073/1‑150 (MQ=255)
                                                      tttGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTAtc                                                                                             >  2:3630508/1‑150 (MQ=255)
                                                      tttGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCGTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACCAAGTAtc                                                                                             >  2:247391/1‑150 (MQ=255)
                                                           ggtgCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGa                                                                                        <  1:1341654/150‑1 (MQ=255)
                                                            gtgCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGAt                                                                                       >  1:2756041/1‑150 (MQ=255)
                                                              gCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGa                                                                                     <  1:226937/150‑1 (MQ=255)
                                                                ccAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAgg                                                                                   <  2:3927923/150‑1 (MQ=255)
                                                                  aaCAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGggc                                                                                 <  1:4129798/150‑1 (MQ=255)
                                                                    cAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGggcgg                                                                               <  1:782462/150‑1 (MQ=255)
                                                                     agagGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGggcggc                                                                              >  2:3749432/1‑150 (MQ=255)
                                                                      gagGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCa                                                                             <  2:1222153/150‑1 (MQ=255)
                                                                           gCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACta                                                                        <  2:1185130/150‑1 (MQ=255)
                                                                           gCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACta                                                                        >  1:4020617/1‑150 (MQ=255)
                                                                            cAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACtat                                                                       >  1:2053704/1‑150 (MQ=255)
                                                                            cAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACtat                                                                       >  2:28991/1‑150 (MQ=255)
                                                                             aaTGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACtata                                                                      <  1:3964451/150‑1 (MQ=255)
                                                                             aaTGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACtata                                                                      <  1:1507064/150‑1 (MQ=255)
                                                                               tGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATAtc                                                                    >  2:3165409/1‑150 (MQ=255)
                                                                                     cggcgCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTg                                                              >  2:3239126/1‑150 (MQ=255)
                                                                                      ggcgCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTgg                                                             <  2:3373789/150‑1 (MQ=255)
                                                                                      ggcgCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTgg                                                             <  2:2243125/150‑1 (MQ=255)
                                                                                        cgcgTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGa                                                           <  1:1178623/150‑1 (MQ=255)
                                                                                              tCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGa                                                     <  2:4328673/150‑1 (MQ=255)
                                                                                              tCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGa                                                     <  1:392041/150‑1 (MQ=255)
                                                                                                   gCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATATCTCTGGGATATCGAGTaaa                                                >  1:100719/1‑150 (MQ=255)
                                                                                                     cAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTACATCTCTCTGGGATATCGAGTAAATc                                              >  2:1900964/1‑150 (MQ=255)
                                                                                                      aaTATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATca                                             <  1:3055123/150‑1 (MQ=255)
                                                                                                          tCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACaaa                                         <  2:602471/150‑1 (MQ=255)
                                                                                                          tCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACaaa                                         <  2:2558099/150‑1 (MQ=255)
                                                                                                              ttATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTg                                     <  1:1309775/150‑1 (MQ=255)
                                                                                                                aTCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGa                                   <  2:3522535/150‑1 (MQ=255)
                                                                                                                  cgcgTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACt                                 <  1:1482530/150‑1 (MQ=255)
                                                                                                                    cgTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCACCAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATCTCGAGTAAATCACAAAGCTGGACTGt                               >  1:2921048/1‑150 (MQ=255)
                                                                                                                    cgTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGt                               >  1:3274066/1‑150 (MQ=255)
                                                                                                                         ttGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCt                          >  1:3777258/1‑150 (MQ=255)
                                                                                                                           gTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTgg                        <  1:2865508/150‑1 (MQ=255)
                                                                                                                            ttACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                       <  2:2084604/150‑1 (MQ=255)
                                                                                                                            ttACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                       >  1:1036798/1‑150 (MQ=255)
                                                                                                                            ttACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTggg                       <  2:1029416/150‑1 (MQ=255)
                                                                                                                             tACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTgggg                      <  1:3817287/150‑1 (MQ=255)
                                                                                                                                aGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGaaa                   >  2:3796062/1‑150 (MQ=255)
                                                                                                                                  tCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGACATc                 >  1:2651563/1‑150 (MQ=255)
                                                                                                                                    aGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAg               <  2:766955/150‑1 (MQ=255)
                                                                                                                                     gCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAgg              <  1:949949/150‑1 (MQ=255)
                                                                                                                                       gATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCa            <  2:2545995/150‑1 (MQ=255)
                                                                                                                                       gATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCa            <  1:3579027/150‑1 (MQ=255)
                                                                                                                                        aTCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAg           >  2:663523/1‑150 (MQ=255)
                                                                                                                                          ggCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGAt         >  2:2732569/2‑150 (MQ=255)
                                                                                                                                          ggCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGAt         >  2:3884612/2‑150 (MQ=255)
                                                                                                                                          ggCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTAAATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGACATCAGGCAGAt         >  2:1809110/2‑150 (MQ=255)
                                                                                                                                          ggCAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGAGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGAt         >  2:654107/2‑150 (MQ=255)
                                                                                                                                           gCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATc        >  1:4348261/1‑150 (MQ=255)
                                                                                                                                           gCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATc        >  2:1555211/1‑150 (MQ=255)
                                                                                                                                            cAACGGTTTACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATca       >  2:1259687/1‑150 (MQ=255)
                                                                                                                                            cAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATca       >  1:1271271/1‑150 (MQ=255)
                                                                                                                                              aCGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGTCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATcaac     <  2:3818267/150‑1 (MQ=255)
                                                                                                                                              aCGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATcaac     <  1:2516058/150‑1 (MQ=255)
                                                                                                                                              aCGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATcaac     <  1:2115709/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGGTAGATACTACGCAGCTGCCCGATAACGTCGATCTTGAACAAACAACACAGTTTGTGGTGCCCAACAGAGGTGCAATGGTAGCGGCGCGTTTCAACGCCAATATCGGTTATCGCGTACTTGTTACAGTCAGCGATCGCAACGGTAAACCGTTGCCCTTTGGCGCTCTTGCCAGCAACGATGATACGGGGCAACAAAGTATCGTCGATGAGGGCGGCATACTATATCTCTCTGGGATATCGAGTAAATCACAAAGCTGGACTGTACGCTGGGGAAATCAGGCAGATCAACAAT  >  NC_000913/3290962‑3291256

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