Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,396,598 A→T 100% E43V (GAA→GTA)  rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,396,5980AT99.1% 712.9 / ‑2.9 215E43V (GAA→GTA) rnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (126/87);  minor base C (0/2);  total (126/89)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 215 total aligned reads displayed.

GACAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATT  >  CP009273/4396453‑4396740
                                                                                                                                                 |                                                                                                                                              
 aCAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACt                                                                                                                                             >  2:896271/1‑148 (MQ=255)
     cGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGcac                                                                                                                                         <  1:299316/148‑1 (MQ=255)
     cGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGcac                                                                                                                                         >  2:101420/1‑148 (MQ=255)
       gaggaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACAt                                                                                                                                       >  1:3085437/1‑148 (MQ=255)
       gaggaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACAt                                                                                                                                       >  1:264674/1‑148 (MQ=255)
          gaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGa                                                                                                                                    >  2:1993170/1‑148 (MQ=255)
          gaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGa                                                                                                                                    <  2:248201/148‑1 (MQ=255)
           aaCCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGaa                                                                                                                                   <  2:1105255/148‑1 (MQ=255)
           aaCCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGaa                                                                                                                                   <  2:2884050/148‑1 (MQ=255)
             ccGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAgg                                                                                                                                 >  1:3210768/1‑148 (MQ=255)
             ccGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACACGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAgg                                                                                                                                 <  1:700444/148‑1 (MQ=255)
                agaTGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGa                                                                                                                              >  2:2499500/1‑148 (MQ=255)
                  aTGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAag                                                                                                                            <  2:467897/148‑1 (MQ=255)
                   tGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAaga                                                                                                                           <  1:2694519/148‑1 (MQ=255)
                    gTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAagag                                                                                                                          <  2:406166/148‑1 (MQ=255)
                           gATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGACTTTATCCTCGAACATCTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTg                                                                                                                   >  2:358209/1‑148 (MQ=255)
                           gATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTg                                                                                                                   >  1:3315727/1‑148 (MQ=255)
                            aTCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGa                                                                                                                  >  2:2678369/1‑148 (MQ=255)
                              ccTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAg                                                                                                                >  1:44869/1‑148 (MQ=255)
                              ccTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAg                                                                                                                >  1:44839/1‑148 (MQ=255)
                              ccTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAg                                                                                                                >  1:2588845/1‑148 (MQ=255)
                               cTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAgg                                                                                                               <  1:1835565/148‑1 (MQ=255)
                                tttCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGc                                                                                                              <  1:2047553/148‑1 (MQ=255)
                                 ttCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGcc                                                                                                             >  2:353257/1‑148 (MQ=255)
                                     aGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGACTTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCg                                                                                                         >  1:919576/1‑148 (MQ=255)
                                     aGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCg                                                                                                         <  2:1302343/148‑1 (MQ=255)
                                       gAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTc                                                                                                       >  2:1512283/1‑148 (MQ=255)
                                       gAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTc                                                                                                       >  2:3120404/1‑148 (MQ=255)
                                        aaCGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcg                                                                                                      >  1:1318003/1‑148 (MQ=255)
                                        aaCGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcg                                                                                                      >  1:3387928/1‑148 (MQ=255)
                                        aaCGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcg                                                                                                      <  2:1086927/148‑1 (MQ=255)
                                          cgcgAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcgcc                                                                                                    >  1:1956469/1‑148 (MQ=255)
                                             gAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCgccgcc                                                                                                 <  2:3062971/148‑1 (MQ=255)
                                              aaGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCt                                                                                                >  2:3445882/1‑148 (MQ=255)
                                                      aaaTACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGa                                                                                        >  2:158520/1‑148 (MQ=255)
                                                       aaTACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGAt                                                                                       <  1:1665608/148‑1 (MQ=255)
                                                         tACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATgg                                                                                     <  2:875442/148‑1 (MQ=255)
                                                          aCGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGa                                                                                    >  2:1948882/1‑148 (MQ=255)
                                                           cgcgAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAg                                                                                   <  2:3435001/148‑1 (MQ=255)
                                                             cgAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgcg                                                                                 >  2:2361212/1‑148 (MQ=255)
                                                                 tcccatcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGAt                                                                             >  1:1011178/1‑148 (MQ=255)
                                                                 tcccatcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGAAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGAt                                                                             >  1:2532862/1‑148 (MQ=255)
                                                                  cccatcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATg                                                                            >  2:3276023/1‑148 (MQ=255)
                                                                    catcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGt                                                                          <  1:2427442/148‑1 (MQ=255)
                                                                       cccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCaa                                                                       <  1:718307/148‑1 (MQ=255)
                                                                       cccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCaa                                                                       <  1:1434225/148‑1 (MQ=255)
                                                                       cccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGCACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCaa                                                                       <  1:2642125/148‑1 (MQ=255)
                                                                        ccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAAc                                                                      >  2:2242121/1‑148 (MQ=255)
                                                                           aGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTgg                                                                   >  2:1982281/1‑148 (MQ=255)
                                                                           aGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTgg                                                                   >  2:2739596/1‑148 (MQ=255)
                                                                            gTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGt                                                                  >  2:1796214/1‑148 (MQ=255)
                                                                             tCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTc                                                                 <  2:803697/148‑1 (MQ=255)
                                                                                  aaTTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCAc                                                            <  2:2442956/148‑1 (MQ=255)
                                                                                   aTTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACt                                                           >  1:2372329/1‑148 (MQ=255)
                                                                                   aTTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACt                                                           <  1:1408358/148‑1 (MQ=255)
                                                                                        tCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACtcgtcg                                                      <  2:593403/148‑1 (MQ=255)
                                                                                        tCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACtcgtcg                                                      >  1:2488903/1‑148 (MQ=255)
                                                                                         ccTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTcgtcgt                                                     >  1:1172566/1‑148 (MQ=255)
                                                                                         ccTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTcgtcgt                                                     >  1:85663/1‑148 (MQ=255)
                                                                                            cGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAg                                                  <  1:2256204/148‑1 (MQ=255)
                                                                                            cGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAg                                                  <  1:462198/148‑1 (MQ=255)
                                                                                              aaCATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTg                                                >  2:313085/1‑148 (MQ=255)
                                                                                               aCATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGc                                               <  1:2297932/148‑1 (MQ=255)
                                                                                                cATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCt                                              <  2:1880612/148‑1 (MQ=255)
                                                                                                cATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCt                                              >  1:1601468/1‑148 (MQ=255)
                                                                                                cATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCt                                              >  2:733856/1‑148 (MQ=255)
                                                                                                 aTTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTa                                             >  1:207055/1‑148 (MQ=255)
                                                                                                 aTTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTa                                             >  1:2913919/1‑148 (MQ=255)
                                                                                                  tttAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTAt                                            <  1:2098172/148‑1 (MQ=255)
                                                                                                   ttAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATg                                           >  1:133484/1‑148 (MQ=255)
                                                                                                     aaCCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATgcg                                         <  1:516916/148‑1 (MQ=255)
                                                                                                     aaCCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATgcg                                         <  1:2068677/148‑1 (MQ=255)
                                                                                                       ccAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCt                                       >  2:1284605/1‑148 (MQ=255)
                                                                                                        cAAACGTGAAAAACCTGCCAGCCTAGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCATGCGTCGCCGACAGCGCGCGATGGAGCGCGATGGTCAACTGGTCTGCTCTCGTCGTCCGTGCTATGCGCTg                                      >  1:1830130/1‑148 (MQ=255)
                                                                                                        cAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTg                                      >  1:1520551/1‑148 (MQ=255)
                                                                                                        cAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTg                                      <  2:196979/148‑1 (MQ=255)
                                                                                                         aaaCGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGc                                     >  2:151497/1‑148 (MQ=255)
                                                                                                           aCGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCg                                   >  1:3484585/1‑148 (MQ=255)
                                                                                                              tGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGaa                                <  1:2255558/148‑1 (MQ=255)
                                                                                                                  aaaCCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGcc                            <  1:611191/148‑1 (MQ=255)
                                                                                                                   aaCCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCt                           >  1:1952863/1‑148 (MQ=255)
                                                                                                                      cGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGa                        <  1:3200082/148‑1 (MQ=255)
                                                                                                                      cGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGa                        <  1:1768555/148‑1 (MQ=255)
                                                                                                                      cGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGa                        <  1:2047666/148‑1 (MQ=255)
                                                                                                                       ggCCAGGCGTGATGATATGGCGGTATTACTGCACATTGATGGCTAAGAGCAGATTGAAGGCCGGCGTAGCCGCCTTCTAGCGATGTAGCGCGTTGGTCTACTGTTCTTCACTAGTCGTCAGTGCTATGCGCTGCCGGAACGCGTAGAc                       <  1:2089836/148‑1 (MQ=255)
                                                                                                                        gCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGAcc                      <  2:919576/148‑1 (MQ=255)
                                                                                                                        gCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGAcc                      >  2:3127960/1‑148 (MQ=255)
                                                                                                                         ccAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCt                     >  2:2953982/1‑148 (MQ=255)
                                                                                                                         ccAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCt                     >  1:3400646/1‑148 (MQ=255)
                                                                                                                              cGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGa                <  2:2052877/148‑1 (MQ=255)
                                                                                                                                tgatgaGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGaaa              >  1:3252187/1‑148 (MQ=255)
                                                                                                                                tgatgaGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGaaa              <  2:1601207/148‑1 (MQ=255)
                                                                                                                                 gatgaGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAg             <  1:2499500/148‑1 (MQ=255)
                                                                                                                                  atgaGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAgg            >  2:1095626/1‑148 (MQ=255)
                                                                                                                                    gaGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTa          <  1:369107/148‑1 (MQ=255)
                                                                                                                                     aGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTAc         <  2:2059483/148‑1 (MQ=255)
                                                                                                                                       cTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCg       >  1:2258285/1‑148 (MQ=255)
                                                                                                                                          gCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTa    >  1:2313069/1‑148 (MQ=255)
                                                                                                                                          gCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTa    >  2:459805/1‑148 (MQ=255)
                                                                                                                                            ggTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTAtt  <  2:256513/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GACAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATT  >  CP009273/4396453‑4396740

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: