Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,138 (A)6→7 coding (15/594 nt) bdcR → transcriptional repressor for divergent bdcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,1321.A98.7% 253.0 / ‑4.1 78coding (9/594 nt)bdcRtranscriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (44/33);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTA  >  NC_000913/4474000‑4474265
                                                                                                                                     |                                                                                                                                      
gATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑aaaaaa                                                                                                                                  >  2:928807/1‑139 (MQ=255)
        gACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAAtctc                                                                                                                          <  2:681941/139‑1 (MQ=255)
         aCTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTcg                                                                                                                         >  2:926852/1‑139 (MQ=255)
                 cGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAg                                                                                                                 >  1:929339/1‑139 (MQ=255)
                            aCTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGa                                                                                                      <  1:928807/139‑1 (MQ=255)
                                tCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGtt                                                                                                  >  2:545704/1‑139 (MQ=255)
                                 cTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                 <  2:975679/139‑1 (MQ=255)
                                     ccTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCt                                                                                             >  1:1084954/1‑139 (MQ=255)
                                                gcCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGc                                                                                  >  1:594156/1‑139 (MQ=255)
                                                 cccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                 <  2:560652/123‑1 (MQ=255)
                                                 cccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                 >  1:560652/1‑123 (MQ=255)
                                                   cATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAt                                                                               <  1:926852/139‑1 (MQ=255)
                                                              cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaa                                                                    >  1:199741/1‑139 (MQ=255)
                                                              cTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaa                                                                    >  2:617539/1‑139 (MQ=255)
                                                                      attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCtt                                                            >  1:616533/1‑139 (MQ=255)
                                                                      attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCtt                                                            >  1:1087998/1‑139 (MQ=255)
                                                                      attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCtt                                                            >  2:347154/1‑139 (MQ=255)
                                                                      attatGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCtt                                                            >  2:1119223/1‑139 (MQ=255)
                                                                               cATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTcg                                                                                                                         >  1:1027847/1‑69 (MQ=255)
                                                                               cATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTcg                                                                                                                         <  2:1027847/69‑1 (MQ=255)
                                                                               cATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaa                                                   >  1:651366/1‑139 (MQ=255)
                                                                                 tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCGGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                 >  1:174174/1‑139 (MQ=255)
                                                                                 tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                 >  2:22327/1‑139 (MQ=255)
                                                                                 tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                 >  1:1156409/1‑139 (MQ=255)
                                                                                 tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                 >  1:500120/1‑139 (MQ=255)
                                                                                 tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                 >  1:460866/1‑139 (MQ=255)
                                                                                 tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                                 >  1:223300/1‑139 (MQ=255)
                                                                                  gCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                            >  1:811833/1‑95 (MQ=255)
                                                                                  gCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                            <  2:811833/95‑1 (MQ=255)
                                                                                   cATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAg                                               >  2:612901/1‑139 (MQ=255)
                                                                                    aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACg                                                                                                    >  2:161329/1‑85 (MQ=255)
                                                                                    aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACg                                                                                                    <  1:161329/85‑1 (MQ=255)
                                                                                      tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                  <  2:924718/133‑1 (MQ=255)
                                                                                      tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaa                                                  >  1:924718/1‑133 (MQ=255)
                                                                                      tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGtt                                            >  1:107998/1‑139 (MQ=255)
                                                                                      tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGtt                                            >  1:1228801/1‑139 (MQ=255)
                                                                                        cAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaa                                                                    <  2:1136054/113‑1 (MQ=255)
                                                                                        cAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCaa                                                                    >  1:1136054/1‑113 (MQ=255)
                                                                                        cAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAg                                                                 >  2:834679/1‑116 (MQ=255)
                                                                                        cAAAACTGTTTTAACTTTCTGTCAACACCTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGATCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAg                                                                 <  1:834679/116‑1 (MQ=255)
                                                                                        cAAAACTGTTTTAAC‑TTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCg                                         <  1:1302117/139‑1 (MQ=255)
                                                                                         aaaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                       >  2:703114/1‑93 (MQ=255)
                                                                                         aaaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                       <  1:703114/93‑1 (MQ=255)
                                                                                          aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                        <  1:617539/139‑1 (MQ=255)
                                                                                          aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                        <  2:199741/139‑1 (MQ=255)
                                                                                          aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                        <  1:1119223/139‑1 (MQ=255)
                                                                                           aaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCTAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAATGTTTCGAt                                       >  1:780313/1‑139 (MQ=255)
                                                                                             cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGAt                                       >  1:620683/1‑137 (MQ=255)
                                                                                             cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGAt                                       <  2:620683/137‑1 (MQ=255)
                                                                                              tGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCt                                    <  2:500120/139‑1 (MQ=255)
                                                                                                ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                            <  2:635098/81‑1 (MQ=255)
                                                                                                ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTc                                                                                            >  1:635098/1‑81 (MQ=255)
                                                                                                ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAtctc                                                                            <  1:172700/97‑1 (MQ=255)
                                                                                                ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAt                                                                               <  2:861596/94‑1 (MQ=255)
                                                                                                ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAt                                                                               >  1:861596/1‑94 (MQ=255)
                                                                                                ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                  >  2:230834/1‑139 (MQ=255)
                                                                                                ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                  >  2:815781/1‑139 (MQ=255)
                                                                                                ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                  >  1:771114/1‑139 (MQ=255)
                                                                                                ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AAAAAACCAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAAtctc                                                                            >  2:172700/1‑97 (MQ=255)
                                                                                                   tAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAg                               <  2:260877/139‑1 (MQ=255)
                                                                                                          cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                       >  1:885929/1‑76 (MQ=255)
                                                                                                          cTGTCAACAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                       <  2:885929/76‑1 (MQ=255)
                                                                                                               aaCAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGt                   <  1:230834/139‑1 (MQ=255)
                                                                                                               aaCAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGt                   <  1:612901/139‑1 (MQ=255)
                                                                                                                 cAGGTTTTGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTa                 <  2:107998/139‑1 (MQ=255)
                                                                                                                     ttttGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa             <  2:594156/139‑1 (MQ=255)
                                                                                                                     ttttGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGa             <  2:771114/139‑1 (MQ=255)
                                                                                                                       ttGTAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAtt           >  1:378429/1‑139 (MQ=255)
                                                                                                                          tAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATc        <  1:22327/139‑1 (MQ=255)
                                                                                                                          tAATGGTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATc        <  2:616533/139‑1 (MQ=255)
                                                                                                                              ggTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                 <  2:1162365/46‑1 (MQ=255)
                                                                                                                              ggTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                                 >  1:1162365/1‑46 (MQ=255)
                                                                                                                              ggTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGt                   >  2:1023202/1‑124 (MQ=255)
                                                                                                                              ggTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGt                   <  1:1023202/124‑1 (MQ=255)
                                                                                                                              ggTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTgg    >  2:1214609/1‑139 (MQ=255)
                                                                                                                              ggTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTgg    >  1:408221/1‑139 (MQ=255)
                                                                                                                              ggTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTgg    >  1:1161857/1‑139 (MQ=255)
                                                                                                                               gTCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGt   >  1:972190/1‑139 (MQ=255)
                                                                                                                                tCACTAAAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTa  <  2:972190/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
GATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AAAAAA‑CAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGGTA  >  NC_000913/4474000‑4474265

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑