Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,812,7920CA40.3% 128.8 / 32.1 72T21N (ACC→AAC) dfpfused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (8/35);  new base A (29/0);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 3.51e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

CGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCC  >  NC_000913/3812642‑3812937
                                                                                                                                                       |                                                                                                                                                   
cGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAAc                                                                                                                                                     <  2:586354/151‑1 (MQ=255)
      gTGAATCTCAGGTGAAGGCACGACAACAAAAGGATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGA‑TTTCTGCCTTTAAAAaccc                                                                                                                                                  >  2:545790/1‑148 (MQ=255)
       tGAATCTCAGGTGAAGGCACGACAAGAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGa                                                                                                                                               <  2:47708/150‑1 (MQ=255)
             tcAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGA‑GTGCTGCCTTTAAAAAACCTGAATTg                                                                                                                                           >  1:292653/1‑148 (MQ=255)
              cAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGt                                                                                                                                         <  1:10854/149‑1 (MQ=255)
                 gTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTTGTGCCTCTAAAAacccccaaaatgtgcg                                                                                                                                      >  1:743019/1‑137 (MQ=255)
                 gTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGTTGCCTATAAAACCCCCGAACTGGTGcg                                                                                                                                      >  1:583604/1‑149 (MQ=255)
                    aaGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTCTAAAAACCCCGAACAGGTGggtcgct                                                                                                                                 >  2:788532/1‑149 (MQ=255)
                     aGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTTCTGCCATTAAACccccggccta                                                                                                                                            >  2:874772/1‑134 (MQ=255)
                        cACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTgcg                                                                                                                             <  2:1200426/151‑1 (MQ=255)
                            acaacaAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGAt                                                                                                                         <  2:809170/151‑1 (MQ=255)
                             caacaaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGCA‑TTGCTGCCTCTAAAAACCCCGAACTGGTGCGGCGCTTGcgcg                                                                                                                           >  1:1003881/1‑148 (MQ=255)
                              aacaaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTTCTGCCCATAAAAAACCTGAACTGGTggggcggt                                                                                                                                 >  1:876661/1‑134 (MQ=255)
                              aacaaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATcg                                                                                                                       <  2:884493/151‑1 (MQ=255)
                              aacaaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAAACCCCGAAATGGGGGGTCGTTTGCGCGATc                                                                                                                        >  2:426474/1‑150 (MQ=255)
                                  aaaGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGTTTCCTATAAAAACCCCCCACTGGTGCGTGGtttttgcggagccggg                                                                                                                   >  1:69508/1‑138 (MQ=255)
                                    aGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTTGCGCCTATAAAAACCCCGAACTTGTGCGTCCGTTGcgcgggcgc                                                                                                                      >  2:237135/1‑141 (MQ=255)
                                         tgtgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGAc                                                                                                             <  2:1139305/150‑1 (MQ=255)
                                          gtgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGt                                                                                                           <  2:618427/151‑1 (MQ=255)
                                           tgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGGTGCCTTTAAAAACCCTGAAATGGTGCGCCGTTTGCTGGATCGCCGGGCCGACg                                                                                                            >  2:664341/1‑149 (MQ=255)
                                            gATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTc                                                                                                          <  1:661962/150‑1 (MQ=255)
                                            gATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTc                                                                                                          <  1:814058/150‑1 (MQ=255)
                                               aaaaTCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGTTGCCTATAAAAACCCCGAACTGGTGCGGCGTTTGCGCGGTTGCGGGGCCGCCGGcccc                                                                                                       >  1:96790/1‑148 (MQ=255)
                                                aaaTCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTTGTGCCCCTAAAAACCCCGAACCGGTGGGTTGGTTTCGCGGTCGCgggg                                                                                                                  >  1:80634/1‑138 (MQ=255)
                                                aaaTCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTTCTGCCTATAAAAACCCCGAACTGGTGCGGCGGTTTCGCGGTCGCGGGGGCGACGTcc                                                                                                         >  2:571729/1‑147 (MQ=255)
                                                 aaTCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTa                                                                                                    <  2:711062/151‑1 (MQ=255)
                                                  aTCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAAACCCTGAACTGGGGCGGCGTTTGCGCGATCGCGGGGCCGACGGCcgcg                                                                                                      >  2:349960/1‑148 (MQ=255)
                                                    cGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGcc                                                                                                 <  1:571711/151‑1 (MQ=255)
                                                       cAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTTCTGCCCATAAAAAACCCTAACTGGGGGGGCGGTTGCGCGATCGCGGGGCCGGCGGCCGCGCAGCCATg                                                                                              >  1:436024/1‑151 (MQ=255)
                                                        aaCTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGa                                                                                             <  1:14134/151‑1 (MQ=255)
                                                           ttCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTTCGTATAAAAACCCCTAACTGCTGGGGGCTTTGTGCGGTTGCGGGGGCCGCGTCcgcgctgcgatctcc                                                                                           >  1:327898/1‑139 (MQ=255)
                                                            tcgggTGTGACCCGACAGGAAAAATCATCTTGAGCCTGTCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCAGTA‑TTGCTGTGTATAAAGCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGTCGAGGTCCGCGTAGCGATGAcc                                                                                           <  1:320327/146‑1 (MQ=255)
                                                             cTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAAACCCCGAACTGGGGCGTCGTTTGCGCGATCGCGGGGCCGACGGCCGCGTAGCCATGAcc                                                                                           >  2:586167/1‑148 (MQ=255)
                                                               ggTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAg                                                                                       <  2:114376/150‑1 (MQ=255)
                                                               ggTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGAA‑TTTGTGCCCAAAAAAACACCCAAAATGTGGGTCGCTTTTGCGGTCGCGGGGGCGGCGACCGCGGAGCCCTGGTCGaa                                                                                        >  1:486090/1‑149 (MQ=255)
                                                               ggTGTAACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGaa                                                                                        <  2:797771/149‑1 (MQ=255)
                                                                  gtCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCgg                                                                                    <  2:772804/150‑1 (MQ=255)
                                                                      cccTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCaaaa                                                                               <  2:1227518/151‑1 (MQ=255)
                                                                       ccTGCAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAg                                                                              <  2:678606/151‑1 (MQ=255)
                                                                             ggAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTa                                                                        <  2:1160907/151‑1 (MQ=255)
                                                                               aaaaaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATc                                                                      <  1:539076/151‑1 (MQ=255)
                                                                               aaaaaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGAGTTATGGCCCAAAAAAACCC‑CAACCAGTTCGTCGTTTATGCGTTCGCGGGGGCGGGGTGCCCGTAGCGCTGTCCGCAGCCGCAAAAGCCttt                                                                         >  1:1216322/1‑148 (MQ=255)
                                                                                  aaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAAACCCCGAACTGGGGGGGCGTTTGCGCGATCGCGGGGGCGACGTCCGCGGAGCCCTGAACGAAGCGGGAAAAGCCCTTATTAcc                                                                   >  2:902929/1‑151 (MQ=255)
                                                                                       tcatGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCCATAAAAACCCCGAACTGGGGCGGCGTTTGCGCGGTCGCGGGGGCGCCGGCCGCGCAGCCCTGACCGCAGCGGGAAAAGCCTTTTTCAcccccc                                                               >  1:372579/1‑148 (MQ=255)
                                                                                        catGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACtt                                                             <  2:897190/151‑1 (MQ=255)
                                                                                         atGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTa                                                            <  2:164968/151‑1 (MQ=255)
                                                                                           gAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAg                                                           <  2:64447/150‑1 (MQ=255)
                                                                                            aGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAAACCCCGAACTGGGGCGGCGGTTTGGCGGTCGCGGGGCCGGCGTCCCCGGAGCCATGAACGGAGCGGGAAAAACCCTTTTCAcccc                                                                 >  2:1131903/1‑143 (MQ=255)
                                                                                            aGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGA‑TTTCTGCCCCTAAAAACCCCTAACTGGGGCGGCGGTTTCGCGATCGCGGGGGCGACGGCCGCGGAGCCCTGACCGGAGCGGGAAAAAACTTTTTCTCCCCCCCTAGCt                                                         >  1:1078988/1‑151 (MQ=255)
                                                                                             gCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCtt                                                        <  2:39309/151‑1 (MQ=255)
                                                                                              ccTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTg                                                       <  2:583623/151‑1 (MQ=255)
                                                                                                  gCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑gg                                                  <  2:111/151‑1 (MQ=255)
                                                                                                     ggTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGGCTATAAAAACCCCTAACTGGGGCGTCGGTTGTGCGATCGCGGGGGCGGCGGCCGCCGAGCCATGACCCAAGCGGCAAAAGCCTTTTTCACCCCCCCTAGCTTGCA‑Ggg‑gg                                              >  1:271640/1‑151 (MQ=255)
                                                                                                       tAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCCATAAAAACCCTGAACTGGGGCGTCGGTTGCGCGATCGCGGGGGCGACGGCCGCGGAGCCATGACCGAAGCGGCAAAAACCCTTTTTACCCCCCCTAGCTTGCA‑GGC‑gggt                                            >  2:741723/1‑149 (MQ=255)
                                                                                                        aaaaaaaTCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGttt                                           <  1:270547/151‑1 (MQ=255)
                                                                                                        aaaaaaaTCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCCATAAAAACCCCGAACTGGTGCGTCGGTTGCGCGATCGCGGGGGCGACGTCCGCGGAGCCCTGACCGCAGCGGCAAAAGCCTTTTTTACCCCCCTTAACTTTCA‑GGC‑GGgtt                                           >  2:678174/1‑151 (MQ=255)
                                                                                                          aaaaaTCGTTCTCGGCGTTAGCGGCGGTA‑TTTCTGCCTATAAAACCCCTGAACTGGGGCCCCGTTTTTGCGGTCGCGGGGGCGGCAACCGCGGGGCCCTGACCGGAGCGGGAAAAACCTTTATTACCCCACCTTGCtttccgggggct                                              >  2:679873/1‑138 (MQ=255)
                                                                                                           aaaaTCGTTCTCGGCGTTAGCGGCGGTA‑TTTGTGCCCATAAAAACCCCCCACTGGGGGGGCGTTTTGGGGATCGCGGGGGCGGCGTCCGCGGCGGCATGACCGAAGCGGGAAAAGACTTTTTTACCCCCCCTAGCCTGCA‑GGC‑GGgtt                                           >  1:761976/1‑148 (MQ=255)
                                                                                                           aaaaTCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTg                                        <  2:366433/151‑1 (MQ=255)
                                                                                                             aaTCGTTCTCGGCGTTAGCGGCGGTA‑TTTCTGCCTATAAAACCCCTTAACTGGTGGGGCGGTTGCGCGATCGCGGGGCCGACGTCCGCGGAGCCATGACCGAAGCGGCAAAAGCCTTTTTCACCCCCCTTTGCTTGCA‑GGC‑GGgtttt                                         >  2:195483/1‑146 (MQ=255)
                                                                                                               tCGTTCTCGGCGTTAGCGGCGGTA‑TTTCTTCCTATAAAAAACCTTAACTGGGGCGGCGTTTGCGCGGTCGCGGGGCCGCCGTCCGCCGAGCCCTGACCGAAGGGGCAAAAACCTTTTTTACCCCACCTAGCTTTGA‑GGCGGGTTTTTGGt                                      >  1:1128871/1‑150 (MQ=255)
                                                                                                                  ttCTCGGCGTTAGCGGCGGTA‑TTTCTGCCTATAAAAACCCTGAAATGGTGGGGCGTTTGCGCGATCGCGGGGGCGACGTCCGCGGAGGCATGACCGCAGCGGGAAAAGGCTTTTTTACCCCCCCTAGCTTGCA‑GGG‑GGGTTTTgggt                                     >  2:992547/1‑145 (MQ=255)
                                                                                                                   tctcGGCGTTAGCGGCGGTA‑TTTCTTCCTATAAAAAACCCTAACTGGGGGGTCGTTTTCGCGGTCGCCGGGGCGGCGTCCGCGTAGCCCTGAACGAAGGGGGAAAAACCTTTATTACCCCACCTAGCTTGCA‑GGCGGGTTTTTGGtttt                                   >  1:810304/1‑147 (MQ=255)
                                                                                                                      cGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCg                               <  2:903270/149‑1 (MQ=255)
                                                                                                                      cGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGt                              <  1:783357/150‑1 (MQ=255)
                                                                                                                         cGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTcc                          <  2:464012/151‑1 (MQ=255)
                                                                                                                             aGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTCCGAc                       <  2:1075903/150‑1 (MQ=255)
                                                                                                                             aGCGGCGGTA‑TTGCTGCCTATAAAAACCCCGGACCGGGGCGGCGTGTTTGCCATCGCGGGGGCGACGGCCGCCGAGCCATGACCGGAGCGGCCAAAGCCCTTTTTACCCCACTTTGCCTCCAGGGG‑GGTTTTTGGTTTTCCCCttttcccc                       >  1:626497/1‑146 (MQ=255)
                                                                                                                                   ggTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTctgc                <  2:768367/151‑1 (MQ=255)
                                                                                                                                    gTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTctgc                >  1:165061/1‑150 (MQ=255)
                                                                                                                                    gTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGCTTCTGGTTTTCCCGTTTTCGGCAGTctgc                >  1:591173/1‑150 (MQ=255)
                                                                                                                                        ttGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTgg             >  2:472945/1‑150 (MQ=255)
                                                                                                                                            tgcCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGAccc         >  2:896046/1‑150 (MQ=255)
                                                                                                                                             gcCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGAccc         <  2:80635/149‑1 (MQ=255)
                                                                                                                                                  taAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTCCGCCCGTTTGCTGGACCCGGCAGcc  >  2:147442/1‑151 (MQ=255)
                                                                                                                                                   aaaaCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGcc  <  2:743053/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTA‑TTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCA‑GGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCC  >  NC_000913/3812642‑3812937

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: ‑