Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,813,2720CT44.2% 58.6 / 42.2 87S181L (TCG→TTG) dfpfused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (12/36);  new base T (38/0);  total (50/37)
Fisher's exact test for biased strand distribution p-value = 3.66e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGT  >  NC_000913/3813135‑3813422
                                                                                                                                          |                                                                                                                                                         
gcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGTGGGCCATTTTTCGccggctaccc                                                                                                                                                >  1:42952/1‑142 (MQ=255)
gcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTTCCCCTCAACGa                                                                                                                                               >  2:406838/1‑150 (MQ=255)
gcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTGCCCCTCCACGa                                                                                                                                               >  2:1468072/1‑150 (MQ=255)
gcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTTCCCgcccaaaa                                                                                                                                               >  2:1234809/1‑143 (MQ=255)
gcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAACGCCTTTTTTGCCCCGCAAcca                                                                                                                                               >  2:661530/1‑148 (MQ=255)
gcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCATTTTTTGCCCCCCAACGa                                                                                                                                               >  1:1446296/1‑150 (MQ=255)
 caTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGAcc                                                                                                                                             <  2:560302/151‑1 (MQ=255)
 caTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGAcc                                                                                                                                             <  2:1249672/151‑1 (MQ=255)
     aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCGTTTTTGCCCCGCAAcacactgga                                                                                                                                         >  1:483788/1‑143 (MQ=255)
     aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCCTTTTTTGCCCCTCAACGACCCGaa                                                                                                                                        >  2:1252613/1‑151 (MQ=255)
     aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCATTTTTTGCCCCTCAAAGAACTCaa                                                                                                                                        >  2:1114789/1‑151 (MQ=255)
       ttAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAAc                                                                                                                                      <  2:272357/151‑1 (MQ=255)
         agcggTGCTTGCTTCCCGTGGTGTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCGTTTTTCTCCCCTCAAAGACCTAAAAcca                                                                                                                                    >  1:1125480/4‑149 (MQ=255)
             gTGCTTGCTTCCCGTGGNNNGCTCATCTNNNNNNNAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTg                                                                                                                                 <  1:377010/150‑1 (MQ=255)
                cttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAAc                                                                                                                              <  2:1249793/150‑1 (MQ=255)
                cttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAAc                                                                                                                              <  2:235381/150‑1 (MQ=255)
                     ttCCCGTGGTTTGCTCATCTGGGGGCCGGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCGTTTTTTCCCCTCAAACACCCGTAACATCTTAACATTTTg                                                                                                                        >  1:1051724/1‑151 (MQ=255)
                     ttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGTTATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCATTTTTTTCCCCTCCAAGACCTGCAACATCTTGACATTATg                                                                                                                        >  1:892521/1‑151 (MQ=255)
                      tCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGa                                                                                                                       <  2:134763/151‑1 (MQ=255)
                      tCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGa                                                                                                                       <  1:718938/151‑1 (MQ=255)
                      tCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGa                                                                                                                       <  2:808421/151‑1 (MQ=255)
                        ccGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGAtt                                                                                                                    <  1:1020368/151‑1 (MQ=255)
                               ttGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTTCCCGTCAACGACCTGAAAAATCTGAACATTATCATTAccgcc                                                                                                              >  1:17010/1‑150 (MQ=255)
                               ttGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCCACGACCTGAAAAATCTGAACATTATGATTTcccccg                                                                                                             >  2:10937/1‑151 (MQ=255)
                                      tCTGGGGGCCATACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGAGCGTTTTTTGCGCGTCTACAAGCTGAAAAAATCTAAAAATTTGTTTACCGGCCGCCCGa                                                                                                       >  1:1356846/1‑150 (MQ=255)
                                      tCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTCGCGCACTTTTTGCCCGCCCACGCCCCTAAAAAACTTAAAAATTTATTTACCGCCGGCCCgg                                                                                                       >  1:563453/1‑149 (MQ=255)
                                      tCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTTCCCCCCAAAAACCCCAAACATCTGAAAATTATGATTTACccccgccccg                                                                                                        >  2:1449420/1‑150 (MQ=255)
                                             gCCAGACAGTGGCAGTCAGGCTTGTGGTGATGTCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCCTTTTTCGCCCCCCCAACACCTGAAAAATCTGAACATTATGATTAACGCCGGCCCCAcgcgctga                                                                                              >  1:947805/1‑147 (MQ=255)
                                              ccggACAGTGGCAGTCAGGCTTGTGGGGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTAAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAAc                                                                                             <  1:136324/148‑1 (MQ=255)
                                              ccAGACGGTGGCAGTCAGGCTTGTGGGGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAAc                                                                                             <  1:401665/150‑1 (MQ=255)
                                              ccAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAAc                                                                                             <  2:17024/151‑1 (MQ=255)
                                              ccAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATGGTGGATATGGCGGTAGGGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAAc                                                                                             <  2:257635/151‑1 (MQ=255)
                                              ccAGACAGTGGCAGGGAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGTTAGCGCATTTTTCCCCCGTCAACGACCTGAAACATCTGACCATTATGATTACCGCCGGCCCGACGCGTGAAc                                                                                             <  2:541266/151‑1 (MQ=255)
                                                   cAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGTCAAAGACCCGAAAAATCTGAACATTATTATTAACGCCGGCCCCACGCGGGaaaccct                                                                                         >  2:44552/1‑145 (MQ=255)
                                                     gTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGCCAACGACCTGAAACATCTGAAAATTTTGATTAACCCCGGCCCCCCCCGGGAACCCCCCGa                                                                                      >  2:379019/1‑151 (MQ=255)
                                                         cAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCACCGACCTGAAACAACTGAAAATTATTATTAACGCCCGCCCGGCGCGTggaacccccgattcg                                                                                  >  1:1409322/1‑137 (MQ=255)
                                                         cAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCCGCAACGACCTGAAAAATCTTAAAATTATGATTACCCCCGGCCCCACGCCGGGAACCCCCGATCCg                                                                                  >  2:532582/1‑151 (MQ=255)
                                                               ggCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTTGTGCCTTTTTCTCCCGTCAACCACCTGAAACAACTGAAAATTAAGATTACCGCCGGCCCCGCGCGGGAAACGCCCCAACCGGGGCgg                                                                            >  1:1319512/1‑150 (MQ=255)
                                                               ggCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGGCCCTCCACGACATGGAACATCTGAACATTATGATTACCGCCCGGCCGACGCGTGGAACGCCCGCTCCGGGGCGt                                                                            >  1:1295120/1‑151 (MQ=255)
                                                               ggCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTCCCCGCCAACCACCCGCAAAATCTTAACCTTATGATTAACGCCGGCCCGCCCCGGGAACCCCTCCATcccgggggg                                                                            >  1:1397856/1‑144 (MQ=255)
                                                               ggCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCCTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCCCCGCGTGAACCCCCCGATCCGGTGCGt                                                                            >  1:43086/1‑151 (MQ=255)
                                                                gCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCTTTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGc                                                                              <  2:1091995/148‑1 (MQ=255)
                                                                gCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAACGCATTTTTTGCCCCGCAACGAACTTGAAAAACTGGACATTATTATTACCGCCCGCCCCACGCCGGAACCCCCCGAttccggggggt                                                                           >  2:37835/1‑141 (MQ=255)
                                                                    gtggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGAAGCCCATTTTTTCCCCGTCCACCACCTGAAACAATTTAACATTATTATTTTCCCCGGCCCCCCCCGTTAACCCCTCGCCCCCGTGCGTTatat                                                                       >  2:389661/1‑151 (MQ=255)
                                                                     tggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTCGCGCGTTTTTTTACCGGCCAAAACCTTAAACATCTGGACATTATCATTCCCCCCGGCCCGGCGCGTGGACCGCCCGCTCCGGGGCGTtatct                                                                       >  2:72661/1‑148 (MQ=255)
                                                                     tggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATTATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATAtc                                                                      <  2:1295756/151‑1 (MQ=255)
                                                                     tggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTatat                                                                       <  1:1467303/150‑1 (MQ=255)
                                                                     tggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTAGGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTatat                                                                       <  2:112346/150‑1 (MQ=255)
                                                                     gggtgTTTTCGTTCCTGGGGAATTGCTCGATCCGTTACCCATTGTGGATATGGCGGTAGCGAATTTTTCGCCCGTCAAGGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTatat                                                                       <  1:1482396/149‑1 (MQ=255)
                                                                       gtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGGGCATTTTTTGCCCCTCCACCACCCGAAAAATCTGAACATTATGTTTAACGCCGGCCCCACGCGGGAACCCCCCCATCCGGTGCgcttaatct                                                                     >  1:201904/1‑142 (MQ=255)
                                                                        tgATATCGGTCCGGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATAtctc                                                                    <  2:1357524/150‑1 (MQ=255)
                                                                        ggATATCGGTCCTGGGGGAATGCCCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGGCAACGACCTGAAACATCTAAACATTATAATAACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATAtctc                                                                    <  2:918524/149‑1 (MQ=255)
                                                                         gATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGTCAACGACCTGAAAAATCTGAACATTATTATTACCCCCCGCCCCACGCGTGGAACGCCCGATACGGTGGGTTTTATCTCTa                                                                  >  2:212022/1‑151 (MQ=255)
                                                                            atCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAACCATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATc                                                               <  1:880284/151‑1 (MQ=255)
                                                                              ctggccagcggataTGCGCGATCGGTTAACCTTTGTGAATGTTGCGGTAGCGTAGTTTTGGGCCGGACAGGCCTCGAGCCAGCTAAGCTGTGTGATTACCGCCGGCCCGACGGGTGAACCGCTCGAGCCGGTGCGTTATATCTCTAATcac                                                             <  1:215831/138‑1 (MQ=255)
                                                                              cGGTCCTGGGCGAATGCTCGATCCTTTACCCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATcac                                                             <  2:590008/151‑1 (MQ=255)
                                                                              cGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATcac                                                             <  2:1202621/151‑1 (MQ=255)
                                                                               ggTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGTCAACGACCCGAAAAATCTGAACATTTTGATTACCCCCCGCCCCACGCGGGAACCGCTCCATCCGGGGCGGTTTATTTCTTATcaca                                                            >  2:204334/1‑151 (MQ=255)
                                                                                ttCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCATTTTTTCCCCCGCCACGACCTGAAACATCTGTACCTTTTGGTTTCCGCCGGCCCCCCCCGCGAACCCGCCGATCCCGTTCGTTATATCTCTCATcacaa                                                           >  1:1233282/2‑150 (MQ=255)
                                                                                gTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTCTGGCCTTTTTTGCCCCTCCACGACCCTAAAAATCTGGCCATTATGATTACCGCCGGCCCGCCCCCGGTACCCCTCCCTTCCGGGCGTTGTATCTCTAATccac                                                            >  1:487621/1‑147 (MQ=255)
                                                                                gTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGTGCGTTTTTCTCCCGTCCACCACCCGGAAAATCTGAAAAATATGTTTATCGCCGGCCCGACGCGTGAACCGCTCGCACCCGTGCGTTAATTCTCTTATcccaa                                                           >  1:24526/1‑150 (MQ=255)
                                                                                gTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAAACACCTGAAAAAACTTAAAATTATGATTTCTCCCTGCCCGCCGACTGAACACCCACATCCCGTGCGGTATATCTCTCAttccag                                                           >  1:1192558/1‑146 (MQ=255)
                                                                                 tCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGc                                                          <  2:1328643/151‑1 (MQ=255)
                                                                                      ggCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTTCCCGCCCACGACCCGAAAAAACTGAAAATTATGATTAACGCCGGCCCGACGCGTGAACCGCCCGATCCGGTGCGTTTTATCTCTAATcacaacctcgg                                                     >  1:650284/1‑144 (MQ=255)
                                                                                        cGAATGCTCGTTCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCGGCCAACGACCTGCACCACCTGACCGTTAGGATTACAGCCGGCCCGACGCGTGAACCGCTCGATCCGGTTTGTTATATCTCTAATCACAGCTCCGGCa                                                   <  1:1110/151‑1 (MQ=255)
                                                                                        cGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCa                                                   <  2:655870/151‑1 (MQ=255)
                                                                                         gAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTAGGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCa                                                   <  1:1334310/150‑1 (MQ=255)
                                                                                                  aTCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCCTCAACGACCTGAAAAATCTGAAAATTATGATTACCCCCGGCCCGACGCGTGAACCGCCCGATCCGGTGCGTTTTATCTCTAATCACAGCTCCCGCAAGATGGGgtt                                         >  2:481283/1‑151 (MQ=255)
                                                                                                    ccGTTAACCATTGTGGATATGGCGGTTGCGCCTTTTTTGCCCCCCCACGACCTGAAACATCTTGACATTATGATTAACGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTTTATCTCTAATCACCGCTCCCGCAAGATGGGTTTTg                                       >  1:242035/1‑151 (MQ=255)
                                                                                                    ccGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGGGAAACGCCCGATCCGGGGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGGTTTg                                       >  2:1174977/1‑151 (MQ=255)
                                                                                                    ccGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCCGCCCGACGCGTGAACCGCCCGATCCGGTGCGTTATATCTCTAATCACAGCCCCGGCAAAATGGGGTTTg                                       >  2:1207942/1‑151 (MQ=255)
                                                                                                    ccGTTAACCATTGTGGATATGGCGGTAACGCATTTTTTGCCCGCCAAAGACCTGAAAAATCTGAACATTATGATTACCGCCCGCCCGACACGGGAAACGCTCGATCCGGGGCGGTATATCCCTAATCACAACCCCGGCCAGAAGGGGTTTg                                       >  2:1167142/1‑151 (MQ=255)
                                                                                                    ccGTTAACCATTGTGGATATGGCGGGAGCGCCTTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCCCCGGCCCGACGCGGGGACCCCTCCATCCCGGGCGTTATATCTCTAATCACAGCCCCCGCCAGATGGGGTTTg                                       >  2:1131425/1‑151 (MQ=255)
                                                                                                     cGTTAACCATTGTGGTTATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGc                                      <  1:185221/151‑1 (MQ=255)
                                                                                                      gTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGTCCACAACCTGAAACATTTGAACATTTTGATTACCGCCGGCCCGACGCGTGGAACGCCCCATCCGGGGCGGTCAATTTCTTATCACAGCCCCGGGCAGATGGGgtttttt                                     >  2:38413/1‑148 (MQ=255)
                                                                                                      gTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTCCCCCCCAACGCCCTGCAACATCTTAACATTATGATTACCGCCCGCCCGACGCGTGGACCGCCCGAGCCGGTGCGTTTTTTCCCTAATCCCAACTCCCGGAAGATGGGGTTTg                                       >  1:1242754/1‑149 (MQ=255)
                                                                                                           ccATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGAGCCGGGGCGTTATATCTCTAATCACACCTCCGGCAAGAGGGGGTTTGCTTTcg                                 >  1:1154747/1‑150 (MQ=255)
                                                                                                           ccATTGTGGATATGGCGGTAGCCCCTTTTTTGCCCCGCAACCACCCTAAACATCTGAACATTATGATTACCCCCCGCCCCACGCGTGAACCCCCCCATCCGGTGCGTTATATCTCTAATCACAACTCCGGCAAGATGGGGTTTGCTATct                                 >  1:608737/1‑149 (MQ=255)
                                                                                                                  ggATATGGCGGTCGCGCATTTTTCGCCCGTGAACGACCTGCAACATCTGAACATTATGGTTTACGCCCGCCCGCCGCGCGAACCGCCCGATCCGGGGCGGTATATCTCTAATCACAGCTCCGGCAAGATGTGGTTTTCTATcgccgccgc                          >  1:17809/1‑150 (MQ=255)
                                                                                                                  ggATATGGCGGTAGCGCATTTTTTGCCCCTCAAAGAcc                                                                                                                                             >  1:443884/1‑38 (MQ=255)
                                                                                                                  ggATATGGCGGTAGCGCATTTTTCTCCCGTCAACGACCTGAAACATCTGAACATTATGATTATCTCCGCCCCCACGCGCGAAACGCTCCATTCCGGGCGTTATATCTTTAATTGCCGCCCTGGCGCGATTGGGTTTGGTCTcgcccgcgc                          >  1:865371/1‑145 (MQ=255)
                                                                                                                  ggATATGGCGGTAGCGCATTTTTCGCCCGTCAACGAcc                                                                                                                                             <  2:444039/38‑1 (MQ=255)
                                                                                                                         gCGGTAGCGCATTTTTCGCCCTTAAACGCCCTGAAACATCTAAACATTATGATTCCCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGcc                   <  1:399382/150‑1 (MQ=255)
                                                                                                                         gCGGTACCGCATTTTGCGCCCGTCAACCACCTGAAACATCTGAACAGTAGGATTACCGCCGGCCAGACGCGTGAACCGCTCGATCCGGGGCGTTATATCTCGAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGcc                   <  1:215862/150‑1 (MQ=255)
                                                                                                                                    cgttttCGCCCGTCAACGACCGGAACCATCTGACCATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGgcg         <  2:467500/148‑1 (MQ=255)
                                                                                                                                        ttCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAAc    <  2:1320165/151‑1 (MQ=255)
                                                                                                                                         tCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTAAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACg   <  2:1171265/151‑1 (MQ=255)
                                                                                                                                          cGCCCGTCAACGCCCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGt  <  2:1033245/151‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                         
GCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCGCGAACGT  >  NC_000913/3813135‑3813422

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