Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,771,144 A→C 5.3% R161R (CGA→CGC) ‡ PP_5599 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,771,1440AC 5.3% 914.9 / 27.1 265R161R (CGA→CGC) ‡PP_5599LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (115/136);  new base C (3/11);  total (118/147)
Fisher's exact test for biased strand distribution p-value = 9.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 281 total aligned reads displayed.

TCAAGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGAATGACAATGACCTGCAG  >  NC_002947/3770995‑3771291
                                                                                                                                                     |                                                                                                                                                   
tCAAGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCga                                                                                                                                                     <  1:1008727/150‑1 (MQ=255)
tCAAGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCga                                                                                                                                                     <  1:6473131/150‑1 (MQ=255)
tCAAGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCga                                                                                                                                                     >  1:6716823/1‑150 (MQ=255)
   aGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGagaa                                                                                                                                                  <  1:5099924/150‑1 (MQ=255)
   aGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGagaa                                                                                                                                                  <  2:4412227/150‑1 (MQ=255)
   aGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGagaa                                                                                                                                                  <  1:2947305/150‑1 (MQ=255)
    gCATTGCCCCGCAACTGAGGCTACGCATCACGCCCGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCACAAGCATTCGTATACAGGTTCGTCgagacg                                                                                                                                                 >  2:5695591/1‑148 (MQ=255)
       ttGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAAc                                                                                                                                              <  1:3949801/150‑1 (MQ=255)
         gCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCTCGAAGAAcct                                                                                                                                            <  1:4187284/150‑2 (MQ=255)
         gCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCTCGAAGAAcct                                                                                                                                            <  1:3362840/150‑2 (MQ=255)
         gCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCTCGAAGAAcct                                                                                                                                            <  1:138375/150‑2 (MQ=255)
         gCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAAccc                                                                                                                                            <  2:2009432/150‑1 (MQ=255)
         gCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTGTCCGAACCAGTAGCCGGCGGAGAGCTCGATGCGGGAGTCATGGTAGGTCCGCCCTTCAGCGTCCCCACAAGCATTCGTATACAGGTTCTTCGAGAAGAAccc                                                                                                                                            <  2:4358793/150‑1 (MQ=255)
          ccccGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAAcccc                                                                                                                                           >  2:1561780/1‑150 (MQ=255)
           cccGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCt                                                                                                                                          <  2:1222602/150‑1 (MQ=255)
           ccagCAACTGAGGCTACGCATCACGACAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCAGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTACAGCGTCCCCAAAAGCAATCGTATACAGGTTCTTCGAGAAGAACCCCt                                                                                                                                          >  1:4611325/4‑150 (MQ=255)
            ccGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTg                                                                                                                                         >  2:3695306/1‑150 (MQ=255)
            ccGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTg                                                                                                                                         >  2:142828/1‑150 (MQ=255)
              gCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCTCGAAGAAcctgtata                                                                                                                                       <  2:3894956/150‑7 (MQ=255)
              gCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGAt                                                                                                                                       <  2:4800266/150‑1 (MQ=255)
                    gAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGTTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCtttt                                                                                                                                 <  1:247811/150‑1 (MQ=255)
                     aGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTg                                                                                                                                >  2:6627251/1‑150 (MQ=255)
                      ggCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGa                                                                                                                               >  1:899768/1‑150 (MQ=255)
                         tACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCa                                                                                                                            >  1:4299630/1‑150 (MQ=255)
                           cGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCaaa                                                                                                                          <  2:1055236/150‑1 (MQ=255)
                               tCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGAGGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGAcc                                                                                                                      >  1:1062117/1‑150 (MQ=255)
                                    ccAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGtt                                                                                                                 <  2:2009233/150‑1 (MQ=255)
                                         ggACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACaa                                                                                                            <  2:4974456/150‑1 (MQ=255)
                                           aCTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAgc                                                                                                          <  2:3075789/150‑1 (MQ=255)
                                            cTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAgcg                                                                                                         <  2:6745366/150‑1 (MQ=255)
                                             tCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAgcgc                                                                                                        <  2:5099095/150‑1 (MQ=255)
                                             tCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAgcgc                                                                                                        <  1:4564409/150‑1 (MQ=255)
                                             tCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAgcgc                                                                                                        >  2:185912/1‑150 (MQ=255)
                                               gcgcAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTTAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGa                                                                                                      >  1:5927720/1‑150 (MQ=255)
                                               gcgcAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGa                                                                                                      <  1:6250635/150‑1 (MQ=255)
                                                  cAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACAt                                                                                                   <  1:1597241/150‑1 (MQ=255)
                                                     tCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCg                                                                                                >  2:5784671/1‑150 (MQ=255)
                                                          aCGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTcga                                                                                           <  1:5250748/150‑1 (MQ=255)
                                                           cGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTcgac                                                                                          <  1:4065157/150‑1 (MQ=255)
                                                            gAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTcgacg                                                                                         >  2:5055798/1‑150 (MQ=255)
                                                            gAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTcgacg                                                                                         >  2:1326239/1‑150 (MQ=255)
                                                               cAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGAtt                                                                                      <  2:528304/150‑1 (MQ=255)
                                                                     gCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGt                                                                                >  2:4163375/1‑150 (MQ=255)
                                                                     gCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGt                                                                                <  1:5599003/150‑1 (MQ=255)
                                                                      ccTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTc                                                                               <  2:4543619/150‑1 (MQ=255)
                                                                        tGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGa                                                                             >  1:1181645/1‑150 (MQ=255)
                                                                          cGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGTGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGaaa                                                                           <  2:248662/150‑1 (MQ=255)
                                                                           ggAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAAc                                                                          >  2:2105769/1‑150 (MQ=255)
                                                                           ggAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAAc                                                                          <  2:3394740/150‑1 (MQ=255)
                                                                               ggCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGa                                                                      <  2:4257925/149‑1 (MQ=255)
                                                                               agCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGATGATTCTTGGTCGAAACCTGa                                                                      <  1:1750697/150‑1 (MQ=255)
                                                                               agCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGa                                                                      <  1:3782910/150‑1 (MQ=255)
                                                                               agCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGa                                                                      <  2:4104122/150‑1 (MQ=255)
                                                                                 cTCGATGCGTGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGAtt                                                                    <  1:3569548/150‑1 (MQ=255)
                                                                                 cTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGAtt                                                                    >  2:3575060/1‑150 (MQ=255)
                                                                                    gATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAATCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGa                                                                 <  1:4328672/150‑1 (MQ=255)
                                                                                    gATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGa                                                                 >  2:5455044/1‑150 (MQ=255)
                                                                                      tGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCCACGACCCCGTTGACAAGCGCGACATCCGCGCGACGATTCTTGGTCGAAACCTGATTCGata                                                               <  2:6389821/150‑1 (MQ=255)
                                                                                      tGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGata                                                               <  2:1746867/150‑1 (MQ=255)
                                                                                      tGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGata                                                               >  1:3124947/1‑150 (MQ=255)
                                                                                         gggAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGa                                                            <  2:6372536/150‑1 (MQ=255)
                                                                                          ggAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGAGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGAc                                                           >  1:714528/1‑150 (MQ=255)
                                                                                              tCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTcac                                                       >  1:2829920/1‑150 (MQ=255)
                                                                                                  ggTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGa                                                                 <  1:5083502/136‑1 (MQ=255)
                                                                                                      tGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTgg                                               <  1:4429797/150‑1 (MQ=255)
                                                                                                       gTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTggg                                              >  1:2134334/1‑150 (MQ=255)
                                                                                                        tCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTgggg                                             >  1:1879621/1‑150 (MQ=255)
                                                                                                        tCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTgggg                                             <  1:4504344/150‑1 (MQ=255)
                                                                                                        tCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTgggg                                             >  1:5424015/1‑150 (MQ=255)
                                                                                                          cGCCCTTCAGCGTCCCCAAAAGCATTAGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCGGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGTAACCTGATTCGATATGACTCACATGCCTGGGGTg                                           <  1:5314136/150‑1 (MQ=255)
                                                                                                           gCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTgg                                          <  2:2802571/150‑1 (MQ=255)
                                                                                                           gCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTgg                                          >  2:4853878/1‑150 (MQ=255)
                                                                                                           gCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTgg                                          >  2:5204562/1‑150 (MQ=255)
                                                                                                           gCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTgg                                          >  2:2680968/1‑150 (MQ=255)
                                                                                                           gCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTgg                                          >  1:2026516/1‑150 (MQ=255)
                                                                                                            cccTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTggg                                         <  1:1439829/150‑1 (MQ=255)
                                                                                                               ttCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAg                                      >  1:3704187/1‑150 (MQ=255)
                                                                                                                     gTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGcc                                >  1:2450232/1‑150 (MQ=255)
                                                                                                                     gTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGcc                                <  1:2037262/150‑1 (MQ=255)
                                                                                                                      tCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGccc                               >  2:6167977/1‑150 (MQ=255)
                                                                                                                        cccAACAGCATTCGTAGACAGGTTCGTCGAGAAGAACCCCTGATGCTTGTGACCAAAGACCCCGTTGACAAGAGCGACATCCGCTCGCCGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCcag                             <  2:5629401/150‑1 (MQ=255)
                                                                                                                         ccAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCcagc                            >  2:3321038/1‑150 (MQ=255)
                                                                                                                         ccAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCcagc                            >  2:6467262/1‑150 (MQ=255)
                                                                                                                         ccAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCcagc                            >  2:4019286/1‑150 (MQ=255)
                                                                                                                          cAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCcagca                           <  2:4933631/150‑1 (MQ=255)
                                                                                                                             aaGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTa                        <  1:2271428/150‑1 (MQ=255)
                                                                                                                               gCAGTCGTATACCGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCACAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACt                      <  2:294530/150‑1 (MQ=255)
                                                                                                                                cATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACtt                     <  1:5833799/150‑1 (MQ=255)
                                                                                                                                 aTTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTg                    >  1:3212170/1‑150 (MQ=255)
                                                                                                                                   tCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaa                  <  1:1561780/150‑1 (MQ=255)
                                                                                                                                    cGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaat                 <  2:3112661/150‑1 (MQ=255)
                                                                                                                                      tataCAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaatga               <  2:2368393/150‑1 (MQ=255)
                                                                                                                                        taCAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaatgaca             >  2:6811920/1‑150 (MQ=255)
                                                                                                                                             gTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaatgacaatgac        >  1:6299063/1‑150 (MQ=255)
                                                                                                                                             gTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGaatgacaatgac        >  1:3570783/1‑150 (MQ=255)
                                                                                                                                               ccgaccagaagaACCCCTGAGGCTTTTCACCCAAGACCCCGTTCACAAGCGCGACCTCCGCTCGACGATTCTCGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACGTGACTGACAATGACCt      <  1:4715958/144‑1 (MQ=255)
                                                                                                                                                 ttcGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGAATGACAATGACCTGc    >  1:1496999/1‑150 (MQ=255)
                                                                                                                                                  tcGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGAATGACAATGACCTGCa   >  2:4590060/1‑150 (MQ=255)
                                                                                                                                                  tcGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGAATGACAATGACCTGCa   >  1:4840149/1‑150 (MQ=255)
                                                                                                                                                  tcGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGAATGACAATGACCTGCa   <  2:899768/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCAAGCATTGCCCCGCAACTGAGGCTACGCATCACGCCAGGGGACTCGCGCAGTCTTTACGAACAAGTAGCCTGCGGAGAGCTCGATGCGGGAGTCATGGTATGTCCGCCCTTCAGCGTCCCCAAAAGCATTCGTATACAGGTTCTTCGAGAAGAACCCCTGATGCTTTTGACCAAAGACCCCGTTGACAAGCGCGACATCCGCTCGACGATTCTTGGTCGAAACCTGATTCGATATGACTCACATGCCTGGGGTGGGCAGATTGCCCAGCAGTACTTGAATGACAATGACCTGCAG  >  NC_002947/3770995‑3771291

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