Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 337,720 G→T 100% E295D (GAG→GAT yahJ → putative metallo‑dependent hydrolase domain deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273337,7200GT98.0% 713.8 / 3.1 233E295D (GAG→GATyahJputative metallo‑dependent hydrolase domain deaminase
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (134/93);  total (135/98)
Fisher's exact test for biased strand distribution p-value = 8.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 236 total aligned reads displayed.

CCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGC  >  CP009273/337574‑337865
                                                                                                                                                  |                                                                                                                                                 
   aCCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAAc                                                                                                                                              >  2:2258351/1‑149 (MQ=255)
   aCCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAAc                                                                                                                                              <  2:2048788/149‑1 (MQ=255)
   aCCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAAc                                                                                                                                              <  1:3790546/149‑1 (MQ=255)
         gttgttGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCcaca                                                                                                                                        <  1:3326612/149‑1 (MQ=255)
         gtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAACCGGTAGATAAAACGCcaca                                                                                                                                        >  2:520004/1‑149 (MQ=255)
           tGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGc                                                                                                                                      <  2:575367/149‑1 (MQ=255)
           tGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGc                                                                                                                                      >  1:4099472/1‑149 (MQ=255)
             aTGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTg                                                                                                                                    >  1:2716227/1‑149 (MQ=255)
               ggCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGaa                                                                                                                                  >  2:3197532/1‑149 (MQ=255)
               ggCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGaa                                                                                                                                  <  2:3902584/149‑1 (MQ=255)
               ggCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGaa                                                                                                                                  <  1:2062536/149‑1 (MQ=255)
                   cgATGTAAAAATCCCTCGACACATTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGGGTGGCAGCCATAAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGc                                                                                                                              <  1:457767/149‑1 (MQ=255)
                   cgATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGGAACGGTAGATAAAACGCCACAGCTGAAGGGc                                                                                                                              <  2:2102999/149‑1 (MQ=255)
                   cgATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGc                                                                                                                              <  2:3340627/149‑1 (MQ=255)
                   cgATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGc                                                                                                                              <  1:2526980/149‑1 (MQ=255)
                    gATGTAAAAATCCCTCGACACATTTTTCCAGATGGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAACACGCCACAGCTGAAGGGCa                                                                                                                             <  1:3136825/149‑1 (MQ=255)
                    gATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCa                                                                                                                             >  2:1334549/1‑149 (MQ=255)
                     aTGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCaa                                                                                                                            <  2:1658084/149‑1 (MQ=255)
                     aTGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCaa                                                                                                                            <  2:2623900/149‑1 (MQ=255)
                        gAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCt                                                                                                                         >  1:2736440/1‑149 (MQ=255)
                        gAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCt                                                                                                                         >  1:2477928/1‑149 (MQ=255)
                          aaaaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGa                                                                                                                       >  1:2594807/1‑149 (MQ=255)
                            aaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGAcc                                                                                                                     >  2:2637738/1‑149 (MQ=255)
                            aaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCAGCACTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGAcc                                                                                                                     >  2:3373463/1‑149 (MQ=255)
                             aTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCa                                                                                                                    >  1:872993/1‑149 (MQ=255)
                             aTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCa                                                                                                                    <  2:2276955/149‑1 (MQ=255)
                             aTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCa                                                                                                                    >  1:3946974/1‑149 (MQ=255)
                              tCCCTCGACACCTTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCAt                                                                                                                   <  1:3691352/149‑1 (MQ=255)
                               cccTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATc                                                                                                                  <  2:1349017/149‑1 (MQ=255)
                                  tCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGt                                                                                                               >  2:1189653/1‑149 (MQ=255)
                                       acCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGc                                                                                                          >  2:186627/1‑149 (MQ=255)
                                       acCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGCAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGc                                                                                                          <  1:318517/149‑1 (MQ=255)
                                        ccATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGcc                                                                                                         >  2:2137635/1‑149 (MQ=255)
                                         cATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                        >  2:1690010/1‑149 (MQ=255)
                                         cATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                        >  2:3979148/1‑149 (MQ=255)
                                         cATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                        >  1:3239539/1‑149 (MQ=255)
                                          aTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCtt                                                                                                       >  1:1712359/1‑149 (MQ=255)
                                          aTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCtt                                                                                                       <  2:2593967/149‑1 (MQ=255)
                                           tGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCttt                                                                                                      >  1:2950336/1‑149 (MQ=255)
                                                   aTTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGc                                                                                              <  1:811123/149‑1 (MQ=255)
                                                          tGGACTACCACAACGGCGTCGATATACACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTc                                                                                       <  1:1403580/149‑1 (MQ=255)
                                                             aCTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAAc                                                                                    <  2:3068583/149‑1 (MQ=255)
                                                             aCTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAAc                                                                                    <  2:2752242/149‑1 (MQ=255)
                                                             aCTACGACAAAGGCGTCGATATTCACATTCACGAAAACACTCCGGCACGCGAGGCAGCCATAAAATATATGGTTGAAACGGAAGATAAAACACAACAGATCAAGTGCAACCAGACAATAATTCACGCCTTTGCGCTTGCAACGCTCaaa                                                                                    >  1:949211/1‑148 (MQ=255)
                                                              cTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACg                                                                                   >  2:966150/1‑149 (MQ=255)
                                                                aCGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAg                                                                                 <  1:2697563/149‑1 (MQ=255)
                                                                   aCAAAGGCGTCGATATTCCCCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCaa                                                                              <  1:3472908/149‑1 (MQ=255)
                                                                     aaaGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACa                                                                            >  1:1865190/1‑149 (MQ=255)
                                                                      aaGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAg                                                                           <  1:2413327/149‑1 (MQ=255)
                                                                      aaGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAg                                                                           >  2:273491/1‑149 (MQ=255)
                                                                         gCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTa                                                                        <  1:211320/149‑1 (MQ=255)
                                                                            tCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGAt                                                                     >  2:62282/1‑149 (MQ=255)
                                                                             cGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATg                                                                    <  2:2379717/149‑1 (MQ=255)
                                                                                 atTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACt                                                                >  2:4079838/1‑149 (MQ=255)
                                                                                  ttCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTg                                                               >  1:3552595/1‑149 (MQ=255)
                                                                                    cACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGc                                                             <  1:2841425/149‑1 (MQ=255)
                                                                                      ccTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGa                                                           >  1:3937698/1‑149 (MQ=255)
                                                                                         gCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAAcc                                                        >  2:3993886/1‑149 (MQ=255)
                                                                                         gCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAAcc                                                        >  1:2059090/1‑149 (MQ=255)
                                                                                            cGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATCTGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGa                                                     <  2:502568/149‑1 (MQ=255)
                                                                                            cGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGa                                                     <  1:126008/149‑1 (MQ=255)
                                                                                              aaaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtg                                                   >  1:2058225/1‑149 (MQ=255)
                                                                                              aaaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtg                                                   >  1:2058292/1‑149 (MQ=255)
                                                                                               aaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtgg                                                  >  2:4087000/1‑149 (MQ=255)
                                                                                               aaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtgg                                                  <  1:601528/149‑1 (MQ=255)
                                                                                                aCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtggt                                                 >  2:3528883/1‑149 (MQ=255)
                                                                                                aCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtggt                                                 >  2:1382868/1‑149 (MQ=255)
                                                                                                   aCTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATggtggt                                              >  2:2308760/1‑149 (MQ=255)
                                                                                                           aGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACaaa                                      >  1:971824/1‑149 (MQ=255)
                                                                                                            ggCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                     >  1:1603111/1‑149 (MQ=255)
                                                                                                                 ggCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTa                                <  1:3781965/149‑1 (MQ=255)
                                                                                                                  gCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTAt                               <  1:1130310/149‑1 (MQ=255)
                                                                                                                  gCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTAt                               >  2:1678026/1‑149 (MQ=255)
                                                                                                                  gCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTAt                               >  1:3044493/1‑149 (MQ=255)
                                                                                                                  gCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTAt                               <  1:1351882/149‑1 (MQ=255)
                                                                                                                   cAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATc                              <  2:87864/149‑1 (MQ=255)
                                                                                                                        aTCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTc                         <  1:60548/149‑1 (MQ=255)
                                                                                                                         tCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCg                        >  1:1132103/1‑149 (MQ=255)
                                                                                                                             ttATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACgg                    <  2:2716227/149‑1 (MQ=255)
                                                                                                                             ttATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACgg                    <  2:596628/149‑1 (MQ=255)
                                                                                                                             ttATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACgg                    >  2:2903579/1‑149 (MQ=255)
                                                                                                                             ttATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACgg                    <  2:1009990/149‑1 (MQ=255)
                                                                                                                              tattttgTTGAAACGGTATATAATACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGt                   <  2:3476321/143‑1 (MQ=255)
                                                                                                                              tataTGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGt                   <  2:1426543/149‑1 (MQ=255)
                                                                                                                               atatGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTg                  <  1:1335339/149‑1 (MQ=255)
                                                                                                                                tatGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGc                 <  2:3631758/149‑1 (MQ=255)
                                                                                                                                tatGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGc                 <  2:1585611/149‑1 (MQ=255)
                                                                                                                                 atGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGcc                <  1:3264996/149‑1 (MQ=255)
                                                                                                                                  tGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCg               >  1:555056/1‑149 (MQ=255)
                                                                                                                                   ggTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACGAATTTCTATCGCCTCGACGGTGCCGa              <  2:2466236/149‑1 (MQ=255)
                                                                                                                                   ggTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGa              >  1:2822122/1‑149 (MQ=255)
                                                                                                                                       gAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTgg          >  2:1981998/1‑149 (MQ=255)
                                                                                                                                        aaaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACTCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGc         <  2:634942/149‑1 (MQ=255)
                                                                                                                                        aaaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGc         <  1:2316526/149‑1 (MQ=255)
                                                                                                                                         aaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCa        >  2:2475875/1‑149 (MQ=255)
                                                                                                                                         aaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCa        <  1:2167501/149‑1 (MQ=255)
                                                                                                                                             gTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCACCAAATTTCTATCGCCTCGACGGTGCCGCTTGGCACGCt    >  2:2687125/1‑149 (MQ=255)
                                                                                                                                             gTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCt    >  1:4076598/1‑149 (MQ=255)
                                                                                                                                               agatAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGc  <  2:310341/149‑1 (MQ=255)
                                                                                                                                               agatAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGc  <  2:636419/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTGC  >  CP009273/337574‑337865

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: