Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,8830CT48.9% 74.5 / 120.0 131intergenic (+52/‑51)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base C (16/51);  new base T (64/0);  total (80/51)
Fisher's exact test for biased strand distribution p-value = 1.49e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

TACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACG  >  NC_000913/2386736‑2387032
                                                                                                                                                       |                                                                                                                                                                   
tACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCt                                                                                                                                                                  <  1:1065716/151‑1 (MQ=255)
    aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                                                                                             <  1:661988/151‑1 (MQ=255)
         cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGttt                                                                                                                                                     <  1:1696234/151‑1 (MQ=255)
           cGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCAAATTTTTTTTTTCCCCCTCttttttt                                                                                                                                                   >  2:1480494/1‑150 (MQ=255)
             gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGGAAAAAAGAAAAAGACAATTTATTTTTTCTCATCTGtttttta                                                                                                                                                 >  1:1047384/1‑148 (MQ=255)
             gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCCGAAAAAAAGAAAAGCCAAATTTTTTTTTTCTCCTCTGTTTTTAg                                                                                                                                                >  1:522680/1‑151 (MQ=255)
             gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTTCGtttttaa                                                                                                                                                 >  1:382935/1‑148 (MQ=255)
                        gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGGACATTTTTTTTTTCCTCATCtttttttaaaagggaaaaat                                                                                                                                      >  2:1371814/1‑137 (MQ=255)
                          cagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGGACATTTTTTTTTTCCTCCTCTTTTTTTAAATACGAAACAAGa                                                                                                                                >  2:614859/1‑151 (MQ=255)
                          cagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAAGACATGCACATTTTATTTTTTCTCTTCTGTTTTTAGATTCGAAAAATGa                                                                                                                                >  1:854863/1‑151 (MQ=255)
                            gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAAAAAAACACATTTTTTTTTTTCTCATCTGTTTTTAGAATCGAAAAAAGaaa                                                                                                                              >  1:532742/1‑151 (MQ=255)
                            gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAATCAAAATTTATTTTTCCTCCTCTGTTTTTAAATggagaaaaataaa                                                                                                                                >  2:1703947/1‑139 (MQ=255)
                            gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAACACATTTTTTTTTTTCCCCTCTTTTTTTAGATGCGgaaaaaaaaa                                                                                                                                >  2:50733/1‑144 (MQ=255)
                             aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGAAAAGCACCTTTTTTTTTTCCCCCTCTTTTTTTAGATGCGgaaaaaaaaat                                                                                                                               >  2:728087/1‑143 (MQ=255)
                                          aaaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaa                                                                                                             <  1:50716/150‑1 (MQ=255)
                                                gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGCAAAGACCATTTTTTTTTTTCTCATCTGTTTTTAGATACGGAAACTGAAATTCAAACGGGGGaaaaaaaa                                                                                                      >  1:1408933/1‑151 (MQ=255)
                                                     ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCAAAAAAAAGACAAAGACAATTTTTTTTTTCTCCTCTGTTTTTAAATACGAAAAAAGAAATTTAACAGAGGGAAAAAAAAATggg                                                                                                 >  1:1313303/1‑151 (MQ=255)
                                                     ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGAAAAGCACATTTTTTTTTTTCTCATCTGTTTTTATATGCGGAAAAAGAAAATTAAAGGAGGGAAAAAAAATtggg                                                                                                >  1:452300/1‑150 (MQ=255)
                                                     ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACATGCCAAATTTTTTTTTCCTCCTCTTTTTTTATAAAGGGAAAAAGAAAATTAAAGGAGGGAAAAAAAAAgggg                                                                                                >  2:260239/1‑150 (MQ=255)
                                                     ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGAAATGCCCATTTTTTTTTTCCCCCTCTTTTTTTAAAAGCGAAAAATGAAAATCAACGGAGGGAAAAAAAATGGGa                                                                                                >  2:683118/1‑151 (MQ=255)
                                                     ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGACAATGACATTTTTTTTTTTCTTCTCTGGTTTTAAATGCGAAACAAGAAAATCAAAGGGGGGAAAATAAATgggg                                                                                                >  1:82968/1‑150 (MQ=255)
                                                      gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGACAAATCCCTCTTTTTTTTCCTCCTTCGTTTTAACAAGCGGAAAAAAGAATTCCAAGGAGGAAAAATCAATtggga                                                                                               >  2:929615/1‑149 (MQ=255)
                                                         gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGCAAACAACAATTTTTTTTTCCTCCTCTGGTTTTAGAAGGGGAAAAAGAAAAATAAAAGAGGAaaaaaaaaatggggaga                                                                                            >  2:672416/1‑142 (MQ=255)
                                                          tGGATGGTGTTATAACTTTGTTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATg                                                                                             <  2:117370/150‑1 (MQ=255)
                                                          tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGa                                                                                            <  1:1479813/150‑1 (MQ=255)
                                                          tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAt                                                                                          <  2:279775/151‑1 (MQ=255)
                                                            gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAtt                                                                                         <  2:1418877/150‑1 (MQ=255)
                                                             aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTg                                                                                        <  2:437391/150‑1 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGAAAAAAACATTTTTTTTTTCCTCCTCTTTTTTTAAATGCGAAAAAATAAATTTAACGGGGGAAAAAAAAAAAGGGAAGATTTGCTa                                                                                    >  2:1277275/1‑151 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCCAAAAAAAAAAAAACACATTTTTTTTTTTCTCCTCTGTTTTTAAAAACGGAAAATGAAAATTTAACGGGGAAAAAAAAAAGGGGATGAtttggca                                                                                    >  1:760449/1‑146 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGAACAGGACATTTTTTTTTTCCTCCTTTTTTTTTAAAATGGAAAAAAAAAAATCAACGGAGGGAAAATAAAATGGGAAGGTTTggtt                                                                                     >  2:776705/1‑149 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAATTACATTTTATTTTTCCTCATCTGTTTTTAAATGCGCAAAATGAAAATTAAAGGAGGAAAAATAAATGGGAAAGAtttggta                                                                                    >  2:1529155/1‑146 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCAAAAAAAAGACATGCACATTTTTTTTTTTCTCCTCTTTTTTTAAATGCGAAAAATGAAAATTAAAGGAGGGAAAAAAAATAGGAAATATTGGGTa                                                                                    >  1:1178027/1‑151 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt                                                                                     <  2:465404/150‑1 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACAATTTATTTTTCCTCCTTTGTTTTTAGATGCGAAAAAAGAAAATTCACGGGGGGAAAAAAAATTGGGATGAtttggta                                                                                    >  1:370815/1‑146 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGAAAAGCACATTTTATTTTTCCTCATCTTTTTTTAAAAGCGAAAAAAAAAAATTAACGGGGGGAAAAAAAAATGGGATGATTTGCtt                                                                                    >  2:197227/1‑150 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCCCCTCTGTTTTTATAATCGGAAAAAGAAATTCAAAGGAGGAAAAAAAAAAGGGGAAGATTgggtt                                                                                    >  2:242336/1‑148 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACAATTTTTTTTTCCTCCTCTGTTTTTAAAAGCGAAAAAAGAAAATTAAAGGAGGGAAAAAAAAAGGGGAAGATTTGCTa                                                                                    >  2:536547/1‑151 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAAAAAAACAAATTTTTTTTTTCCCCCTCTGTTTTTAAATGGGAAAAATAAAAATTAAAGGGGGAAAAAAAAAATGGGaggatttggtt                                                                                      >  2:1666113/1‑141 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGACAAGCAAATTTTATTTTTCCTCCACTTTGTTTAGATGAGGAACATGGAAATTAATCGAGGGAAaaaaaaatagagtgcattggca                                                                                      >  1:1189736/1‑137 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAATTCACAATTTTTTTTTCCCTCTTTGTTTTAAGAATCGAAAAATGAAAATTAAAGGGGGGAAAAAAAAATGGGatgaatgggtt                                                                                      >  1:1257792/1‑141 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAGCAAATTTTATTTTTTCTCATCTGGTTTTAAAAGGGGAAAAATAAAATTAAAAGAGGAAAAATAAATGGGGATGATTggggt                                                                                    >  1:991024/1‑148 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGGACATTTTATTTTTCCTCCTCTGTTTTTAGAAGCGGAAAATGAAAATCAAAGGAGGAAAAATAAAATGGGATGAtttggta                                                                                    >  1:53815/1‑146 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGGACATTTTATTTTTCCTCATTTGTTTTTAGAAGGGGAAAAAGAAAATTTAAGGGGGGAAAAAAAATTGGGAAGAATGGCTa                                                                                    >  2:807671/1‑151 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCCCATTTTTTTTTTCCCCATCTGTTTTTAGATGCGGAAAAAGAAAATTAAAGGGGGGAAAAAAAATGGGGATGAtttggta                                                                                    >  2:1024842/1‑146 (MQ=255)
                                                                gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAGAAGCGGAACATGAAATTCAACGGGGGGAAAAAAAATGGGGATGATTGGCTa                                                                                    >  2:110580/1‑151 (MQ=255)
                                                                 tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa                                                                                    <  2:1178425/150‑1 (MQ=255)
                                                                  gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGACAAACCAAATTTATTTTTTCTCCTCTGTTTTTAAATAAGAAACATGAAATTCAACCGAGGGAAAAAAAATTGTGATGATTGGGTAtt                                                                                  >  1:1182684/1‑151 (MQ=255)
                                                                    tataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTAtttt                                                                                <  2:625458/151‑1 (MQ=255)
                                                                           ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAACGAAATGCACATTTTTTTTTTTCTTATTTGTTTTTATAATCGAAAAAAGAAAATTAACTGAGGGAAAAAAAATTGGAAAGATTTGCTATTTTTCAaaaa                                                                         >  1:666775/1‑151 (MQ=255)
                                                                             gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAGGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAAAAGGGGAAAAAAAAAATTAAAGGAGGGAAAATAAATTGGGATGATTGGGttttttttcaaaaat                                                                       >  2:1688894/1‑142 (MQ=255)
                                                                               ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACCATAAGTCCATTTTTTTTTTTCCTCTCCGTTTTTAGAGTTGCAAAACAAAATTTTAACGAGGGAAAATAAATGGGGATGGATGGGTTTTTTTCAGaaaatt                                                                      >  1:1143926/1‑148 (MQ=255)
                                                                               ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACAATTTTTTTTTCCTCATCTGTTTTTAAAAGGGAAAAATGAAAATTAACGGGGGGAAAAAAAATAGGGATGATTGGCTAttttttaaaaaatt                                                                      >  2:1129124/1‑141 (MQ=255)
                                                                                 ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATTATTGGCTATTTTGCAGAAATTGAc                                                                    <  1:612228/151‑1 (MQ=255)
                                                                                   acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc                                                                  <  1:545213/151‑1 (MQ=255)
                                                                                    caacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCCCCTTTTTTTTTTCCCCCTCCTTTTTTAAATGCGAAAAATGAAATTCAACGGGGGAAAAATAAAATGGGATGGTTGGGTATTTTTGAGaaaattaatca                                                                 >  2:45600/1‑143 (MQ=255)
                                                                                    caacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTTCCTCCTCTGTTTTAAGATGCGAAACAAGAAATTCAAAGGAGGAAAAATAAATGGGAATGTTTTGCTATTTTTTAGAAATTAAccta                                                                 >  2:414174/1‑148 (MQ=255)
                                                                                    caacaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAACACATTTTATTTTTTCCCTTCTTGTTTTAAAAGCGGAAAAAGAAAATCTAAGGAGGGAAAAAAAAAAGGGAAGAATTGGTAttttttcgaaaatgaacta                                                                 >  1:1166955/1‑136 (MQ=255)
                                                                                     aacaGTAAATAATAAAAATTTCTCGAGGCGGGGGCAAAAAAAAGGAAAAGCACATTTTATTTTTCCTCCTCTGTTTTTAGATGCGGAACAAGAAAATTAAACGAGGGGAAAAAAATAGGGATGATTTGGTATTTTTTAGAAAATTACTCAg                                                                >  1:1613238/1‑151 (MQ=255)
                                                                                      acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACACGCACATTTTTTTTTTCCCCCCCTTTTTTTAGATGGGAAAAAAAAAAATTAAAGGAGGAAAAAAAAAATGGAAAGATTTGGTTTTTTTCAGAAAATTACTCaaa                                                               >  2:1535162/1‑149 (MQ=255)
                                                                                      acaGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAGAACAATTTATTTTTCCCCATCTGTTTTTAGAAGCGAAAAATGAAAATTAACGGAGGAAAAAAAAATTGGTATGATTTGGTAATTTTGTAGAAAtgacatcaa                                                                >  1:129857/1‑143 (MQ=255)
                                                                                       caGTAAATAATAAAAATTTCTCGAGGCGGGCACAGAAAAAAACAAAAGATCATTTTTTTTTTTCTTTACTGGTTTTAGAAGAGGAACAATAAAATCACAAGGGGGAAAAAAAAAGGGGATGATTGGGTATTTTGCATAAAATGACATAGaa                                                              >  1:589434/1‑151 (MQ=255)
                                                                                           aaaTAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGa                                                          <  2:1550047/151‑1 (MQ=255)
                                                                                            aataatAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGAt                                                         <  2:903161/151‑1 (MQ=255)
                                                                                                 taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                                                    <  2:1230660/151‑1 (MQ=255)
                                                                                                  aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa                                                   <  2:275522/151‑1 (MQ=255)
                                                                                                  aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATGTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa                                                   <  1:577171/151‑1 (MQ=255)
                                                                                                       tttCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGc                                              <  1:1024795/151‑1 (MQ=255)
                                                                                                            cGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAg                                         <  2:926201/151‑1 (MQ=255)
                                                                                                                    ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa                                     >  2:748310/1‑148 (MQ=255)
                                                                                                                       ccAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaa                               <  1:1243150/151‑1 (MQ=255)
                                                                                                                       ccAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaa                               <  1:739145/151‑1 (MQ=255)
                                                                                                                          gAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaaa                             >  1:812550/1‑150 (MQ=255)
                                                                                                                          gAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaaa                             >  2:1771933/1‑150 (MQ=255)
                                                                                                                           aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAcc                           <  1:827227/151‑1 (MQ=255)
                                                                                                                              aaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCttt                        >  1:190467/1‑151 (MQ=255)
                                                                                                                                 cGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt                     <  1:747948/151‑1 (MQ=255)
                                                                                                                                 cGACATGCACATTTTATTTTCACTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGAGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAGCCAACTGCATGAGTCTACTGAAAAGCCTTTGAt                     <  1:864937/151‑1 (MQ=255)
                                                                                                                                       gCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAtatgt                 >  2:1458680/1‑145 (MQ=255)
                                                                                                                                          caTTTTATTTTCCCTCATCTGTTTTGAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  1:275325/151‑1 (MQ=255)
                                                                                                                                          caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAGAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  2:1047912/151‑1 (MQ=255)
                                                                                                                                          caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  1:5729/151‑1 (MQ=255)
                                                                                                                                          caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  1:536347/151‑1 (MQ=255)
                                                                                                                                          caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            >  1:329772/1‑151 (MQ=255)
                                                                                                                                          caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            >  1:1043535/1‑151 (MQ=255)
                                                                                                                                          caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  1:297209/151‑1 (MQ=255)
                                                                                                                                          caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            >  2:1320604/1‑151 (MQ=255)
                                                                                                                                          caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  1:1749437/151‑1 (MQ=255)
                                                                                                                                          caTTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAatat            <  1:1073874/151‑1 (MQ=255)
                                                                                                                                            attttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAAtata           <  2:823436/151‑1 (MQ=255)
                                                                                                                                              ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  2:1130770/151‑1 (MQ=255)
                                                                                                                                              ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATAc          <  1:1720887/151‑1 (MQ=255)
                                                                                                                                               tttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACa         <  2:812968/151‑1 (MQ=255)
                                                                                                                                                ttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACAt        >  1:1019259/1‑151 (MQ=255)
                                                                                                                                                 tattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATg       <  2:783231/151‑1 (MQ=255)
                                                                                                                                                 tattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATg       >  1:1437380/1‑151 (MQ=255)
                                                                                                                                                      tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACCCAGAAAAGATTAACCAACTGCTTGAATCTACTGAAAAACCTTTGATGGATAATATACATGACACg  >  2:620442/1‑151 (MQ=255)
                                                                                                                                                      tCCCTCATCCGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAACTTTGATGGATAATATACATGACACg  >  1:1323426/1‑151 (MQ=255)
                                                                                                                                                       cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACg  >  1:156654/1‑150 (MQ=255)
                                                                                                                                                       cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACg  >  2:1054302/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                                   
TACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACG  >  NC_000913/2386736‑2387032

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