Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,478,971 G→A 100% E552K (GAG→AAG)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,478,9710GA96.0% 702.6 / 15.2 226E552K (GAG→AAG) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base A (114/103);  total (119/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 227 total aligned reads displayed.

CGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACCC  >  CP009273/1478827‑1479115
                                                                                                                                                |                                                                                                                                                
cgTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGt                                                                                                                                               <  1:1066444/148‑1 (MQ=255)
cgTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGt                                                                                                                                               <  1:1081254/148‑1 (MQ=255)
cgTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGt                                                                                                                                               <  1:2022935/148‑1 (MQ=255)
   aTGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTACCACGACAAAAAGTCTg                                                                                                                                            >  2:2374918/1‑148 (MQ=255)
         cTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTtc                                                                                                                                      <  1:2074490/148‑1 (MQ=255)
         cTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtc                                                                                                                                      <  1:1883855/148‑1 (MQ=255)
          tGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTtct                                                                                                                                     <  2:383652/148‑1 (MQ=255)
              aGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTcgcg                                                                                                                                 >  2:253760/1‑148 (MQ=255)
               gcgcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGt                                                                                                                                <  2:738356/148‑1 (MQ=255)
               gcgcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGt                                                                                                                                <  2:1595126/148‑1 (MQ=255)
                cgcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGtt                                                                                                                               <  2:3149218/148‑1 (MQ=255)
                  cAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTtg                                                                                                                             >  2:2320190/1‑148 (MQ=255)
                  cAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTtg                                                                                                                             <  2:1171069/148‑1 (MQ=255)
                   aaaaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCgtgtgt                                                                                                                            >  1:528555/1‑148 (MQ=255)
                    aaaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTtgtg                                                                                                                           >  1:2778458/1‑148 (MQ=255)
                      aaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGaa                                                                                                                         <  2:2402072/148‑1 (MQ=255)
                       aCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAAt                                                                                                                        >  2:881191/1‑148 (MQ=255)
                       aCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAAt                                                                                                                        >  2:2662591/1‑148 (MQ=255)
                        cATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATc                                                                                                                       >  1:1018451/1‑148 (MQ=255)
                           ggCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCtgt                                                                                                                    >  1:3497113/1‑148 (MQ=255)
                            gCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCtgtg                                                                                                                   <  1:1014254/148‑1 (MQ=255)
                               gcgtgcgtGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGaa                                                                                                                <  2:822457/148‑1 (MQ=255)
                                   gcgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTAt                                                                                                            <  1:1815558/148‑1 (MQ=255)
                                   gcgtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTAt                                                                                                            <  2:2640676/148‑1 (MQ=255)
                                          gCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCg                                                                                                     >  1:3183671/1‑148 (MQ=255)
                                             atgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagc                                                                                                  >  2:546547/1‑148 (MQ=255)
                                             atgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagc                                                                                                  <  1:1341299/148‑1 (MQ=255)
                                             atgatTATCACGTCCGCGCTCTCCAGTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagc                                                                                                  <  2:1820261/148‑1 (MQ=255)
                                              tgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagca                                                                                                 >  2:2242049/1‑148 (MQ=255)
                                              tgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagca                                                                                                 >  2:425151/1‑148 (MQ=255)
                                              tgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagca                                                                                                 >  1:3035559/1‑148 (MQ=255)
                                               gatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagcag                                                                                                <  2:1483977/148‑1 (MQ=255)
                                               gatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTAGCGTTTGTGAATCTGTGGAATTATCTTGGCGtgcag                                                                                                >  1:800889/1‑148 (MQ=255)
                                                 ttATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAgcagca                                                                                              >  2:1761267/1‑148 (MQ=255)
                                                 ttATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAgcagca                                                                                              >  2:499187/1‑148 (MQ=255)
                                                   aTCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCaaa                                                                                            <  1:261266/148‑1 (MQ=255)
                                                   aTCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCaaa                                                                                            <  2:114124/148‑1 (MQ=255)
                                                    tCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAATCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAATGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCaaaa                                                                                           <  1:2756864/148‑1 (MQ=255)
                                                     cACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAg                                                                                          <  1:3305653/148‑1 (MQ=255)
                                                     cACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAg                                                                                          >  2:2269904/1‑148 (MQ=255)
                                                       cGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGgc                                                                                        >  2:1950891/1‑148 (MQ=255)
                                                       cGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGgc                                                                                        >  2:154134/1‑148 (MQ=255)
                                                        gTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGgcg                                                                                       <  1:2346372/148‑1 (MQ=255)
                                                          ccGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCt                                                                                     >  2:1855316/1‑148 (MQ=255)
                                                              gcTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttct                                                                                 <  1:2269425/148‑1 (MQ=255)
                                                                 ctcCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTACACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttgttac                                                                              >  1:1210599/1‑146 (MQ=255)
                                                                      ttCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGc                                                                         >  2:3379718/1‑148 (MQ=255)
                                                                       tCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGcc                                                                        >  2:2748217/1‑148 (MQ=255)
                                                                        cAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGACTTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCt                                                                       >  1:3510122/1‑148 (MQ=255)
                                                                         aGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCtt                                                                      <  2:3198483/148‑1 (MQ=255)
                                                                           gATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTcc                                                                    >  2:1113839/1‑148 (MQ=255)
                                                                             tCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGt                                                                  >  1:2910049/1‑148 (MQ=255)
                                                                                  gcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCt                                                             <  1:49862/148‑1 (MQ=255)
                                                                                  gcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCt                                                             <  1:1850997/148‑1 (MQ=255)
                                                                                  gcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCt                                                             <  1:2866531/148‑1 (MQ=255)
                                                                                      aCGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTgtgt                                                         >  2:1875697/1‑148 (MQ=255)
                                                                                          ccGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACACCTTCCGTCGCCTGTGTCGTa                                                     <  1:1702693/148‑1 (MQ=255)
                                                                                          ccGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACACCTTCCGTCGCCTGTGTCGTa                                                     <  1:3460411/148‑1 (MQ=255)
                                                                                             aTGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCg                                                  <  1:3427470/148‑1 (MQ=255)
                                                                                             aTGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCg                                                  >  2:226747/1‑148 (MQ=255)
                                                                                             aTGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCg                                                  <  1:1240090/148‑1 (MQ=255)
                                                                                              tGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGa                                                 >  1:3480136/1‑148 (MQ=255)
                                                                                              tGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGa                                                 <  1:2118300/148‑1 (MQ=255)
                                                                                                gACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGAtt                                               <  1:854299/148‑1 (MQ=255)
                                                                                                   aaaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtc                                            <  1:2531644/148‑1 (MQ=255)
                                                                                                   aaaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtc                                            >  1:283154/1‑148 (MQ=255)
                                                                                                   aaaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtc                                            <  2:3492023/148‑1 (MQ=255)
                                                                                                    aaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtct                                           >  2:453863/1‑148 (MQ=255)
                                                                                                    aaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtct                                           >  2:3026993/1‑148 (MQ=255)
                                                                                                    aaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtct                                           >  2:2070466/1‑148 (MQ=255)
                                                                                                    aaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtct                                           >  1:1382350/1‑148 (MQ=255)
                                                                                                    aaCAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtct                                           >  1:900777/1‑148 (MQ=255)
                                                                                                    aaCAGCAGACATCGGACGAAAAACATCGTCGCTTCCACGACAAAAATTCTGACTTTCTCGCGTGTGTGAATCTGTGGAAATATCTTTGCGAGCAGCAAAAGGCGCTTTATTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTAtct                                           >  1:3470409/1‑148 (MQ=255)
                                                                                                        gcagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAAc                                       >  1:2917054/1‑148 (MQ=255)
                                                                                                        gcagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAAc                                       >  1:3260085/1‑148 (MQ=255)
                                                                                                         cagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                      >  2:2918387/1‑148 (MQ=255)
                                                                                                         cagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                      <  1:2338405/148‑1 (MQ=255)
                                                                                                         cagGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACt                                      <  1:1761267/148‑1 (MQ=255)
                                                                                                            gCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATc                                   <  2:3265534/148‑1 (MQ=255)
                                                                                                            gCATCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATc                                   <  1:1516435/148‑1 (MQ=255)
                                                                                                               tCGGACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgc                                >  1:1184293/1‑148 (MQ=255)
                                                                                                                  gACGAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcg                             >  1:2468489/1‑148 (MQ=255)
                                                                                                                     gAAAAACATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgc                          >  1:387287/1‑148 (MQ=255)
                                                                                                                      aaaaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcg                         >  1:2643247/1‑148 (MQ=255)
                                                                                                                      aaaaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCtgcgcgtgcg                         >  1:61151/1‑148 (MQ=255)
                                                                                                                       aaaaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTAGCGTTTTAGAAGCTGTGGAATTATCTTGGCGAGCATCAAAAGGCGCTTTCTTCCAAAGCCTACCGTCGCCTGTGTCGTAACGATTATCTCAACAATCtgcgcgtgcgc                        >  1:2966953/1‑148 (MQ=255)
                                                                                                                         aaCATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGa                      <  1:711004/148‑1 (MQ=255)
                                                                                                                           cATCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAAt                    <  2:926434/148‑1 (MQ=255)
                                                                                                                            aTCGTCGCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATg                   <  2:1018451/148‑1 (MQ=255)
                                                                                                                             tcgtcgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATgg                  <  2:1355616/148‑1 (MQ=255)
                                                                                                                                 cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGACTGGCCgg              >  1:1956545/1‑148 (MQ=255)
                                                                                                                                 cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAgg              <  2:2099374/148‑1 (MQ=255)
                                                                                                                                 cgCTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAgg              >  1:2224592/1‑148 (MQ=255)
                                                                                                                                   cTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGat            <  2:3520588/148‑1 (MQ=255)
                                                                                                                                   cTTCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGat            <  1:392166/148‑1 (MQ=255)
                                                                                                                                    ttCCACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGata           >  1:442342/1‑148 (MQ=255)
                                                                                                                                      ccACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATc         >  2:204041/1‑148 (MQ=255)
                                                                                                                                       cACGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCt        >  1:852367/1‑148 (MQ=255)
                                                                                                                                         cGACAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTac      >  2:3414400/1‑148 (MQ=255)
                                                                                                                                            cAAAAAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACAcc   >  2:2194818/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGCCTGTGTCGTACCGATTATCTCAACTATCTGCGCGTGCGCGAATGGCAGGATATCTACACCC  >  CP009273/1478827‑1479115

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: ‑