Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,643,1550GA46.2% 47.5 / inf 276Y676Y (TAC→TATprlColigopeptidase A
Reads supporting (aligned to +/- strand):  ref base G (89/59);  new base A (66/61);  total (155/121)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-08

Only 100 of 280 total aligned reads displayed.

CCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGC  >  NC_000913/3643006‑3643305
                                                                                                                                                     |                                                                                                                                                      
 cAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATa                                                                                                                                                      <  4:242757/151‑1 (MQ=255)
    aTGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATg                                                                                                                                                   >  1:1753538/1‑151 (MQ=255)
    aTGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATg                                                                                                                                                   <  1:1678889/151‑1 (MQ=255)
          cGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTcc                                                                                                                                               >  3:384252/1‑149 (MQ=255)
            gCGTTAACACCAGCGCCATCAGGATGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCATGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAgc                                                                                                                                            >  3:231075/1‑150 (MQ=255)
              gTTAACCCCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAgca                                                                                                                                           <  3:41012/149‑1 (MQ=255)
              gTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAgcatc                                                                                                                                         >  4:647672/1‑151 (MQ=255)
                    acCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATcc                                                                                                                                   <  4:1604638/151‑1 (MQ=255)
                      cAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATcc                                                                                                                                   >  4:1787967/1‑149 (MQ=255)
                      cAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAg                                                                                                                                 <  4:576874/151‑1 (MQ=255)
                       aGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGc                                                                                                                                >  1:488134/1‑151 (MQ=255)
                         cgcCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTg                                                                                                                              <  1:383688/151‑1 (MQ=255)
                          gcCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGc                                                                                                                             >  4:1382085/1‑151 (MQ=255)
                           ccATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCg                                                                                                                            >  1:1195222/1‑151 (MQ=255)
                                aGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTc                                                                                                                        >  1:1910321/1‑150 (MQ=255)
                                 ggTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCAc                                                                                                                      >  4:1738530/1‑151 (MQ=255)
                                    tGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGAc                                                                                                                   >  1:857846/1‑151 (MQ=255)
                                    tGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGa                                                                                                                    <  1:225281/150‑1 (MQ=255)
                                     gTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGAcc                                                                                                                  >  4:1310829/1‑151 (MQ=255)
                                          cATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACgg                                                                                                              >  4:1141557/1‑150 (MQ=255)
                                          cATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACgg                                                                                                              >  1:842140/1‑150 (MQ=255)
                                           aTCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGa                                                                                                             >  4:545366/1‑150 (MQ=255)
                                             cGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGc                                                                                                          <  4:1540126/151‑1 (MQ=255)
                                                gCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGtt                                                                                                       >  1:1001234/1‑151 (MQ=255)
                                                      gCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACa                                                                                                 <  4:847786/151‑1 (MQ=255)
                                                        ccccAGTGGGCCGCCAGAACAGATAAGGCCCCGTTTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGa                                                                                                    >  3:305921/1‑146 (MQ=255)
                                                        ccccAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAAc                                                                                                  >  4:1068052/1‑148 (MQ=255)
                                                         cccAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGAt                                                                                              <  2:355767/151‑1 (MQ=255)
                                                         cccAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGAt                                                                                              <  1:515002/151‑1 (MQ=255)
                                                          ccAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAAc                                                                                                  >  3:470742/1‑146 (MQ=255)
                                                           cAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATcc                                                                                            >  4:1842901/1‑151 (MQ=255)
                                                            aGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATcc                                                                                            >  4:12418/1‑150 (MQ=255)
                                                              cGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATc                                                                                         >  3:278410/1‑151 (MQ=255)
                                                                  ccgccAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGct                                                                                     <  3:347777/151‑1 (MQ=255)
                                                                  ccgccAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGct                                                                                     >  4:349961/1‑151 (MQ=255)
                                                                      cAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCtcttc                                                                                 <  1:1621467/151‑1 (MQ=255)
                                                                       aGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCtcttct                                                                                <  4:1005908/151‑1 (MQ=255)
                                                                         aaCAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGa                                                                              <  4:871090/151‑1 (MQ=255)
                                                                           cAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAAc                                                                            >  2:221566/1‑151 (MQ=255)
                                                                           cAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAAc                                                                            >  2:167308/1‑151 (MQ=255)
                                                                            aGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAAcc                                                                           <  4:1423663/151‑1 (MQ=255)
                                                                             gATAAGGCACAGTCTCCGGTGTCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGTCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCg                                                                          <  2:381500/151‑1 (MQ=255)
                                                                                aaGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGcc                                                                        <  1:1479880/150‑1 (MQ=255)
                                                                                 aGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCAc                                                                      >  4:1089960/1‑151 (MQ=255)
                                                                                    cACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGc                                                                   <  4:1804250/151‑1 (MQ=255)
                                                                                    cACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGc                                                                   >  3:352226/1‑151 (MQ=255)
                                                                                      ccGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCt                                                                  <  4:1261502/150‑1 (MQ=255)
                                                                                      ccGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTc                                                                 >  2:446825/1‑151 (MQ=255)
                                                                                      ccGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTc                                                                 >  4:1863458/1‑151 (MQ=255)
                                                                                           tcCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAAt                                                            >  1:602133/1‑151 (MQ=255)
                                                                                             cGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAAtgt                                                          >  4:27119/1‑151 (MQ=255)
                                                                                              ggTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAAtgt                                                          <  4:1338043/150‑1 (MQ=255)
                                                                                               gTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAAtgttg                                                        <  1:1613378/151‑1 (MQ=255)
                                                                                                   cTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCg                                                     >  1:1879411/1‑150 (MQ=255)
                                                                                                   cTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCg                                                     >  1:1346520/1‑150 (MQ=255)
                                                                                                      tttCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAgg                                                  <  1:1305902/150‑1 (MQ=255)
                                                                                                      tttCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAgg                                                  <  2:113416/150‑1 (MQ=255)
                                                                                                       ttCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGa                                                 >  1:1057701/1‑150 (MQ=255)
                                                                                                         catcaATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAAc                                               <  1:1274215/150‑1 (MQ=255)
                                                                                                         catcaATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAAc                                               <  4:325322/150‑1 (MQ=255)
                                                                                                          atcaATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGa                                             <  1:1122766/151‑1 (MQ=255)
                                                                                                            caATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGAc                                            <  3:34045/150‑1 (MQ=255)
                                                                                                            caATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGAc                                            >  1:62748/1‑150 (MQ=255)
                                                                                                            caATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACt                                           >  1:881583/1‑151 (MQ=255)
                                                                                                             aaTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTg                                          >  4:1620469/1‑151 (MQ=255)
                                                                                                              aTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGc                                         >  4:1476684/1‑151 (MQ=255)
                                                                                                               ttAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGcc                                        >  4:739785/1‑151 (MQ=255)
                                                                                                                tAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGcc                                        <  4:776657/150‑1 (MQ=255)
                                                                                                                    cAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGtt                                   >  2:496486/1‑151 (MQ=255)
                                                                                                                    cAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGtt                                   <  1:38354/151‑1 (MQ=255)
                                                                                                                     aGATTTTCCCTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGATGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCAGGttt                                  >  4:1251095/1‑151 (MQ=255)
                                                                                                                     aGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGttt                                  >  1:272281/1‑151 (MQ=255)
                                                                                                                       aTTTTCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTc                                 >  4:923268/1‑150 (MQ=255)
                                                                                                                        ttttCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCa                                <  4:705701/150‑1 (MQ=255)
                                                                                                                           tCACTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGt                            >  4:397480/1‑151 (MQ=255)
                                                                                                                             aCTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTg                          >  1:303454/1‑151 (MQ=255)
                                                                                                                              cTGAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTg                          <  1:1840635/150‑1 (MQ=255)
                                                                                                                                gAATGATCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGaaa                       <  1:423531/151‑1 (MQ=255)
                                                                                                                                  aTGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAAt                     >  1:750977/1‑151 (MQ=255)
                                                                                                                                   tGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATg                    >  1:93260/1‑151 (MQ=255)
                                                                                                                                    gATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGc                   >  2:323777/1‑151 (MQ=255)
                                                                                                                                     aTCAGCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGcc                  <  4:1109694/151‑1 (MQ=255)
                                                                                                                                      tCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGccc                 >  2:13209/1‑151 (MQ=255)
                                                                                                                                      tCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGccc                 >  4:1530594/1‑151 (MQ=255)
                                                                                                                                         gCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCctct              >  4:1301902/1‑151 (MQ=255)
                                                                                                                                         gCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCctct              >  4:500587/1‑151 (MQ=255)
                                                                                                                                         gCCCTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCctct              <  4:1186735/151‑1 (MQ=255)
                                                                                                                                          cccTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCtt             >  4:1058877/1‑151 (MQ=255)
                                                                                                                                           ccTTAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTc            >  2:425104/1‑151 (MQ=255)
                                                                                                                                            cTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTcc           >  1:707715/1‑151 (MQ=255)
                                                                                                                                             ttAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCt          >  2:202132/1‑151 (MQ=255)
                                                                                                                                              tAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTc         >  3:509293/1‑151 (MQ=255)
                                                                                                                                              tAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCt          <  3:91801/150‑1 (MQ=255)
                                                                                                                                              tAATGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTc         >  3:170571/1‑151 (MQ=255)
                                                                                                                                               aaTGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTc         <  1:1829181/150‑1 (MQ=255)
                                                                                                                                               aaTGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCa        <  1:1404327/151‑1 (MQ=255)
                                                                                                                                                aTGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCa        >  3:157016/1‑150 (MQ=255)
                                                                                                                                                aTGCCATAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCa        >  3:459532/1‑150 (MQ=255)
                                                                                                                                                   ccGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAgc    <  3:194523/151‑1 (MQ=255)
                                                                                                                                                     gTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAgcgc  <  4:1925490/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCAGATGTTCCGGCGTTAACACCAGCGCCATCAGGTTGTCTTCATCGTGCTCCAGCCCCCAGCGGGCCGCCAGAACAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGC  >  NC_000913/3643006‑3643305

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