Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,358,641 G→T 100% I153I (ATC→ATA ais ← putative LPS core heptose(II)‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,358,6410GT98.6% 680.6 / 0.5 215I153I (ATC→ATAaisputative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (110/102);  total (111/104)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
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 220 total aligned reads displayed.

TTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCGGACGGTATTACTG  >  CP009273/2358496‑2358785
                                                                                                                                                 |                                                                                                                                                
tttAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTat                                                                                                                                                <  1:560570/148‑1 (MQ=255)
 ttAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACtata                                                                                                                                               <  2:771753/148‑1 (MQ=255)
 ttAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACtata                                                                                                                                               <  2:3516791/148‑1 (MQ=255)
 ttAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACtata                                                                                                                                               >  1:1637360/1‑148 (MQ=255)
   aGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATAtt                                                                                                                                             <  2:2162328/148‑1 (MQ=255)
      ggTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATAttttt                                                                                                                                          <  2:3403653/148‑1 (MQ=255)
        ttAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTAt                                                                                                                                        >  1:2749209/1‑148 (MQ=255)
         tAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATc                                                                                                                                       <  2:3562097/148‑1 (MQ=255)
           aCGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAg                                                                                                                                     <  1:2860391/148‑1 (MQ=255)
            cGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAgg                                                                                                                                    <  1:183794/148‑1 (MQ=255)
              aaTTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGgcg                                                                                                                                  <  2:1423939/148‑1 (MQ=255)
                ttCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCgct                                                                                                                                >  1:3240088/1‑148 (MQ=255)
                ttCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCgct                                                                                                                                <  2:1529151/148‑1 (MQ=255)
                ttCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCgct                                                                                                                                >  1:1250892/1‑148 (MQ=255)
                   cccATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCgctttg                                                                                                                             <  1:2067146/148‑1 (MQ=255)
                   cccATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCgctttg                                                                                                                             >  1:981423/1‑148 (MQ=255)
                    ccATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCgctttgc                                                                                                                            >  2:1713670/1‑148 (MQ=255)
                        ccAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATGACTATATTTTTATCAGGCGctttgctttg                                                                                                                        >  1:91582/1‑148 (MQ=255)
                        ccAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGctttgctttg                                                                                                                        >  2:376134/1‑148 (MQ=255)
                         cAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGt                                                                                                                       <  1:2565308/148‑1 (MQ=255)
                          aGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTa                                                                                                                      >  1:1624375/1‑148 (MQ=255)
                             tAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGt                                                                                                                   <  2:1728692/148‑1 (MQ=255)
                               aaCTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTcc                                                                                                                 <  1:1308838/148‑1 (MQ=255)
                               aaCTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTcc                                                                                                                 <  1:1931763/148‑1 (MQ=255)
                                aCTTTTCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCTGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCt                                                                                                                <  2:3078115/148‑1 (MQ=255)
                                aCTTTGCCTTTTTCCACATGGATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATAAGGCGCTTTGCTTTGTAAGTCCt                                                                                                                >  2:2059092/1‑148 (MQ=255)
                                aCTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAGCTATATTTTTATCAGGCGCTGTGGTTTGTAAGTCCt                                                                                                                >  2:2525442/1‑148 (MQ=255)
                                aCTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCt                                                                                                                <  2:3077808/148‑1 (MQ=255)
                                  tttGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTa                                                                                                              >  1:3338825/1‑148 (MQ=255)
                                  tttGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTa                                                                                                              <  2:2442401/148‑1 (MQ=255)
                                  tttGCCTTTGTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTa                                                                                                              <  2:2140418/148‑1 (MQ=255)
                                    tGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAAt                                                                                                            <  2:897397/148‑1 (MQ=255)
                                       cTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGc                                                                                                         <  2:2899407/148‑1 (MQ=255)
                                       cTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGc                                                                                                         >  1:2701593/1‑148 (MQ=255)
                                         ttttCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCAc                                                                                                       >  1:1190585/1‑148 (MQ=255)
                                          tttCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACt                                                                                                      <  1:135435/148‑1 (MQ=255)
                                           ttCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACta                                                                                                     <  1:877472/148‑1 (MQ=255)
                                            tCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACtat                                                                                                    >  2:2209218/1‑148 (MQ=255)
                                              cacaTGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATaa                                                                                                  <  2:2926606/148‑1 (MQ=255)
                                                catgcatgACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAAt                                                                                                >  1:782692/1‑148 (MQ=255)
                                                  tgcatgACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAAtct                                                                                              <  1:2831191/148‑1 (MQ=255)
                                                    catgACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCa                                                                                            <  2:2808793/148‑1 (MQ=255)
                                                    catgACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCa                                                                                            <  2:1479998/148‑1 (MQ=255)
                                                    catgACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCa                                                                                            <  2:1387680/148‑1 (MQ=255)
                                                      tgACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCAtt                                                                                          >  1:2936225/1‑148 (MQ=255)
                                                       gACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTa                                                                                         >  1:1076717/1‑148 (MQ=255)
                                                             aCCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCAc                                                                                   >  2:3462639/1‑148 (MQ=255)
                                                             aCCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCAc                                                                                   >  2:866985/1‑148 (MQ=255)
                                                              ccATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACt                                                                                  >  2:3317981/1‑148 (MQ=255)
                                                               cATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTg                                                                                 <  2:1365282/148‑1 (MQ=255)
                                                                aTCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTga                                                                                <  1:159533/148‑1 (MQ=255)
                                                                  ccAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTgaag                                                                              >  2:3028016/1‑148 (MQ=255)
                                                                    aGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCATGCAATGATTTTGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTgaagaa                                                                            >  1:2381025/1‑148 (MQ=255)
                                                                      aTAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGt                                                                          >  2:3221251/1‑148 (MQ=255)
                                                                      aTAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGt                                                                          >  1:1337561/1‑148 (MQ=255)
                                                                      aTAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGt                                                                          >  1:29489/1‑148 (MQ=255)
                                                                           cAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGttt                                                                     <  1:2850605/148‑1 (MQ=255)
                                                                               tttAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCt                                                                 >  1:1501186/1‑148 (MQ=255)
                                                                                ttAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTa                                                                <  1:186389/148‑1 (MQ=255)
                                                                                ttAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTa                                                                <  1:1792099/148‑1 (MQ=255)
                                                                                    aTGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGt                                                            <  2:509507/148‑1 (MQ=255)
                                                                                      gTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCa                                                          >  1:1442740/1‑148 (MQ=255)
                                                                                       tCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCaa                                                         >  1:1734428/1‑148 (MQ=255)
                                                                                           gTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAAtttt                                                     >  2:1796475/1‑148 (MQ=255)
                                                                                           gTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAAtttt                                                     <  2:1910447/148‑1 (MQ=255)
                                                                                            tCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAAttttt                                                    >  1:2110998/1‑148 (MQ=255)
                                                                                             cACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAAtttttt                                                   >  2:778128/1‑148 (MQ=255)
                                                                                             cACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAAtttttt                                                   >  1:2680835/1‑148 (MQ=255)
                                                                                             cACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAAtttttt                                                   >  2:1183636/1‑148 (MQ=255)
                                                                                                 cTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTAccc                                               >  1:2359745/1‑148 (MQ=255)
                                                                                                    aTCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCt                                            >  2:3406477/1‑148 (MQ=255)
                                                                                                      cTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGa                                          >  2:3041092/1‑148 (MQ=255)
                                                                                                        ttAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGaaa                                        >  2:3124435/1‑148 (MQ=255)
                                                                                                           gCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAAcc                                     >  2:2960044/1‑148 (MQ=255)
                                                                                                           gCAATATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAAcc                                     >  1:2628498/1‑148 (MQ=255)
                                                                                                            cAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTAACCGCTGAAAACCa                                    >  1:3460758/1‑148 (MQ=255)
                                                                                                             aaTATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAg                                   <  2:29489/148‑1 (MQ=255)
                                                                                                             aaTATATGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAg                                   <  2:2701593/148‑1 (MQ=255)
                                                                                                               tatatGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGt                                 <  2:800359/148‑1 (MQ=255)
                                                                                                               tatatGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGt                                 <  2:2304685/148‑1 (MQ=255)
                                                                                                                 tatGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAg                               <  1:2970359/148‑1 (MQ=255)
                                                                                                                  atGTCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAGTCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGc                              <  2:1850585/148‑1 (MQ=255)
                                                                                                                     tCAGGCAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGa                           <  1:2925085/148‑1 (MQ=255)
                                                                                                                          cAATGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGaa                      <  2:2749209/148‑1 (MQ=255)
                                                                                                                             tGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGTAGAAGTCGTTTATCTGCCGTCAATTGTTTTCCCGCTGAAACCCAGGTAGCCGGCTGCATgg                   >  2:3349011/1‑148 (MQ=255)
                                                                                                                             tGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATgg                   <  2:290956/148‑1 (MQ=255)
                                                                                                                             tGATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATgg                   <  1:2209218/148‑1 (MQ=255)
                                                                                                                              gATTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGt                  <  2:1082667/148‑1 (MQ=255)
                                                                                                                               aTTATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTc                 <  2:2645967/148‑1 (MQ=255)
                                                                                                                                ttATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTcc                >  2:1624392/1‑148 (MQ=255)
                                                                                                                                 tATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCg               >  2:2045816/1‑148 (MQ=255)
                                                                                                                                 tATGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCg               >  1:3079368/1‑148 (MQ=255)
                                                                                                                                  aTGGGTGAAAATAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCgg              <  1:1767250/148‑1 (MQ=255)
                                                                                                                                         aaaaTAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCGGACGGTAt       <  2:1442740/148‑1 (MQ=255)
                                                                                                                                         aaaaTAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTAACGTCAATTTTTTAACCGCTGAAAACCAGGTAGCCGACTGAATGGTCCGGACGGTAt       >  1:708656/1‑148 (MQ=255)
                                                                                                                                          aaaTAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCGGACGGTAtt      >  2:2433108/1‑148 (MQ=255)
                                                                                                                                           aaTAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCGGACGGTATTa     >  2:618901/1‑148 (MQ=255)
                                                                                                                                           aaTAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCGGACGGTATTa     <  1:3510885/148‑1 (MQ=255)
                                                                                                                                           aaTAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCGGACGGTATTa     >  2:247167/1‑148 (MQ=255)
                                                                                                                                             tAACTATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCGGACGGTATTACt   >  1:2523254/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATGATTATGGGTGAAAATAACGATATTTTTATCAGGCGCTTTGCTTTGTAAGTCCTTAATTGCACTATAAATCTCATTACCGCACTGAAGAAGTCGTTTATCTACCGTCAATTTTTTACCCGCTGAAAACCAGGTAGCCGACTGAATGGTCCGGACGGTATTACTG  >  CP009273/2358496‑2358785

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: