Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6050TC47.3% 95.7 / 24.5 75G187G (GGA→GGGycfTinner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (30/9);  new base C (0/35);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 5.87e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGT  >  NC_000913/1174458‑1174751
                                                                                                                                                      |                                                                                                                                                    
caacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTccc                                                                                                                                                   >  1:453675/1‑151 (MQ=255)
 aacCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTcccc                                                                                                                                                  >  1:130501/1‑151 (MQ=255)
     aGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacc                                                                                                                                               >  1:358264/1‑150 (MQ=255)
     aGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacc                                                                                                                                               >  1:628426/1‑150 (MQ=255)
      ggCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCcaccac                                                                                                                                             >  1:1063670/1‑151 (MQ=255)
       gCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACg                                                                                                                                            <  2:822793/151‑1 (MQ=255)
       gCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACg                                                                                                                                            <  1:347753/151‑1 (MQ=255)
        cACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTGCCCACCACgg                                                                                                                                           <  1:1013509/151‑1 (MQ=255)
         accgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGc                                                                                                                                          <  2:719953/151‑1 (MQ=255)
         accgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGc                                                                                                                                          >  1:664103/1‑151 (MQ=255)
         accgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGc                                                                                                                                          >  1:645215/1‑151 (MQ=255)
           cgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGt                                                                                                                                        <  1:950293/151‑1 (MQ=255)
             accgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTgg                                                                                                                                      <  1:881621/151‑1 (MQ=255)
                  cAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGTGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAAc                                                                                                                                  >  2:863337/1‑150 (MQ=255)
                    gAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaa                                                                                                                               <  2:851539/151‑1 (MQ=255)
                    gAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGa                                                                                                                                >  2:612394/1‑150 (MQ=255)
                     aaCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGaa                                                                                                                               >  2:874165/1‑150 (MQ=255)
                      aCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCACACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAt                                                                                                                            <  2:787140/151‑1 (MQ=255)
                            aGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCg                                                                                                                     <  2:125094/151‑1 (MQ=255)
                             gCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGc                                                                                                                    >  2:9159/1‑151 (MQ=255)
                               aaaGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGcc                                                                                                                   >  1:1117284/1‑150 (MQ=255)
                                      ccaccaGCAGATGGCGGCGCAACGGCCCCTCTTAAACACAGGTCATTATTGTTGCGCCAAACAATGCGCAAAGGCTGTAAAAAGGCAAATGGCGGATCACACTGTTCTTCCCCCCCCACGGCGTGGAAACGAAATTAACCGCCACACTca                                                                                                            <  1:190797/150‑1 (MQ=255)
                                        aTCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGCCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcagc                                                                                                          >  2:111812/1‑150 (MQ=255)
                                           agcggATGGCGGCGCACCGGCCCCTCTTAAACACAGATCTTTATTTTTCCGCAAAACCATGCGCAAAGGCTGTAAAAAGGAAATTTGGGGATCCCCCTGTTCATCCCCCCCCGCGGGGGGGAAACAAAATAAACCGCCACACTCAGCAGTa                                                                                                      <  1:546903/151‑1 (MQ=255)
                                            gcagATGGCGGCGCAACGGCCCCTCTTTAACACAGGTTATTATTGTTGCGCAAAACCAGGCGCCAAGGTTGTAAAAAGCCAATTGGCGGATCACACTTTTCATCCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTa                                                                                                      <  2:743101/150‑1 (MQ=255)
                                                  ggcggcGAAACGCCACCTTTTTAACACAGGTAATTATTGTGGCGCCAAACCATGCGCAAAGGCTGAAAAAAGGCAAATGGCGGATCACATTGTTCATCCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTa                                                                                                      <  2:1057091/144‑1 (MQ=255)
                                                        gcAAGGGCACTTTTTTAACAGAGGTTATTATTGTGGCCCAAACCCTGGCGCCACGTATGTAAAAGGGAAATTGGGGGATCCCATTGTTTATCCCCCCCCGCGGCGGGGAAGTAAATTTACCCGCCACACTCAGCAGTACAAACAAGGCg                                                                                           <  2:110933/149‑1 (MQ=255)
                                                          aaCGCCCCCTTTTAAACACAGGTTATTATTTTTCCCCCAACCCATGCCCAAGGGGTGTAAAAAGAAAAATTGGGGATCCCCCTTTTCATCCCCCCCCACGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATa                                                                                       <  1:253175/151‑1 (MQ=255)
                                                           aCGGCACCTCTTAAACACAGGTCATTTTTGTTGCGCCAACCAGTGCGCCAAGGCTAAAAGAAGGAAAAGTGCGGATCACCCTGTTATTCCCCCCCCACGGCGGGGAAAGAAAATAAGCCGCCACACTCAGCAGTACAAACAAGGCGAATa                                                                                       <  1:838878/150‑1 (MQ=255)
                                                            cGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGt                                                                                     >  2:807476/1‑151 (MQ=255)
                                                               cACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                  >  2:627421/1‑151 (MQ=255)
                                                               cACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                  >  1:1099891/1‑151 (MQ=255)
                                                                 cctcttaaaacaaGGTTATTATTGTGGCGCAAAACAAGGCGCCAGGGTTGAAAAAAGGAAAATGGCGAACCACTCTGTTCATCCCCCCCCACGGGGGGGGAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCa                                                                                 <  1:197845/138‑1 (MQ=255)
                                                                 ccTTTTAAACACGGGTCATTATTGTTGCGCAAACCCAGGCGCCAGGGCTTAAAAAAGGAAAATTGGGAACCACCCTTTTCTTCCCCCCCCACGGCGGGGAAAAGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCa                                                                                 <  1:1181736/150‑1 (MQ=255)
                                                                 ccTCTTAAACCAAGTTAATTATTGTGGCGCCAACCCAGGCGCAAAGGCTGTAATAAGGAAAATGGCGGATCACATTGTTCATCCCCCCCCACGGCGTGGGAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGt                                                                                     <  2:106968/146‑1 (MQ=255)
                                                                        aaCACGGTTCTTTTTTTTTCCGCAAAACCTGGCGCCAGGGTTATAATAAGGAAATTTGCTGATCACATTTTTCATCCCCCCCCCGGGGGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCa                                                                           <  2:711308/149‑1 (MQ=255)
                                                                        aaCACGGCTCATTATGTTTGCGAAAACCCCGGCGCCAAGGCTGTAAAAGGGAAATTACAGGATCCTGCTGTCGTTACCCCACCAGGGGGTGGAGACAAAATTAGCCGCCACCCTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGACAGa                                                                         <  1:334975/151‑1 (MQ=255)
                                                                        aaCACAGGTCATTTTTGTTGCGCAAAACCATGCGCCAAGGCTGAAAAAAGGAAAATTGCGGATCACCCTTTTCATCCCCCCCCCCGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGa                                                                         <  2:858452/151‑1 (MQ=255)
                                                                        aaCAAAGTTCATTTTTTTGGCGCAAACCAAGGCCCCAGGGTTTAAAAAAGGAAAATTGCGAATCACTCTTTTCATCCCCCCCCCCGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAg                                                                          <  2:710733/150‑1 (MQ=255)
                                                                         acacAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGAGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGa                                                                         >  1:548011/1‑150 (MQ=255)
                                                                                  attattGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaa                                                                >  1:453756/1‑150 (MQ=255)
                                                                                   tttttgtggCCCAAAACCAGGCGCCAAGGCTGTAATAAGAAAAATTGGGGATCCCCCTTTTCATCCCCCACCACGGCGGGGGAACAAAATTACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                    <  1:915627/142‑1 (MQ=255)
                                                                                      tttttgCGCAAACCCTTGCGCAAGGGCTTAAATAAGGCAATTGGCGATCCCCATTGTTCTTCCCCCCCCAGGGCGGGGGAACAAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAt                                                              <  2:390192/145‑1 (MQ=255)
                                                                                      ttgttgCCCAAAACCATGCCCCAGGGCTTAAAAAAGGAAATTTGGGGATCACCCTTTTCATCCCCCCCCCCGGCGAGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                             <  2:379405/149‑1 (MQ=255)
                                                                                         tgccgacAACCCATCCCCCAGGGCTTAAATAAGAAAAATGGCGCATCACACTGTTCATCCCCCCCCCCGGGGGGAAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaa                                                                 <  1:176822/136‑1 (MQ=255)
                                                                                             gcCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                    >  2:448327/1‑151 (MQ=255)
                                                                                              ccAACCCATGCGCCAGGGCTGTAATAGGGCAAATTGCGGATCACACTTTTCTTCCCCCCCCACGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCa                                                   <  2:586179/151‑1 (MQ=255)
                                                                                                aaaCCATGCCCCAGGGTTGAAAAAAGGCAAATGGCGGACCACCCTGTTCCTCCCCCCCCCCGGCGGGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCAc                                                  <  1:158387/150‑1 (MQ=255)
                                                                                                aaaCCATGCCCCAAGGATCAAATAGGAAAATTTGCGTATCCCTCTGTTCATCCCCCCCCGCGGCGTGGGAAAGAAAGTATCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                                 <  2:408458/151‑1 (MQ=255)
                                                                                                  aCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGa                                                >  2:408346/1‑150 (MQ=255)
                                                                                                   ccAGGCCCCAGGGCTTAAATAGGGAAATTTGCGGACCCCACTTTTCATCCCCCACCAGGGCGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCa                                                   <  1:878515/146‑1 (MQ=255)
                                                                                                    cATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                              >  1:1188765/1‑150 (MQ=255)
                                                                                                      aggcgcAAGGGCTGTAAAAAGGAAAATGGCGGATCCCACTGTTCTCCCCCCACCACGGGGTGGGAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                             <  1:136901/148‑1 (MQ=255)
                                                                                                           caaaGGCTGAAATAAGGCAAATGGCGGATCACACTGTTCATCCCCCACCCGGGGGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACa                                        <  2:418701/148‑1 (MQ=255)
                                                                                                            cAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAt                                      >  1:663800/1‑151 (MQ=255)
                                                                                                             aaGGCTGAAAAAAGCAAAATTGCGAACCACACTTTTCATCCCCCCCCACGGCGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGATGGCTAAAAATTTTACATACCACGAAATAGACaa                                       <  2:1063722/149‑1 (MQ=255)
                                                                                                               ggCTGTACTAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCGGTACAAACAAGGCGAACAGTGGCAGGGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTa                                    >  1:1089325/1‑150 (MQ=255)
                                                                                                               gcttgaaaaagggaaaaTGGCGGATCCCCCTGTTCTTCCCCCCCCCCGGGGGGGGAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTaa                                   <  2:166352/137‑1 (MQ=255)
                                                                                                                     aataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                              >  2:502390/1‑150 (MQ=255)
                                                                                                                     aataaGGAAAATTGCGGATCACACTTTTCTTCCCCCCCCACGGCGGGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCa                               <  2:845147/149‑1 (MQ=255)
                                                                                                                        aaGGAAAATTGCGGATCACCCTTTTCATCCCCCCCCCGGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAg                          <  2:191703/151‑1 (MQ=255)
                                                                                                                         aGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                            >  1:138380/1‑148 (MQ=255)
                                                                                                                          ggaaaaTTGCGGATCACACTTTTCATCCCCCACCACGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGa                         <  2:1075532/147‑1 (MQ=255)
                                                                                                                          ggaaaaTGGCGGTTCACACTGTCCATCCCCCACCCCGGCGGGGAAACAAAATTACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGa                         <  2:3147/147‑1 (MQ=255)
                                                                                                                           gCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAt                        >  2:576092/1‑150 (MQ=255)
                                                                                                                           gCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAt                        >  1:1075129/1‑150 (MQ=255)
                                                                                                                             aaaTTGCGGCTCACCCTGTTCATCCCCCCCCAGGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAc                      <  1:688309/150‑1 (MQ=255)
                                                                                                                                    gcgaataccacTTTTCTTCCCCCCCCCGGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACaa                 <  2:1193899/142‑1 (MQ=255)
                                                                                                                                       gATCACCCTTTTCATCCCCCCCCACGGCGTGGAAAAAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTg             <  1:425493/151‑1 (MQ=255)
                                                                                                                                          cacaCTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTg          >  2:861342/1‑151 (MQ=255)
                                                                                                                                               tGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTgg      >  1:608583/1‑150 (MQ=255)
                                                                                                                                                 ttCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCgg   >  1:144957/1‑151 (MQ=255)
                                                                                                                                                 ttCATCCCCCACCACGGGGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCg    <  1:996765/150‑1 (MQ=255)
                                                                                                                                                   cATCCCCCCCCACGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCg    <  2:152071/148‑1 (MQ=255)
                                                                                                                                                   cATCCCCCACCGCGGGGGGGAAACGAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt  <  2:516496/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGT  >  NC_000913/1174458‑1174751

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