Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,813,2720CT42.6% 69.7 / 28.2 69S181L (TCG→TTG) dfpfused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (9/30);  new base T (29/0);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCG  >  NC_000913/3813125‑3813415
                                                                                                                                                   |                                                                                                                                               
cTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGcc                                                                                                                                              >  1:982020/1‑151 (MQ=255)
cTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGcc                                                                                                                                              >  2:1095336/1‑151 (MQ=255)
  gCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGccc                                                                                                                                             <  2:240787/150‑1 (MQ=255)
     aCGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAACGCATTTTTTGCCCGTc                                                                                                                                          >  2:698923/1‑150 (MQ=255)
         agcaTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGa                                                                                                                                     <  1:574287/151‑1 (MQ=255)
            aTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTCGCCCGTCAACGAcc                                                                                                                                   >  2:661167/1‑150 (MQ=255)
            aTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGa                                                                                                                                     >  2:631825/1‑148 (MQ=255)
            aTAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAACGAcc                                                                                                                                   >  1:324231/1‑150 (MQ=255)
             tAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAACGAcc                                                                                                                                   >  2:968464/1‑149 (MQ=255)
              aaTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTg                                                                                                                                 <  1:124265/150‑1 (MQ=255)
               aTTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGCAGCGCATTTTTTGCCCGTCCACGAcc                                                                                                                                   >  2:1065162/1‑147 (MQ=255)
                   agagGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAAc                                                                                                                             <  1:737413/149‑1 (MQ=255)
                    gagGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGTGGGCCTTTTTTGCCCCTCCACAACCTGaaaaa                                                                                                                            >  1:683532/1‑147 (MQ=255)
                     agGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTTCCCGTCAACGACCCGaaa                                                                                                                              >  2:77043/1‑146 (MQ=255)
                      ggTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTTCCCCTCAACGACCCGaaa                                                                                                                              >  1:388426/1‑145 (MQ=255)
                        tgcttgctTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAAGGGAATTTTGCCCCCCGACGACCCGAAAAATccgta                                                                                                                      >  1:1122194/1‑146 (MQ=255)
                         gcttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCGTTTTCCCCCGCCCACCACCCGAAACATTTGaa                                                                                                                      >  1:772998/1‑150 (MQ=255)
                          cttgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCCCATTTTTTGCCCGTCAACGACCTGCAACATCTGaa                                                                                                                      >  2:70132/1‑149 (MQ=255)
                            tgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACa                                                                                                                    <  2:951998/149‑1 (MQ=255)
                            tgcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACAt                                                                                                                   <  1:1083484/150‑1 (MQ=255)
                             gcttCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACAtt                                                                                                                  <  1:110778/150‑1 (MQ=255)
                                tCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGTCCACCACCCGAAACATCTTaa                                                                                                                      >  2:379646/1‑143 (MQ=255)
                                 cccGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGAt                                                                                                             <  2:794342/151‑1 (MQ=255)
                                    gTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACCACCTGAAACATCCGAACATTATGAtt                                                                                                            >  1:277311/1‑149 (MQ=255)
                                      ggTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTAccg                                                                                                        <  2:475023/151‑1 (MQ=255)
                                       gTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTAccgc                                                                                                       <  1:517721/151‑1 (MQ=255)
                                       gTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTAccgc                                                                                                       <  2:733252/151‑1 (MQ=255)
                                          tGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAAAATTTGAAAATTATGATTAccgcc                                                                                                      >  2:136024/1‑149 (MQ=255)
                                           gCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCCTCCAAGACCTGAAAAATCTGAACATTATGATTaacccc                                                                                                      >  2:827748/1‑143 (MQ=255)
                                           gCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCgg                                                                                                    <  2:673585/150‑1 (MQ=255)
                                             tCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTTCCCCTCCACGACCTGAAACATCTGAACATTATTATTACCGCCGGccc                                                                                                 >  1:1018229/1‑151 (MQ=255)
                                               aTCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTTGTGCGTTTTTTTCCCGCCCACGGCCCCAAAAATTTGTACATTATGATTACAgcccgccc                                                                                                 >  2:619710/1‑149 (MQ=255)
                                               aTCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTTGCGCAGTTTTTTCCCCTCACCCACCTCCAAAAACTTAAAAATTTGATTTACGCCGGccc                                                                                                 >  1:223088/1‑149 (MQ=255)
                                               aTCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCCTTTTTTGCCCGGCAACGACCTGAAACATCTGAACATTATGATTACCgcccgccc                                                                                                 >  2:642648/1‑149 (MQ=255)
                                               aTCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAACGCATTTTTTGCCCGTCAAAGGACTGTAACATCTTAAAATTTAAATTACCgcccgccc                                                                                                 >  2:59237/1‑149 (MQ=255)
                                                tCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGGGCATTTTTTGCCCGCCAACCACCTGAAAAAACTGAAAATTATGATTTACCCCGGccc                                                                                                 >  1:1121225/1‑148 (MQ=255)
                                                   gggggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTACCCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGAcg                                                                                             <  2:502122/149‑1 (MQ=255)
                                                   gggggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAAAGACCTGAAACATCTGAACATTATTATTACCGCCCGCCCGAcgcg                                                                                           >  2:470168/1‑151 (MQ=255)
                                                   gggggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGAcgcg                                                                                           <  2:926786/151‑1 (MQ=255)
                                                   gggggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGAcgcg                                                                                           <  2:632375/151‑1 (MQ=255)
                                                     gggCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGAGGATTTTTTTCCCGTTAACCAGCCGGAACATCTGTACAATAAGTATACCGACCGccccgc                                                                                              >  1:714756/1‑143 (MQ=255)
                                                             cAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTGGGGCATTTTTCTCCCCCCCACCAACCTAAACATCTCAACAATATGATTAACCCCCGccccgcccctgcacccc                                                                                   >  1:827349/1‑135 (MQ=255)
                                                             cAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTc                                                                                 <  1:320546/151‑1 (MQ=255)
                                                              aGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAAAGACCCGAAACATCTGAAAATTATGATTACCGCCGGCCCGACGCGTGAACCGccc                                                                                 >  2:683800/1‑148 (MQ=255)
                                                                tGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAAAATTATGATTTCCGCCCGCCCGGCGCGTGAACCGccc                                                                                 >  1:224275/1‑146 (MQ=255)
                                                                       cAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCCTCAACGACCTGAAAAATCTGAACATTATGATTACCCCCGGCCCGACGCGTGAACCGCCCGATCCgggg                                                                        >  1:528224/1‑148 (MQ=255)
                                                                         ggCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATGGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGCCCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCg                                                                      <  1:339526/150‑1 (MQ=255)
                                                                          gCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGGAGCGCCTTTTTTTCCCGTCCACGACCTGGAACATCTTAACATTTTGATTACCGCCCGCCCGACGCGTGAACCCCTCGCTCCGGGGCg                                                                      >  1:687300/1‑149 (MQ=255)
                                                                               tggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAAACGCCCGATCCGGTGCGGTatat                                                                >  2:819051/1‑150 (MQ=255)
                                                                               tggtgATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTatat                                                                <  1:357654/150‑1 (MQ=255)
                                                                                ggtgATATCGGTCCTGGGCAAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATAtct                                                              <  2:470327/151‑1 (MQ=255)
                                                                                   gATATCGGTCCGGGGCGAATGCTGGATCCGTTAACCATTGTGGATATGGCGGTGGCGCATTTTGCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATctct                                                            <  2:288302/150‑1 (MQ=255)
                                                                                        cGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCCTCAACGACCTTCAACATCTGAACATTATGATTATCGCCGGCCCGCCGCGGGAACCGCTCGGACCGTTGCGTTATATCTCTAATca                                                       >  2:762517/1‑150 (MQ=255)
                                                                                                  cGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCa                                            <  2:212933/151‑1 (MQ=255)
                                                                                                       gCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTTCCCCTCAACGACCTGAAACATCTGAACATTTTGATTACCGCCGGCCCGACGCGGGAACCCCCCCATCCCGGGCGGTATTTCTCTAATCACAGCTCCCGCCAGATg                                       >  1:346046/1‑151 (MQ=255)
                                                                                                            aTCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTTGCCCGGCAACGACCTGAAAAATCTGAACATTATGATTACCGCCGGCCCCACCCGTGAAACGCTCGATCCGGGGGGTTATATCTCTAATCACAACTCCGGCCAGATGGGgtt                                  >  2:99469/1‑151 (MQ=255)
                                                                                                                 ttAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAAAGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCCCGATCCGGTGCGTTATATCTCTAATCACAACTCCGGCAAGATGGGTTTTGCt                              >  2:423113/1‑150 (MQ=255)
                                                                                                                    aCCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcg                          <  2:514523/151‑1 (MQ=255)
                                                                                                                         tgtgGATATGGCGGTAGCGCATTTTTTGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCCCGATCCCGGGCGTTATATCTCTAATCACAGCCCCGGCAAGATGGGGTTTTCTATcgccgcc                     >  1:676160/1‑151 (MQ=255)
                                                                                                                         tgtgGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgccgc                      <  2:533123/150‑1 (MQ=255)
                                                                                                                            ggATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACTTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATcgccgccgc                   <  2:658149/150‑1 (MQ=255)
                                                                                                                                  ggCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGc             <  1:185500/150‑1 (MQ=255)
                                                                                                                                     ggTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCcgt         <  2:800861/151‑1 (MQ=255)
                                                                                                                                       tAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCcgtc        <  1:432859/150‑1 (MQ=255)
                                                                                                                                        aGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCcgtcgt      <  2:127373/151‑1 (MQ=255)
                                                                                                                                         gcgcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTg     >  1:457210/1‑151 (MQ=255)
                                                                                                                                         gcgcATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTg     >  1:175935/1‑151 (MQ=255)
                                                                                                                                             aTTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGgcg  <  1:1078874/150‑1 (MQ=255)
                                                                                                                                                  tCGCCCGTCACCGTCCTGAATCTCTTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCTCTGCCCGTCGTg     <  1:59439/142‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CTGCCACGCAGCATAATTTAGAGGTGCTTGCTTCCCGTGGTTTGCTCATCTGGGGGCCAGACAGTGGCAGTCAGGCTTGTGGTGATATCGGTCCTGGGCGAATGCTCGATCCGTTAACCATTGTGGATATGGCGGTAGCGCATTTTTCGCCCGTCAACGACCTGAAACATCTGAACATTATGATTACCGCCGGCCCGACGCGTGAACCGCTCGATCCGGTGCGTTATATCTCTAATCACAGCTCCGGCAAGATGGGTTTTGCTATCGCCGCCGCCGCTGCCCGTCGTGGCG  >  NC_000913/3813125‑3813415

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: