Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,044,791 C→A 15.8% R29S (CGT→AGT)  yihI → activator of Der GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,044,7910CA15.8% 631.1 / 64.8 213R29S (CGT→AGT) yihIactivator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base C (116/63);  new base A (23/11);  total (139/74)
Fisher's exact test for biased strand distribution p-value = 8.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 215 total aligned reads displayed.

TTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTG  >  CP009273/4044648‑4044938
                                                                                                                                               |                                                                                                                                                   
     gTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCa                                                                                                                                            <  2:3059546/148‑1 (MQ=255)
      ttAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCaa                                                                                                                                           <  2:2204616/148‑1 (MQ=255)
      ttAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCaa                                                                                                                                           >  1:2436327/1‑148 (MQ=255)
       gaGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAg                                                                                                                                          <  2:2880298/147‑1 (MQ=255)
         gATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCg                                                                                                                                        >  2:2761125/1‑148 (MQ=255)
          aTTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGt                                                                                                                                       >  2:574293/1‑148 (MQ=255)
            tCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCa                                                                                                                                     >  1:915104/1‑148 (MQ=255)
            tCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCa                                                                                                                                     >  1:973982/1‑148 (MQ=255)
                  gcAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGaaaaa                                                                                                                               >  2:1791614/1‑148 (MQ=255)
                   cAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGaaaaaa                                                                                                                              <  2:1495651/148‑1 (MQ=255)
                    aaaTCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAAc                                                                                                                             >  2:414700/1‑148 (MQ=255)
                     aaTCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAAcg                                                                                                                            >  1:2094072/1‑148 (MQ=255)
                       tCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGACGCTCGTGATCGCAAGCGTCAGAAAAAAcgtc                                                                                                                          >  2:464444/1‑148 (MQ=255)
                       tCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAAcgtc                                                                                                                          >  2:860027/1‑148 (MQ=255)
                          ataCTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTg                                                                                                                       >  1:2949429/1‑148 (MQ=255)
                             cTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTc                                                                                                                    >  1:2053363/1‑148 (MQ=255)
                              ttaaTAAAACGGCTAAAAGACGAACAATTATTAAACCATCATCTTAAAACTCACGCAGCAAAGGTCATGAAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCa                                                                                                                   <  1:405990/146‑1 (MQ=255)
                              tGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCa                                                                                                                   <  1:406468/148‑1 (MQ=255)
                               gAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCAt                                                                                                                  >  2:118939/1‑148 (MQ=255)
                                   aaaaCGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATgcgc                                                                                                              >  1:742935/1‑148 (MQ=255)
                                     aaCGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCg                                                                                                            >  2:2554054/1‑148 (MQ=255)
                                     aaCGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCg                                                                                                            >  2:388861/1‑148 (MQ=255)
                                          cTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAg                                                                                                       >  2:2456360/1‑148 (MQ=255)
                                            aaaaGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGcc                                                                                                     >  2:2991272/1‑148 (MQ=255)
                                            aaaaGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGcc                                                                                                     >  1:3210556/1‑148 (MQ=255)
                                            aaaaGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGcc                                                                                                     >  2:828948/1‑148 (MQ=255)
                                              aaGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCcgc                                                                                                   >  2:3149706/1‑148 (MQ=255)
                                               aGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCcgcg                                                                                                  >  2:928833/1‑148 (MQ=255)
                                                 aCGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCa                                                                                                >  1:3032193/1‑148 (MQ=255)
                                                 aCGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCa                                                                                                >  2:1973610/1‑148 (MQ=255)
                                                  cGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAg                                                                                               <  1:874388/148‑1 (MQ=255)
                                                    aaCCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAgcg                                                                                             <  2:110510/148‑1 (MQ=255)
                                                     aCCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAgcgg                                                                                            >  1:2278750/1‑148 (MQ=255)
                                                      ccATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAgcggg                                                                                           >  1:1243280/1‑148 (MQ=255)
                                                         ttATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAgcgggcgg                                                                                        >  1:2360110/1‑148 (MQ=255)
                                                         ttATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAgcgggcgg                                                                                        >  1:2322108/1‑148 (MQ=255)
                                                         ttATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAgcgggcgg                                                                                        >  1:661691/1‑148 (MQ=255)
                                                          tATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGa                                                                                       >  1:383500/1‑148 (MQ=255)
                                                          tATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGa                                                                                       >  1:2606796/1‑148 (MQ=255)
                                                            tGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGaca                                                                                     >  2:2146514/1‑146 (MQ=255)
                                                            tGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGaaa                                                                                     <  2:170409/148‑1 (MQ=255)
                                                              aaaCCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAaca                                                                                   >  2:2261784/1‑148 (MQ=255)
                                                              aaaCCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAaca                                                                                   >  1:767652/1‑148 (MQ=255)
                                                              aaaCCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAaca                                                                                   >  2:451388/1‑148 (MQ=255)
                                                                aCCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAAcacc                                                                                 >  1:2585800/1‑148 (MQ=255)
                                                                 ccATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAAcacca                                                                                >  1:156940/1‑148 (MQ=255)
                                                                 ccATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAAcacca                                                                                <  2:1348710/148‑1 (MQ=255)
                                                                  catcatCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAAcaccac                                                                               >  1:1598872/1‑148 (MQ=255)
                                                                   atcatcTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACg                                                                              >  1:968061/1‑148 (MQ=255)
                                                                          tCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGgca                                                                       >  2:3247638/1‑148 (MQ=255)
                                                                                  aCGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGc                                                               >  2:1544622/1‑148 (MQ=255)
                                                                                  aCGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGc                                                               <  1:3011129/148‑1 (MQ=255)
                                                                                          aaGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGc                                                       >  1:1841803/1‑148 (MQ=255)
                                                                                          aaGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGc                                                       <  2:1035793/148‑1 (MQ=255)
                                                                                          aaGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGc                                                       <  1:483749/148‑1 (MQ=255)
                                                                                           aGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCa                                                      >  2:329668/1‑148 (MQ=255)
                                                                                           aGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCa                                                      >  1:1187731/1‑148 (MQ=255)
                                                                                           aGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCa                                                      <  2:2045174/148‑1 (MQ=255)
                                                                                            ggTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCAc                                                     >  1:2205205/1‑148 (MQ=255)
                                                                                            ggTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCAc                                                     >  1:3212613/1‑148 (MQ=255)
                                                                                             gTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCAcc                                                    >  2:1904911/1‑148 (MQ=255)
                                                                                                aTGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCaaa                                                 <  2:2608700/148‑1 (MQ=255)
                                                                                                aTGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCaaa                                                 <  2:303866/148‑1 (MQ=255)
                                                                                                   cAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGACAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGa                                              >  2:2774129/1‑148 (MQ=255)
                                                                                                    aaaaGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGAt                                             >  1:2716499/1‑148 (MQ=255)
                                                                                                        gcgcGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCAc                                         >  2:1869627/1‑148 (MQ=255)
                                                                                                           cgTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTa                                      <  2:915104/148‑1 (MQ=255)
                                                                                                                   aaCACGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGt                              >  1:1761144/1‑148 (MQ=255)
                                                                                                                    acacGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTa                             >  1:2396641/1‑148 (MQ=255)
                                                                                                                    acacGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTa                             >  1:2396687/1‑148 (MQ=255)
                                                                                                                     cacGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTaa                            >  2:3073211/1‑148 (MQ=255)
                                                                                                                     cacGCGAGGAGCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTaa                            >  1:1011631/1‑148 (MQ=255)
                                                                                                                      acGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTaaa                           <  2:2278750/148‑1 (MQ=255)
                                                                                                                       cgcgAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTaaaa                          >  2:647089/1‑148 (MQ=255)
                                                                                                                        gcgAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAAc                         >  1:1102912/1‑148 (MQ=255)
                                                                                                                            ggagCTGGATCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCt                     <  1:828948/148‑1 (MQ=255)
                                                                                                                             gtgCTGGATGAGGAAGCTAGTGATCGAAAGAGTCATAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTTAAACCCCTa                    <  1:2662995/146‑1 (MQ=255)
                                                                                                                             gagCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTa                    >  2:2631997/1‑148 (MQ=255)
                                                                                                                              agCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTAt                   <  2:801949/148‑1 (MQ=255)
                                                                                                                               gCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTAtt                  <  1:2125708/148‑1 (MQ=255)
                                                                                                                                cTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTc                 >  2:2720756/1‑148 (MQ=255)
                                                                                                                                 tGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTcc                <  2:1018966/148‑1 (MQ=255)
                                                                                                                                  ggATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCa               <  1:1171857/148‑1 (MQ=255)
                                                                                                                                   gATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCAt              >  2:1508569/1‑148 (MQ=255)
                                                                                                                                     tCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTg            >  2:1136916/1‑148 (MQ=255)
                                                                                                                                     tCAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTg            <  1:3150533/148‑1 (MQ=255)
                                                                                                                                      cAGGAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTgg           >  1:980384/1‑148 (MQ=255)
                                                                                                                                       aGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAACACCCCTATTCCATTggg          <  1:2609960/148‑1 (MQ=255)
                                                                                                                                        ggAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGc         >  2:2742305/1‑148 (MQ=255)
                                                                                                                                        ggAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGc         >  1:595702/1‑148 (MQ=255)
                                                                                                                                        ggAAGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGc         >  1:3060240/1‑148 (MQ=255)
                                                                                                                                         gAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCg        <  1:2902910/148‑1 (MQ=255)
                                                                                                                                         gAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCg        >  1:2484575/1‑148 (MQ=255)
                                                                                                                                          aaGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGt       >  2:783412/1‑148 (MQ=255)
                                                                                                                                           aGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTg      >  2:2416308/1‑148 (MQ=255)
                                                                                                                                           aGCTAGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTg      <  1:1991715/148‑1 (MQ=255)
                                                                                                                                             cTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGAc    >  1:2647219/1‑148 (MQ=255)
                                                                                                                                              tCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACt   >  1:2156534/1‑148 (MQ=255)
                                                                                                                                               cGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTg  >  1:3218943/1‑148 (MQ=255)
                                                                                                                                               cGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCAATTGGGCGTGACTg  >  1:627343/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGAGGAGCTGGATCAGGAAGCTCGTGATCGCAAGCGTCAGAAAAAACGTCGTGGTCATGCGCCGGGCAGCCGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTG  >  CP009273/4044648‑4044938

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: