Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,941,341 G→A 100% L343L (CTC→CTT amiC ← N‑acetylmuramoyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,941,3410GA96.4% 792.1 / 13.1 250L343L (CTC→CTTamiCN‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base A (135/106);  total (140/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 250 total aligned reads displayed.

GGAAAGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTTGTGC  >  CP009273/2941196‑2941484
                                                                                                                                                 |                                                                                                                                               
ggAAAGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGGTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAg                                                                                                                                               <  1:3184576/148‑1 (MQ=255)
  aaaGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCt                                                                                                                                             >  2:2377901/1‑148 (MQ=255)
    aGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTa                                                                                                                                           <  2:233768/148‑1 (MQ=255)
        gCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCa                                                                                                                                       >  2:2119665/1‑148 (MQ=255)
        gCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCa                                                                                                                                       <  2:3542792/148‑1 (MQ=255)
         ccGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGACAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAg                                                                                                                                      <  1:3483594/148‑1 (MQ=255)
           gTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTa                                                                                                                                    <  1:2487329/148‑1 (MQ=255)
           gTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTa                                                                                                                                    <  1:39904/148‑1 (MQ=255)
            ttttCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGACTATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTAc                                                                                                                                   >  1:675468/1‑148 (MQ=255)
             tttCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTAcc                                                                                                                                  >  1:1884202/1‑148 (MQ=255)
              ttCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCg                                                                                                                                 <  1:2229203/148‑1 (MQ=255)
                cAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCt                                                                                                                               <  1:1652558/148‑1 (MQ=255)
                   tttACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTa                                                                                                                            >  2:1847106/1‑148 (MQ=255)
                   tttACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTa                                                                                                                            >  1:2749187/1‑148 (MQ=255)
                   tttACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTa                                                                                                                            >  1:2519707/1‑148 (MQ=255)
                   tttACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTa                                                                                                                            >  1:1314480/1‑148 (MQ=255)
                       cGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCaa                                                                                                                        <  1:601795/148‑1 (MQ=255)
                          tctTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACt                                                                                                                     >  1:1830119/1‑148 (MQ=255)
                          tctTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACt                                                                                                                     >  1:1367745/1‑148 (MQ=255)
                             tCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCa                                                                                                                  <  1:1173331/148‑1 (MQ=255)
                             tCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCa                                                                                                                  <  1:1721532/148‑1 (MQ=255)
                                tCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGc                                                                                                               >  2:3476527/1‑148 (MQ=255)
                                 cAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCt                                                                                                              >  2:3134479/1‑148 (MQ=255)
                                 cAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCt                                                                                                              >  2:2853247/1‑148 (MQ=255)
                                   aCGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGt                                                                                                            >  1:3489266/1‑148 (MQ=255)
                                   aCGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGt                                                                                                            >  1:1771680/1‑148 (MQ=255)
                                    cGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTc                                                                                                           >  2:3376721/1‑148 (MQ=255)
                                       tACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGc                                                                                                        >  2:1725893/1‑148 (MQ=255)
                                         cTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCaa                                                                                                      >  2:378816/1‑148 (MQ=255)
                                          ttaTAAACTCCGTTTCGACCTGAATGGAGGGAATATCTGGTGCCTTTAGTACTTCAAACCCTGCATGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAAt                                                                                                     <  2:1654718/146‑1 (MQ=255)
                                          tGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAAt                                                                                                     >  1:3105655/1‑148 (MQ=255)
                                           gATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATg                                                                                                    >  2:2856606/1‑148 (MQ=255)
                                             tAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGt                                                                                                  <  1:2811319/148‑1 (MQ=255)
                                               aaCGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCa                                                                                                >  2:2507519/1‑148 (MQ=255)
                                                gcGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAg                                                                                               >  2:839589/2‑148 (MQ=255)
                                                   ccGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGa                                                                                            <  2:2760321/148‑1 (MQ=255)
                                                           aCCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCa                                                                                    >  1:751202/1‑148 (MQ=255)
                                                            caaGATTGGAGGTAATATCTGGTGCCTTTATTACGGCAAACCCGGCCTTTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCat                                                                                   <  1:2913440/146‑1 (MQ=255)
                                                              aGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCatat                                                                                 >  1:1343797/1‑148 (MQ=255)
                                                              aGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCatat                                                                                 <  2:1587453/148‑1 (MQ=255)
                                                              aGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCatat                                                                                 <  2:2897198/148‑1 (MQ=255)
                                                                 aTGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGa                                                                              <  1:2651379/148‑1 (MQ=255)
                                                                  tGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGaa                                                                             >  1:1940859/1‑148 (MQ=255)
                                                                   ggAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAAc                                                                            >  1:458992/1‑148 (MQ=255)
                                                                    gAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACa                                                                           >  2:1798609/1‑148 (MQ=255)
                                                                       ggAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATgg                                                                        <  1:839583/148‑1 (MQ=255)
                                                                         aaTATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGgtg                                                                      >  2:1089530/1‑148 (MQ=255)
                                                                            atCTGTTGCCTTTTGTACGGCAAACCCGGCCTTTTCAACTTTATTTTTATGCAGGTTTTTGATTTTACCAAGGTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGAATTGGTCAGCGACTGTACCATATCGAACATGGTGTgg                                                                   <  1:786631/148‑1 (MQ=255)
                                                                            atCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTgg                                                                   >  2:2083619/1‑148 (MQ=255)
                                                                            atCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTgg                                                                   >  1:888537/1‑148 (MQ=255)
                                                                            atCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTgg                                                                   >  1:766057/1‑148 (MQ=255)
                                                                                 gTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTAACAAACTTCAGGCTGTCGGCAATGGGCAGCGACTGTACCATATCGAACATCGTGTGGTCGAc                                                              >  1:2423995/1‑148 (MQ=255)
                                                                                 gTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACATGATTTATATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAAAGGTCAGCGACTGTACCATATCGAACATGGTTTGGTCGAc                                                              >  1:303676/1‑148 (MQ=255)
                                                                                   gCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACAt                                                            >  1:2114642/1‑148 (MQ=255)
                                                                                    ccTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATa                                                           >  2:869412/1‑148 (MQ=255)
                                                                                    ccTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATa                                                           >  2:1982091/1‑148 (MQ=255)
                                                                                          gTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTc                                                     >  1:1466598/1‑148 (MQ=255)
                                                                                           tACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCa                                                    >  2:1400758/1‑148 (MQ=255)
                                                                                           tACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCa                                                    >  1:3040839/1‑148 (MQ=255)
                                                                                                  aaCCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCttt                                             >  2:1019676/1‑148 (MQ=255)
                                                                                                   aCCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCtttt                                            <  2:1314533/148‑1 (MQ=255)
                                                                                                   aCCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCtttt                                            >  2:497086/1‑148 (MQ=255)
                                                                                                   aCCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCtttt                                            >  1:2344598/1‑148 (MQ=255)
                                                                                                    cccGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTg                                           >  1:223451/1‑148 (MQ=255)
                                                                                                               cAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGccacc                                >  2:2983288/1‑148 (MQ=255)
                                                                                                               cAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTATACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACAGTACCATATCGAACATGGAGTGGTCGACATAGCGGTCACCGCTTATGCTCACGccacc                                >  2:1204215/1‑148 (MQ=255)
                                                                                                                aaCTTTTTTTTTATGCTGCTTGTTGTTTTTACCACGATTATTCAGTACCGATTTACCAAAGTTCGGGCTGTCGGAAATGGTCCGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGccacca                               <  1:3278908/148‑1 (MQ=255)
                                                                                                                aaCTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGccacca                               <  1:2377901/148‑1 (MQ=255)
                                                                                                                 aCTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTACAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCaa                              <  1:716570/148‑1 (MQ=255)
                                                                                                                  cTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAAt                             >  1:2184921/1‑148 (MQ=255)
                                                                                                                   ttGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATc                            <  1:879352/148‑1 (MQ=255)
                                                                                                                    tGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCa                           >  1:1330860/1‑148 (MQ=255)
                                                                                                                      aTTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAg                         <  1:3435582/148‑1 (MQ=255)
                                                                                                                      aTTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAg                         >  2:3219091/1‑148 (MQ=255)
                                                                                                                       tttttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGt                        <  2:1771680/148‑1 (MQ=255)
                                                                                                                       tttttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGt                        >  1:573239/1‑148 (MQ=255)
                                                                                                                       tttttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGt                        <  2:3015369/148‑1 (MQ=255)
                                                                                                                        ttttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTc                       <  1:1188267/148‑1 (MQ=255)
                                                                                                                        ttttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTc                       <  2:3076448/148‑1 (MQ=255)
                                                                                                                          ttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCg                     <  1:557622/148‑1 (MQ=255)
                                                                                                                           tATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGa                    <  2:2749187/148‑1 (MQ=255)
                                                                                                                             tGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAgg                  <  1:553139/148‑1 (MQ=255)
                                                                                                                              gCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGc                 <  1:2948566/148‑1 (MQ=255)
                                                                                                                              gCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGc                 <  2:3230695/148‑1 (MQ=255)
                                                                                                                                 gCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGtt              <  1:1847106/148‑1 (MQ=255)
                                                                                                                                 gCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGtt              <  2:676135/148‑1 (MQ=255)
                                                                                                                                 gCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGtt              <  2:3451064/148‑1 (MQ=255)
                                                                                                                                    tgttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTg           <  1:1759928/148‑1 (MQ=255)
                                                                                                                                    tgttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTg           <  2:2552428/148‑1 (MQ=255)
                                                                                                                                     gttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTgg          <  2:3290335/148‑1 (MQ=255)
                                                                                                                                      ttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTggg         <  1:259826/148‑1 (MQ=255)
                                                                                                                                      ttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTggg         <  1:2062616/148‑1 (MQ=255)
                                                                                                                                         aTTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGttt      <  2:2253273/148‑1 (MQ=255)
                                                                                                                                          ttttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtg     <  2:3040839/148‑1 (MQ=255)
                                                                                                                                          ttttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtg     <  1:2119642/148‑1 (MQ=255)
                                                                                                                                          ttttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtg     >  2:1541906/1‑148 (MQ=255)
                                                                                                                                           tttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtgt    >  1:1626921/1‑148 (MQ=255)
                                                                                                                                           tttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtgt    >  1:1463194/1‑148 (MQ=255)
                                                                                                                                            ttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtgtg   >  1:1075143/1‑148 (MQ=255)
                                                                                                                                            ttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtgtg   >  1:1074886/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GGAAAGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTTGTGC  >  CP009273/2941196‑2941484

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: