Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,218,022 C→A 89.9% E110* (GAA→TAA)  pbpG ← D‑alanyl‑D‑alanine endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,218,0220CA89.9% 462.9 / 45.3 179E110* (GAA→TAA) pbpGD‑alanyl‑D‑alanine endopeptidase
Reads supporting (aligned to +/- strand):  ref base C (11/7);  new base A (90/71);  total (101/78)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 179 total aligned reads displayed.

GTTGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCA  >  CP009273/2217876‑2218165
                                                                                                                                                  |                                                                                                                                               
 ttgttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGc                                                                                                                                               >  2:3181089/1‑148 (MQ=255)
   gttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTa                                                                                                                                             >  1:1021709/1‑148 (MQ=255)
   gttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTa                                                                                                                                             >  1:1021671/1‑148 (MQ=255)
       aTTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCa                                                                                                                                         >  1:2324416/1‑148 (MQ=255)
       aTTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCa                                                                                                                                         >  1:707275/1‑148 (MQ=255)
        ttCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAg                                                                                                                                        <  2:627096/148‑1 (MQ=255)
         tCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGt                                                                                                                                       >  2:2586357/1‑148 (MQ=255)
             aGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTa                                                                                                                                   <  1:149184/148‑1 (MQ=255)
             aGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTa                                                                                                                                   >  1:1836132/1‑148 (MQ=255)
             aGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTa                                                                                                                                   <  1:2537894/148‑1 (MQ=255)
             aGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTa                                                                                                                                   <  1:1661349/148‑1 (MQ=255)
              gCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAc                                                                                                                                  >  2:3190314/1‑148 (MQ=255)
                  tttcgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTCTTCAGTCGTAcgcgc                                                                                                                              >  2:2547062/1‑148 (MQ=255)
                  tttcgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAcgcgc                                                                                                                              >  2:434896/1‑148 (MQ=255)
                  tttcgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAcgcgc                                                                                                                              >  2:605151/1‑148 (MQ=255)
                  tttcgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAcgcgc                                                                                                                              >  1:2193163/1‑148 (MQ=255)
                    tcgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGa                                                                                                                            <  1:1342623/148‑1 (MQ=255)
                    tcgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGa                                                                                                                            >  2:2161319/1‑148 (MQ=255)
                     cgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCAGATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGaa                                                                                                                           >  1:2582095/1‑148 (MQ=255)
                      gctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAAt                                                                                                                          >  1:1717733/1‑148 (MQ=255)
                        tttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAg                                                                                                                        <  2:2070363/148‑1 (MQ=255)
                         ttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCCGCACCAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGa                                                                                                                       >  2:3003907/1‑148 (MQ=255)
                         ttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGa                                                                                                                       <  1:1585315/148‑1 (MQ=255)
                         ttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGa                                                                                                                       >  1:831830/1‑148 (MQ=255)
                         ttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGa                                                                                                                       <  1:3224216/148‑1 (MQ=255)
                            gcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGAccc                                                                                                                    <  2:1536317/148‑1 (MQ=255)
                             cattcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGAcccc                                                                                                                   >  2:533744/1‑148 (MQ=255)
                             cattcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGAcccc                                                                                                                   >  1:2961253/1‑148 (MQ=255)
                              attcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCt                                                                                                                  <  2:1046683/148‑1 (MQ=255)
                              attcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCt                                                                                                                  <  1:388933/148‑1 (MQ=255)
                              attcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCt                                                                                                                  >  1:771824/1‑148 (MQ=255)
                                tcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCttt                                                                                                                >  2:1989377/1‑148 (MQ=255)
                                tcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCttt                                                                                                                >  2:2629074/1‑148 (MQ=255)
                                tcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCttt                                                                                                                <  1:2436728/148‑1 (MQ=255)
                                 cattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTc                                                                                                               <  1:166760/148‑1 (MQ=255)
                                 cattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTc                                                                                                               <  1:520491/148‑1 (MQ=255)
                                   ttGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAt                                                                                                             >  2:2143585/1‑148 (MQ=255)
                                    tGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtc                                                                                                            >  1:2854767/1‑148 (MQ=255)
                                    tGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtc                                                                                                            >  1:2550132/1‑148 (MQ=255)
                                    tGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtc                                                                                                            >  1:2576018/1‑148 (MQ=255)
                                     gCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtct                                                                                                           <  1:605151/148‑1 (MQ=255)
                                      ccTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACAACAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTAACAtctc                                                                                                          >  1:867714/1‑148 (MQ=255)
                                        ttAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCg                                                                                                        >  1:2150528/1‑148 (MQ=255)
                                                ggCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGAGTCTgg                                                                                                >  2:2070387/1‑148 (MQ=255)
                                                  cTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCt                                                                                              >  2:3251159/1‑148 (MQ=255)
                                                  cTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCt                                                                                              >  1:1290326/1‑148 (MQ=255)
                                                  cTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCt                                                                                              >  2:107975/1‑148 (MQ=255)
                                                    ttGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGa                                                                                            >  1:1156363/1‑148 (MQ=255)
                                                    ttGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGa                                                                                            >  1:2045862/1‑148 (MQ=255)
                                                    ttGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGa                                                                                            <  1:2682317/148‑1 (MQ=255)
                                                    ttGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGa                                                                                            >  1:986067/1‑148 (MQ=255)
                                                        aaCCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATc                                                                                        <  1:3157123/148‑1 (MQ=255)
                                                         accaccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATcc                                                                                       <  1:541765/148‑1 (MQ=255)
                                                         accaccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGAGATcc                                                                                       >  1:2760303/1‑148 (MQ=255)
                                                          ccaccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGACTAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCa                                                                                      >  2:148304/1‑148 (MQ=255)
                                                          ccaccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCa                                                                                      <  2:2693050/148‑1 (MQ=255)
                                                          ccaccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCa                                                                                      <  2:594864/148‑1 (MQ=255)
                                                                ggATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTa                                                                                <  2:546167/148‑1 (MQ=255)
                                                                  aTAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGt                                                                              <  2:1616998/148‑1 (MQ=255)
                                                                   tAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGtt                                                                             <  2:3017823/148‑1 (MQ=255)
                                                                    aaTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttt                                                                            <  2:3044374/148‑1 (MQ=255)
                                                                    aaTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttt                                                                            <  2:301223/148‑1 (MQ=255)
                                                                     aTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTCCGCGCGAATAGACCCCTATCATCTCCGGCGTCTGGCTGATATCCACTTTTAGtttt                                                                           >  2:2849028/1‑148 (MQ=255)
                                                                      tGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttttt                                                                          >  2:2907336/1‑148 (MQ=255)
                                                                             cAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAg                                                                   <  1:2194341/148‑1 (MQ=255)
                                                                                ggCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCgg                                                                >  2:2001994/1‑148 (MQ=255)
                                                                                ggCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCgg                                                                >  2:1559866/1‑148 (MQ=255)
                                                                                      ccgccgCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCg                                                          <  2:2371879/148‑1 (MQ=255)
                                                                                              cgGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCa                                                  >  2:957821/1‑148 (MQ=255)
                                                                                              cgGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCa                                                  >  2:1834859/1‑148 (MQ=255)
                                                                                                   ttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGcaca                                             <  2:3230525/148‑1 (MQ=255)
                                                                                                          gACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATcg                                      >  2:1886539/1‑148 (MQ=255)
                                                                                                           aCATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATcgc                                     >  2:1185532/1‑148 (MQ=255)
                                                                                                             aTCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCgg                                   >  1:3021908/1‑148 (MQ=255)
                                                                                                              tCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGt                                  <  2:864900/148‑1 (MQ=255)
                                                                                                               cAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTc                                 >  2:2435214/1‑148 (MQ=255)
                                                                                                               cAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTc                                 >  2:776125/1‑148 (MQ=255)
                                                                                                               cAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTc                                 >  2:608457/1‑148 (MQ=255)
                                                                                                               cAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTc                                 >  1:2049062/1‑148 (MQ=255)
                                                                                                                  cgcCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCATCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCAtt                              >  2:1734946/1‑148 (MQ=255)
                                                                                                                          aCAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCt                      >  1:158057/1‑148 (MQ=255)
                                                                                                                          aCAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCt                      >  1:2321890/1‑148 (MQ=255)
                                                                                                                           cAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTg                     >  2:3255971/1‑148 (MQ=255)
                                                                                                                            aGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGa                    <  2:2468334/148‑1 (MQ=255)
                                                                                                                            aGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGa                    >  1:1307041/1‑148 (MQ=255)
                                                                                                                              cATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATACACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATa                  >  2:3166723/1‑148 (MQ=255)
                                                                                                                                tatCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGa                <  2:697460/148‑1 (MQ=255)
                                                                                                                                tatCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGa                >  2:840065/1‑148 (MQ=255)
                                                                                                                                tatCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGa                >  2:842420/1‑148 (MQ=255)
                                                                                                                                tatCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGa                >  2:842435/1‑148 (MQ=255)
                                                                                                                                tatCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGa                >  2:2722092/1‑148 (MQ=255)
                                                                                                                                  tCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACg              >  1:2435887/1‑148 (MQ=255)
                                                                                                                                  tCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACg              >  2:1444305/1‑148 (MQ=255)
                                                                                                                                     ttACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCaa           <  2:882564/148‑1 (MQ=255)
                                                                                                                                      tACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAAt          <  1:2685576/148‑1 (MQ=255)
                                                                                                                                      tACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAAt          <  1:2685596/148‑1 (MQ=255)
                                                                                                                                         ggCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCgg       <  1:755937/148‑1 (MQ=255)
                                                                                                                                             gATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGc   <  2:1836132/148‑1 (MQ=255)
                                                                                                                                              aTTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCa  >  1:2723166/1‑148 (MQ=255)
                                                                                                                                              aTTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCa  >  1:1069944/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GTTGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCA  >  CP009273/2217876‑2218165

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