Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,804,946 +T coding (132/1203 nt) proV → glycine betaine/proline ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,804,9461.T89.1% 453.6 / 46.2 129coding (132/1203 nt)proVglycine betaine/proline ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base . (8/6);  new base T (51/64);  total (59/70)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 132 total aligned reads displayed.

AGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAA  >  NC_000913/2804798‑2805095
                                                                                                                                                     |                                                                                                                                                     
   aaTCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTataa                                                                                                                                                    >  3:191423/1‑149 (MQ=255)
        ttCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGAc                                                                                                                                                 <  2:1045832/148‑1 (MQ=255)
        ttCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACg                                                                                                                                                <  3:528865/149‑1 (MQ=255)
        ttCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACg                                                                                                                                                <  3:862596/149‑1 (MQ=255)
          cTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGcc                                                                                                                                              >  2:783389/1‑149 (MQ=255)
          cTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGt                                                                                                                                           >  2:1322542/1‑151 (MQ=255)
            aTTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGt                                                                                                                                           <  3:1615054/150‑1 (MQ=255)
               gCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTgg                                                                                                                                       <  2:457143/151‑1 (MQ=255)
               gCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTgg                                                                                                                                       >  2:444279/1‑151 (MQ=255)
                cATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGc                                                                                                                                      <  2:263236/151‑1 (MQ=255)
                  tGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGcc                                                                                                                                     >  3:1060931/1‑150 (MQ=255)
                  tGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGcc                                                                                                                                     >  2:1048599/1‑150 (MQ=255)
                   ggCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCAt                                                                                                                                   <  2:1585098/151‑1 (MQ=255)
                   ggCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCAt                                                                                                                                   <  2:15642/151‑1 (MQ=255)
                     cAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTg                                                                                                                                 >  3:976591/1‑151 (MQ=255)
                     cAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTg                                                                                                                                 >  2:771769/1‑151 (MQ=255)
                       attaaattaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTga                                                                                                                                >  2:883587/1‑150 (MQ=255)
                         taaattaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTga                                                                                                                                >  2:39267/1‑148 (MQ=255)
                          aaattaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTgaaga                                                                                                                             <  2:1565888/150‑1 (MQ=255)
                          aaattaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTgaaga                                                                                                                             >  2:131437/1‑150 (MQ=255)
                           aattaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGaagaag                                                                                                                           <  3:110705/151‑1 (MQ=255)
                            attaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAgg                                                                                                                          >  3:968330/1‑151 (MQ=255)
                             ttaGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCg                                                                                                                        >  2:477990/1‑151 (MQ=255)
                              taGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGc                                                                                                                         <  2:746610/150‑1 (MQ=255)
                               aGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCg                                                                                                                        <  3:1383010/150‑1 (MQ=255)
                                gAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCgag                                                                                                                      <  2:1003372/151‑1 (MQ=255)
                                   aTTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGata                                                                                                                   <  3:45412/151‑1 (MQ=255)
                                    ttAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGatat                                                                                                                  >  3:240069/1‑151 (MQ=255)
                                      aaaaaTCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATAttt                                                                                                                <  3:919419/151‑1 (MQ=255)
                                       aaaaTCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTg                                                                                                               >  3:250309/1‑151 (MQ=255)
                                          aTCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGtc                                                                                                             >  2:643885/1‑150 (MQ=255)
                                          aTCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGtc                                                                                                             >  3:1222517/1‑150 (MQ=255)
                                              ttATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGtca                                                                                                            <  2:731386/147‑1 (MQ=255)
                                              ttATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGtca                                                                                                            <  2:205231/147‑1 (MQ=255)
                                               tataAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATg                                                                                                       >  2:244924/1‑151 (MQ=255)
                                                 taAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATggg                                                                                                     >  3:412110/1‑151 (MQ=255)
                                                  aaaaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGATTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGATATATTTGTCATCATGGGa                                                                                                    <  3:897941/151‑1 (MQ=255)
                                                   aaaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAt                                                                                                   >  2:1234018/1‑151 (MQ=255)
                                                   aaaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAt                                                                                                   >  2:52401/1‑151 (MQ=255)
                                                   aaaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAt                                                                                                   <  2:903782/151‑1 (MQ=255)
                                                   aaaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAt                                                                                                   <  2:340476/151‑1 (MQ=255)
                                                    aaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAtt                                                                                                  <  2:1132225/151‑1 (MQ=255)
                                                    aaTATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAtt                                                                                                  <  3:572827/151‑1 (MQ=255)
                                                     atatTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGAtt                                                                                                  >  3:287661/1‑150 (MQ=255)
                                                      tatTTGGCGAGCATCCACAGCGAGCGTTAAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCTATTGAAGAAGGCGAGATATTTGTCATCATGGGATTa                                                                                                 >  3:523078/1‑150 (MQ=255)
                                                           ggCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCt                                                                                          <  3:1190598/151‑1 (MQ=255)
                                                            gCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCt                                                                                          >  2:694191/1‑151 (MQ=255)
                                                             cGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCg                                                                                        >  3:245514/1‑151 (MQ=255)
                                                                 cATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGt                                                                                     <  3:1056888/151‑1 (MQ=255)
                                                                  aTCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTa                                                                                    <  3:505037/151‑1 (MQ=255)
                                                                     cacaGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAAt                                                                                 >  3:856342/1‑151 (MQ=255)
                                                                     cacaGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAAt                                                                                 >  3:1106450/1‑151 (MQ=255)
                                                                      acaGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATc                                                                                >  3:767847/1‑151 (MQ=255)
                                                                      acaGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATc                                                                                >  3:212552/1‑151 (MQ=255)
                                                                        agcgagcgTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATcc                                                                               >  2:1432469/1‑150 (MQ=255)
                                                                        agcgagcgTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATcc                                                                               <  3:778319/150‑1 (MQ=255)
                                                                        agcgagcgTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCca                                                                              <  2:570241/151‑1 (MQ=255)
                                                                         gcgagcgTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCcac                                                                             <  3:777267/151‑1 (MQ=255)
                                                                               gTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGt                                                                       >  2:696132/1‑151 (MQ=255)
                                                                                     atatatCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCt                                                                 <  2:951637/151‑1 (MQ=255)
                                                                                      tatatCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTtc                                                               >  2:645206/1‑151 (MQ=255)
                                                                                       atatCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCtt                                                                <  3:577391/150‑1 (MQ=255)
                                                                                       atatCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCtt                                                                <  2:670329/150‑1 (MQ=255)
                                                                                       atatCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCtt                                                                <  3:1632710/150‑1 (MQ=255)
                                                                                        tatCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTtct                                                              >  3:783247/1‑151 (MQ=255)
                                                                                         atCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTtctc                                                             <  3:698184/151‑1 (MQ=255)
                                                                                         atCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTtct                                                              <  3:675061/150‑1 (MQ=255)
                                                                                           cGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCaa                                                           <  3:1240711/151‑1 (MQ=255)
                                                                                            gatcAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCAATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATATGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCa                                                            <  3:1206655/146‑1 (MQ=255)
                                                                                            gAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCGTGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAAt                                                          <  2:1536629/151‑1 (MQ=255)
                                                                                                  ggACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTg                                                    <  3:1371602/151‑1 (MQ=255)
                                                                                                     cTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTg                                                >  3:1097285/1‑151 (MQ=255)
                                                                                                      tttCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGTCTGATTg                                                >  3:339870/1‑151 (MQ=255)
                                                                                                       ttCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGAtt                                                 <  3:1508517/149‑1 (MQ=255)
                                                                                                       ttCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTg                                                >  2:368854/1‑150 (MQ=255)
                                                                                                          aaaaGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAacc                                            >  2:112490/1‑151 (MQ=255)
                                                                                                          aaaaGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAacc                                            >  2:334640/1‑151 (MQ=255)
                                                                                                          aaaaGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAaccc                                           >  2:510930/1‑151 (MQ=255)
                                                                                                              gAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAacccacc                                        >  3:210353/1‑151 (MQ=255)
                                                                                                              gAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAacccacc                                        <  3:558557/151‑1 (MQ=255)
                                                                                                              gAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAacccacc                                        >  3:482434/1‑151 (MQ=255)
                                                                                                                aCAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCcg                                      <  2:1334377/151‑1 (MQ=255)
                                                                                                                aCAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCcg                                      <  3:1222337/151‑1 (MQ=255)
                                                                                                                 cAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCcgc                                     >  2:962962/1‑151 (MQ=255)
                                                                                                                 cAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCcg                                      >  3:1209656/1‑150 (MQ=255)
                                                                                                                      tCTGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGc                                 >  2:601288/1‑150 (MQ=255)
                                                                                                                        tGGAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCa                                >  3:325484/1‑149 (MQ=255)
                                                                                                                         ggAAAAAACTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGt                             <  2:1042742/151‑1 (MQ=255)
                                                                                                                               aaCTGGGCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCtgat                       <  2:1171277/151‑1 (MQ=255)
                                                                                                                                aCTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCtgattg                     <  2:706212/151‑1 (MQ=255)
                                                                                                                                   gggCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCtgattga                    <  3:342245/149‑1 (MQ=255)
                                                                                                                                    ggCTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATg                  <  3:1374992/151‑1 (MQ=255)
                                                                                                                                      cTATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGgt                >  3:1080102/1‑151 (MQ=255)
                                                                                                                                       tATCGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCTGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGgtg               >  2:203738/1‑151 (MQ=255)
                                                                                                                                       tATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGgtgt              >  2:392886/1‑151 (MQ=255)
                                                                                                                                          cGCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTgg            <  2:1437024/151‑1 (MQ=255)
                                                                                                                                           gCTTGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGa           <  2:514319/151‑1 (MQ=255)
                                                                                                                                              tGGCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGat          <  3:1048521/149‑1 (MQ=255)
                                                                                                                                               ggCGTATAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATAtt       <  2:1142638/151‑1 (MQ=255)
                                                                                                                                                    aaaaGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCaa  >  2:121331/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGAATCTTTCTATTGCATGGCAATTAAATTAGAAATTAAAAATCTTTATAAAATATTTGGCGAGCATCCACAGCGAGCGTTCAAATATATCGAACAAGGACTTTCAAAAGAACAAATTCTGGAAAAAACTGGGCTATCGCTTGGCGTAAAAGACGCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAA  >  NC_000913/2804798‑2805095

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