Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,638,982 C→A 6.2% E60* (GAG→TAG)  rem ← Qin prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,638,9820CA 6.2% 669.8 / 14.0 176E60* (GAG→TAG) remQin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (99/66);  new base A (9/2);  total (108/68)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 176 total aligned reads displayed.

TGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAA  >  CP009273/1638836‑1639128
                                                                                                                                                  |                                                                                                                                                  
tGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCa                                                                                                                                                   >  1:1719332/1‑148 (MQ=255)
     tCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTc                                                                                                                                              >  1:338413/1‑148 (MQ=255)
      cGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTcg                                                                                                                                             >  1:3760226/1‑148 (MQ=255)
         ttttCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTAACCTTCgcgc                                                                                                                                          >  2:97042/1‑148 (MQ=255)
            tctcCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTAGCGCAAg                                                                                                                                       <  1:2019200/148‑1 (MQ=255)
               ccGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATa                                                                                                                                    <  2:547064/148‑1 (MQ=255)
               ccGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTAACCTTCGCGCAAGATa                                                                                                                                    >  2:1523944/1‑148 (MQ=255)
               ccGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTAACCTTCGCGCAAGATa                                                                                                                                    >  2:1523899/1‑148 (MQ=255)
                  gTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAAt                                                                                                                                 <  1:1366904/148‑1 (MQ=255)
                  gTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAAt                                                                                                                                 <  1:2666077/148‑1 (MQ=255)
                  gTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAAt                                                                                                                                 <  2:2190195/148‑1 (MQ=255)
                     gCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAAtccg                                                                                                                              >  1:2604160/1‑148 (MQ=255)
                      cACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAAtccgt                                                                                                                             <  2:2930119/148‑1 (MQ=255)
                         ggTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAAtccgtccg                                                                                                                          <  1:2447782/148‑1 (MQ=255)
                            gtTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGa                                                                                                                       >  1:3285702/1‑148 (MQ=255)
                              tACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGAcc                                                                                                                     <  1:2263211/148‑1 (MQ=255)
                               aCTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGAccc                                                                                                                    >  1:217300/1‑148 (MQ=255)
                                  cAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAg                                                                                                                 >  1:2138307/1‑148 (MQ=255)
                                     cGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCAt                                                                                                              >  2:3095789/1‑148 (MQ=255)
                                      gggAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATc                                                                                                             <  2:2977842/148‑1 (MQ=255)
                                           ttCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCAGCCACaa                                                                                                        >  2:2112403/1‑148 (MQ=255)
                                                ccccGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAaatta                                                                                                   <  1:2512974/148‑1 (MQ=255)
                                                  ccGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAaattaat                                                                                                 >  1:2139090/1‑148 (MQ=255)
                                                  ccGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAaattaat                                                                                                 <  2:3307706/148‑1 (MQ=255)
                                                    gcgcAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAattaatta                                                                                               <  1:3439355/148‑1 (MQ=255)
                                                    gcgcAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAattaatta                                                                                               <  1:3453492/148‑1 (MQ=255)
                                                       cAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACg                                                                                            >  1:1753229/1‑148 (MQ=255)
                                                        aaGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGc                                                                                           <  1:3679072/148‑1 (MQ=255)
                                                        aaGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCGTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGc                                                                                           <  1:1493191/148‑1 (MQ=255)
                                                         aGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCa                                                                                          >  1:1235738/1‑148 (MQ=255)
                                                         aGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTAACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCa                                                                                          >  1:44183/1‑148 (MQ=255)
                                                          gATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAg                                                                                         >  1:3236502/1‑148 (MQ=255)
                                                          gATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAg                                                                                         >  2:567840/1‑148 (MQ=255)
                                                            ttGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGcc                                                                                       <  1:1017993/148‑1 (MQ=255)
                                                             tGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGccc                                                                                      >  2:2604568/1‑148 (MQ=255)
                                                               tAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCa                                                                                    >  2:2864199/1‑148 (MQ=255)
                                                               tAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCa                                                                                    >  2:2725207/1‑148 (MQ=255)
                                                                  aTGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTc                                                                                 >  2:1659478/1‑148 (MQ=255)
                                                                  aTGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCACTTTCAAGGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTc                                                                                 <  1:857851/148‑1 (MQ=255)
                                                                   tGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCt                                                                                >  1:3757925/1‑148 (MQ=255)
                                                                   tGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCt                                                                                >  2:761714/1‑148 (MQ=255)
                                                                   tGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCt                                                                                <  1:1745064/148‑1 (MQ=255)
                                                                   tGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCt                                                                                >  1:2881581/1‑148 (MQ=255)
                                                                    gAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCtt                                                                               >  1:1408703/1‑148 (MQ=255)
                                                                    gAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCtt                                                                               >  1:1409573/1‑148 (MQ=255)
                                                                     aGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCttt                                                                              >  1:374382/1‑148 (MQ=255)
                                                                        ttATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGc                                                                           >  1:3237690/1‑148 (MQ=255)
                                                                           ttCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCaaa                                                                        >  1:1468134/1‑148 (MQ=255)
                                                                               ccTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAAcaca                                                                    >  1:3463956/1‑148 (MQ=255)
                                                                                 tGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATa                                                                  >  1:100988/1‑148 (MQ=255)
                                                                                  gAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATaa                                                                 >  1:1948778/1‑148 (MQ=255)
                                                                                  gAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATaa                                                                 <  1:1671063/148‑1 (MQ=255)
                                                                                  gAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATaa                                                                 >  1:3714243/1‑148 (MQ=255)
                                                                                  gAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATaa                                                                 <  1:749222/148‑1 (MQ=255)
                                                                                   aaaaaGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAAt                                                                <  2:1730094/148‑1 (MQ=255)
                                                                                     aaaGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGc                                                              <  2:3168602/148‑1 (MQ=255)
                                                                                     aaaGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGc                                                              <  1:3096375/148‑1 (MQ=255)
                                                                                     aaaGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGc                                                              <  1:1194371/148‑1 (MQ=255)
                                                                                       aGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCgg                                                            <  1:2405664/148‑1 (MQ=255)
                                                                                            aaaaGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATg                                                       >  1:1501326/1‑148 (MQ=255)
                                                                                             aaaGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGa                                                      <  2:3539386/148‑1 (MQ=255)
                                                                                             aaaGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGa                                                      <  2:667934/148‑1 (MQ=255)
                                                                                               aGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACt                                                    >  1:2292084/1‑148 (MQ=255)
                                                                                               aGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACt                                                    >  1:2626476/1‑148 (MQ=255)
                                                                                                gCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTa                                                   >  2:1087338/1‑148 (MQ=255)
                                                                                                 ccAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTaa                                                  >  1:1207981/1‑148 (MQ=255)
                                                                                                 ccAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTaa                                                  <  2:264645/148‑1 (MQ=255)
                                                                                                 ccAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTaa                                                  >  1:1208162/1‑148 (MQ=255)
                                                                                                   aGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACa                                                <  2:1210775/148‑1 (MQ=255)
                                                                                                      ttttATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGa                                             >  1:3514879/1‑148 (MQ=255)
                                                                                                      ttttATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGa                                             >  2:1234858/1‑148 (MQ=255)
                                                                                                         tATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGAttt                                          <  2:170351/148‑1 (MQ=255)
                                                                                                               gATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCAttt                                    <  1:1659478/148‑1 (MQ=255)
                                                                                                                aTCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCAtttt                                   >  1:1143574/1‑148 (MQ=255)
                                                                                                                  ctctTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCa                                 <  1:3532669/148‑1 (MQ=255)
                                                                                                                  ctctTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGAAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCa                                 <  2:670219/148‑1 (MQ=255)
                                                                                                                  ctctTTCAATGCCTGTAATGAAGTGACAAACTAACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCa                                 >  2:3323984/1‑148 (MQ=255)
                                                                                                                    ctTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAAc                               >  1:3128966/1‑148 (MQ=255)
                                                                                                                    ctTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAAc                               >  2:1401062/1‑148 (MQ=255)
                                                                                                                      ttCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGc                             >  2:2438891/1‑148 (MQ=255)
                                                                                                                      ttCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGc                             >  1:407427/1‑148 (MQ=255)
                                                                                                                        cAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGcac                           >  1:66031/1‑148 (MQ=255)
                                                                                                                          aTGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGcacac                         >  2:2391486/1‑148 (MQ=255)
                                                                                                                           tGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGcacaca                        >  2:2957803/1‑148 (MQ=255)
                                                                                                                            gCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAg                       >  2:971857/1‑148 (MQ=255)
                                                                                                                               tGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCAt                    <  2:2332347/148‑1 (MQ=255)
                                                                                                                               tGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCAt                    <  2:1719332/148‑1 (MQ=255)
                                                                                                                                 tAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATag                  >  1:2697633/1‑148 (MQ=255)
                                                                                                                                  aaTGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATaga                 <  2:3153953/148‑1 (MQ=255)
                                                                                                                                  aaTGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATaga                 >  1:662315/1‑148 (MQ=255)
                                                                                                                                   aTGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATagag                <  2:2734647/148‑1 (MQ=255)
                                                                                                                                   aTGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATagag                <  2:2734629/148‑1 (MQ=255)
                                                                                                                                   aTGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATagag                >  2:1783346/1‑148 (MQ=255)
                                                                                                                                      aaGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGcc             >  2:255027/1‑148 (MQ=255)
                                                                                                                                        gTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCAt           >  1:3399245/1‑148 (MQ=255)
                                                                                                                                            cAAACTCACCTTCTCTCAAGATAAATCCGTCCGTGACCAGAGCATCCACAAAATTAATTAACGCAGCCACATTCTTTCGCAAACACATAATGCGGTAATGACTAACAATATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTc       <  2:389106/148‑1 (MQ=255)
                                                                                                                                            cAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTc       <  2:389735/148‑1 (MQ=255)
                                                                                                                                                cTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACaa   >  2:2016721/1‑148 (MQ=255)
                                                                                                                                                 tCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACaaa  >  1:2997703/1‑148 (MQ=255)
                                                                                                                                                 tCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACaaa  >  1:2998047/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGCCATCGGTTTTCTCCGGTGGCACGGTGTTACTCAGCGGGAGTTCAGCCCCGCGCAAGATTGTAGATGAGTTTATTCTCCTGAAAAAGCAGAAAAGCCAGCTTTTATTCCGATCTCTTTCAATGCCTGTAATGAAGTGACAAACTCACCTTCGCGCAAGATAAATCCGTCCGTGACCCGAGCATCCACAAAATTAATTAACGCAGCCCCATTCTTTCGCAAACACATAATGCGGTAATGACTAACAAGATTTCCATTTTCAACGCACACAGCATAGAGGCCATCTTCACAAA  >  CP009273/1638836‑1639128

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: