Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6050TC32.6% 100.3 / 14.0 50G187G (GGA→GGGycfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (31/2);  new base C (0/16);  total (31/19)
Fisher's exact test for biased strand distribution p-value = 4.57e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCAC  >  NC_000913/1174375‑1174836
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                       
aTACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacca                                                                                                                                                                                                                                 <  1:1178352/239‑1 (MQ=255)
 tACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCcaccac                                                                                                                                                                                                                                >  1:1036576/1‑239 (MQ=255)
                                 ttCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGGAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCACCCCAATca                                                                                                                                                                                                 >  2:923010/1‑238 (MQ=255)
                                     taaccaccaaAACCGACACAAGCAAAAACAGCAGCGAAATAAAAACACACAGGCACCCGCCCGCAAGAGCGGAAAGAAAAGAAGCCACCAGGAGAGGGCGGCGCAAGGGCCCCTTTTAAACACGGGTCATTATTGGTGCGCAAAACCAGGCGCAAAGGCTGTAAAAAGGAAAATGGAGGACCACCCTGTTCCTCCCCCCCCACGGGGTGGGAACAAAATTAACCGCCACACTCAGCAGt                                                                                                                                                                                            <  2:112633/235‑1 (MQ=255)
                                          accacAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGTCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaa                                                                                                                                                                                        >  1:1175243/1‑238 (MQ=255)
                                                 tCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcagca                                                                                                                                                                                              >  1:1103295/1‑225 (MQ=255)
                                                    aTACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAACAAAGGCGAAt                                                                                                                                                                             >  2:926887/1‑239 (MQ=255)
                                                      aCCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCAAACTCAGCAGTACAAACAAGGCGAATAg                                                                                                                                                                           >  2:244534/1‑239 (MQ=255)
                                                      aCCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCCCCACACTCCGCCGTACAAACAAGGCGAagag                                                                                                                                                                           >  2:233494/1‑239 (MQ=255)
                                                                         agaTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTa                                                                                                                                                        >  1:440726/1‑239 (MQ=255)
                                                                                      ccAGGCACCGACAGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCCCGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGAAGCGACAGACGGCCaaaaacataccaa                                                                                                                                            >  2:629265/1‑230 (MQ=255)
                                                                                         ggCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAAAcc                                                                                                                                         >  2:475042/1‑238 (MQ=255)
                                                                                                 accgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCCACCGCACCCCCTTAACACAGGTCATTATTGCTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGCTCACACTGTTCATTCCCCACCACGGCGTGGGCACGAAATTAACCGGCCCACTCAGCAGCACAACAACTGCGAACAGAGGCAGCGCCAAACCGCGAAAACACTTACAAACCaacaaaaa                                                                                                                                 >  2:516451/1‑232 (MQ=255)
                                                                                                        aaccGGCAAACAAAGACGCCATCAGCGTGGGGCGGCGCACCGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAACCCATGCGCCAGGGCTGAAAAAAGAAAAATTGCGGATCACCCTGTTCATCCCCCACCACGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAtt                                                                                                                          <  2:440726/236‑1 (MQ=255)
                                                                                                           aggTAAGCAAAAGAACCCATGAGAAGGGGGCGGCGCACCGGCCCCTCTAAACCCCAGTTCTTTTTTGTGGCGCAAACCCGTGCGCCAGGGGTAAAAAAAGAAAATTGGCGACTCCCATTGTCCTTCCCCCCCCCCGGGGGGGAAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTaa                                                                                                                        <  2:378380/238‑1 (MQ=255)
                                                                                                             gtaaaaaaaaGAAGCCATCGGCAGATGGCGACACAGCGGCACTTCTTAAAAAGAGGAAATTTTGGTGCGCCAAACCAACGCACAAAGGGTAAAAAAAGGAAAAATGCGGAACACACTGTTCAGCCCCCCCCAGGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACg                                                                                                                      <  2:1175243/232‑1 (MQ=255)
                                                                                                                      agaagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACAACGAAATAGACAATTAACGCATACAGAt                                                                                                             >  1:1085706/1‑239 (MQ=255)
                                                                                                                                   gtggggcggcGCACCGACACTCTGTTAACACAGGTAATTGTTGGTACGCAAACCCATGCGCAAAGGCTGAAATAGGACAAATGGCGGATCACACTGTTCCTCCCCACCCAGGGGGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCt                                                                                                   <  2:938544/232‑1 (MQ=255)
                                                                                                                                     ctGGCGGCGCAACGACATCTCGTAAACCCATGTCATGATTTACGCGACAGACCGCGCTCCAGGGTTGTAATAAGGAAAATGCGGAACCCCCGTTTTCTCCCCCCCCCACGGGGGGGAAAAAAAATTAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATCCCACAAGCTGGTg                                                                                               <  2:112/237‑1 (MQ=255)
                                                                                                                                                       atctTCACCACAGGTAATTGTTGTGCCGCCAAACCAGGCACCAAGGTTATAAGAAGCAAAATGGCGCACCACATTTTTCTCCCCCCCCCGCGGGGGGGAAAAGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTg                                                                                               <  2:671970/219‑1 (MQ=255)
                                                                                                                                                          tttAACACGGGTATTTTTTGTTGCCCCAACCCAGGCGCCAAGGGTGTAATAAGGCAAATTGCGGATCACCCTGTTCATCCCCCACCCCGGCGTGGAAACGAAATTAACCg                                                                                                                                                                                                          <  1:778405/110‑1 (MQ=255)
                                                                                                                                                                  caGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACa                                                                                                                >  1:1321983/1‑193 (MQ=255)
                                                                                                                                                                     gTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTccc                                                                                                                                                                                                                                      <  1:1333946/72‑1 (MQ=255)
                                                                                                                                                                     gTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTccc                                                                                                                                                                                                                                      >  2:1333946/1‑72 (MQ=255)
                                                                                                                                                                     gTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcag                                                                                                                                                                                                >  2:778405/1‑110 (MQ=255)
                                                                                                                                                                     gTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaac                                                                                                                                                                                      >  1:102388/1‑120 (MQ=255)
                                                                                                                                                                     gTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTa                                                                                                                                                                                           >  1:303077/1‑115 (MQ=255)
                                                                                                                                                                     gTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGTGGGGGTCAACCCCGGCAGGAAcacccc                                                               >  2:867623/1‑234 (MQ=255)
                                                                                                                                                                     gTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGGGat                                                                                    >  1:27405/1‑218 (MQ=255)
                                                                                                                                                                     gTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGGGCTGGCGGGGATCATCCCCTGCAGGCAATCAcc                                                               >  2:376775/1‑239 (MQ=255)
                                                                                                                                                                          tattGTTGCGCCAACCCATGCGCCAAGGCTGTAATAAGGAAAATGGCGGATCACACTGTTCATCCCCCCCCCCGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCa                                                                 <  1:64407/232‑1 (MQ=255)
                                                                                                                                                                               ttgCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTa                                                                                                                                                                                           >  1:1261673/1‑105 (MQ=255)
                                                                                                                                                                                 gcgcCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATa                                                                                                            >  1:1404514/1‑182 (MQ=255)
                                                                                                                                                                                    ccAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGAGCAATCCCCGCAGGAACACACCCGTgggaacccgt                                                  >  2:64407/1‑229 (MQ=255)
                                                                                                                                                                                             tGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACt                                                                                                                                                                                                   >  1:917594/1‑84 (MQ=255)
                                                                                                                                                                                              gcgcCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGAGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAAAAAGATACCACAAGCTGGTGCTGGGGGGGAAAATCCCCTGCAGGAAACCACCGGGGGAAAAGGCAAAACCGGCAtt                                       >  1:228374/1‑239 (MQ=255)
                                                                                                                                                                                                cccAAGGCTGAAATAAGGAAAATTGCGGCTCCCCCTTTTCATCCCCCCCCCCGGCGGGGAAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAgg                                                                                                                                                                                  <  1:473332/97‑1 (MQ=255)
                                                                                                                                                                                                     ggCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGa                                                                                                                               >  1:340101/1‑144 (MQ=255)
                                                                                                                                                                                                       gtggaataaGGAAAAGGGGGGCTCACATTGTTCATCCCCCCCCGCGGAGTGAGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCa                                <  2:53570/233‑1 (MQ=255)
                                                                                                                                                                                                          taataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTgg                                                                                                                                                                        >  2:473332/1‑98 (MQ=255)
                                                                                                                                                                                                          taataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGGGGTCAACCCCTGCAGGGACCCCACGGTGGAAACCGCCTAAGCGGCATTGGCGGCCTTGCt                           >  1:53570/1‑239 (MQ=255)
                                                                                                                                                                                                          taataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGAAGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCCCCGGGGGAAACGGGAAAAGCGGCATTGGAgg                                  >  1:1307926/1‑232 (MQ=255)
                                                                                                                                                                                                            ataaGGCAAATTGCGGATCACCCTGTTCATCCCCCACCACGGGGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACca                                                                                                         <  2:114182/159‑1 (MQ=255)
                                                                                                                                                                                                                     attggagATCACCCTGTGCATCCCCCCCCAGGGGGGGGAAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCt                <  1:244534/233‑1 (MQ=255)
                                                                                                                                                                                                                      ttGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTg                                                                                               >  1:114182/1‑159 (MQ=255)
                                                                                                                                                                                                                                 cccttttCATCCCCCCCCACGGGGGGGGAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGc                                                   <  1:664666/188‑1 (MQ=255)
                                                                                                                                                                                                                                  attGTTCATACCCCCCCGGGGGGTGGGAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgccagcca   <  1:923010/237‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGGGatcatc                                                                                >  1:1135437/1‑159 (MQ=255)
                                                                                                                                                                                                                                      ttCATCCCCCCCCGCGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCAc  <  1:324282/237‑1 (MQ=255)
                                                                                                                                                                                                                                        cATCCCCCCCCAGGGGGTGGAAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAg                                                                                                                                <  2:607386/109‑1 (MQ=255)
                                                                                                                                                                                                                                           tCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACa                                                                                                                >  1:86934/1‑122 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                       
ATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCAC  >  NC_000913/1174375‑1174836

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