Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,218,022 C→A 94.6% 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,0220CA94.6% 600.3 / 23.9 204E110* (GAA→TAA) pbpGD‑alanyl‑D‑alanine endopeptidase
Reads supporting (aligned to +/- strand):  ref base C (6/5);  new base A (81/112);  total (87/117)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 204 total aligned reads displayed.

TGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCAC  >  CP009273/2217878‑2218166
                                                                                                                                                |                                                                                                                                                
tgttCAGTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCt                                                                                                                                               <  2:2413868/148‑1 (MQ=255)
 gttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTa                                                                                                                                              <  1:727576/148‑1 (MQ=255)
  ttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTAt                                                                                                                                             >  1:2620794/1‑148 (MQ=255)
   tCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTAtt                                                                                                                                            <  1:1054262/148‑1 (MQ=255)
   tCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTAtt                                                                                                                                            <  2:2780641/148‑1 (MQ=255)
   tCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTAtt                                                                                                                                            <  1:3377305/148‑1 (MQ=255)
        ccGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTc                                                                                                                                       >  1:3553452/1‑148 (MQ=255)
          gAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTGACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGt                                                                                                                                     <  1:1965913/148‑1 (MQ=255)
          gAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGt                                                                                                                                     <  2:864537/148‑1 (MQ=255)
          gAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGt                                                                                                                                     <  2:1136988/148‑1 (MQ=255)
          gAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGt                                                                                                                                     <  1:1483033/148‑1 (MQ=255)
           aGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTa                                                                                                                                    <  2:163710/148‑1 (MQ=255)
           aGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTa                                                                                                                                    <  2:1868880/148‑1 (MQ=255)
             cGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAcg                                                                                                                                  >  1:3045265/1‑148 (MQ=255)
               aTTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAcgcg                                                                                                                                <  1:703334/148‑1 (MQ=255)
               aTTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAcgcg                                                                                                                                >  1:795364/1‑148 (MQ=255)
                 ttcgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAcgcgcg                                                                                                                              <  1:2499220/148‑1 (MQ=255)
                  tcgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGa                                                                                                                             <  1:658467/148‑1 (MQ=255)
                        tcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGAc                                                                                                                       <  1:1672080/148‑1 (MQ=255)
                        tcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGAc                                                                                                                       >  2:1958310/1‑148 (MQ=255)
                            attcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCt                                                                                                                   <  2:582659/148‑1 (MQ=255)
                              tcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCttt                                                                                                                 >  2:1776369/1‑148 (MQ=255)
                               cattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTc                                                                                                                <  1:2210050/148‑1 (MQ=255)
                                  tGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtc                                                                                                             <  2:3195898/148‑1 (MQ=255)
                                   gCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtct                                                                                                            <  2:1930400/148‑1 (MQ=255)
                                   gCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtct                                                                                                            <  2:1930485/148‑1 (MQ=255)
                                   gCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtct                                                                                                            <  2:889953/148‑1 (MQ=255)
                                      ttAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCg                                                                                                         <  2:1324934/148‑1 (MQ=255)
                                      ttAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCg                                                                                                         >  1:2941174/1‑148 (MQ=255)
                                        aataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGc                                                                                                       <  2:1071734/148‑1 (MQ=255)
                                         ataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCg                                                                                                      <  2:3419275/148‑1 (MQ=255)
                                         ataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCg                                                                                                      >  1:266934/1‑148 (MQ=255)
                                         ataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCg                                                                                                      >  2:1440381/1‑148 (MQ=255)
                                          ttaaGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCAATAGACCCCTTTCATCTCCGGCGTc                                                                                                    <  1:280312/146‑1 (MQ=255)
                                          taaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGt                                                                                                     >  1:906722/1‑148 (MQ=255)
                                           aaaGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTc                                                                                                    <  1:913582/148‑1 (MQ=255)
                                           aaaGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTc                                                                                                    >  1:1288509/1‑148 (MQ=255)
                                             aGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTg                                                                                                  <  2:2030142/148‑1 (MQ=255)
                                             aGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTg                                                                                                  <  2:2512665/148‑1 (MQ=255)
                                               gCTTTGTAACCACCGTGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGc                                                                                                <  1:452343/148‑1 (MQ=255)
                                               gCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGc                                                                                                <  2:2395644/148‑1 (MQ=255)
                                                 tttGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTg                                                                                              >  1:2688824/1‑148 (MQ=255)
                                                         caccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCAc                                                                                      <  2:1724372/148‑1 (MQ=255)
                                                         caccGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCAc                                                                                      <  1:652909/148‑1 (MQ=255)
                                                              ggATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTa                                                                                 >  1:2590945/1‑148 (MQ=255)
                                                              ggATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTa                                                                                 <  2:2236432/148‑1 (MQ=255)
                                                                aTAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGt                                                                               <  1:3433568/148‑1 (MQ=255)
                                                                aTAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGt                                                                               <  1:2999543/148‑1 (MQ=255)
                                                                    tGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttttt                                                                           <  1:443121/148‑1 (MQ=255)
                                                                    tGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttttt                                                                           >  1:389161/1‑148 (MQ=255)
                                                                    tGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttttt                                                                           <  1:3265337/148‑1 (MQ=255)
                                                                     ggTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTc                                                                          <  1:2385998/148‑1 (MQ=255)
                                                                             aGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCg                                                                  >  1:28375/1‑148 (MQ=255)
                                                                             aGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCg                                                                  >  1:28183/1‑148 (MQ=255)
                                                                                cTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCa                                                               <  1:2161070/148‑1 (MQ=255)
                                                                                cTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCa                                                               >  1:2526746/1‑148 (MQ=255)
                                                                                 ttGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAg                                                              <  1:306433/148‑1 (MQ=255)
                                                                                    ccgccgCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCg                                                           <  2:2790264/148‑1 (MQ=255)
                                                                                      gccgcGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTg                                                         <  2:8428/148‑1 (MQ=255)
                                                                                          cgcgGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATc                                                     <  2:2943815/148‑1 (MQ=255)
                                                                                                tttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGcac                                               >  2:2884398/1‑148 (MQ=255)
                                                                                                tttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGcac                                               >  1:179794/1‑148 (MQ=255)
                                                                                                 ttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGcaca                                              >  1:2732746/1‑148 (MQ=255)
                                                                                                 ttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCAATTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCAAGcaca                                              >  1:992673/1‑148 (MQ=255)
                                                                                                     gAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCa                                          >  2:2608034/1‑148 (MQ=255)
                                                                                                      aaGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCAt                                         <  2:634520/148‑1 (MQ=255)
                                                                                                      aaGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGATGATATCAACTTTTAGTTTTTCATCCATCGGAAGTCGTGCATCAAGCACAACCAt                                         >  1:2283236/1‑148 (MQ=255)
                                                                                                       aGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATc                                        >  2:3064676/1‑148 (MQ=255)
                                                                                                       aGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATc                                        >  2:1289430/1‑148 (MQ=255)
                                                                                                       aGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATc                                        >  1:548804/1‑148 (MQ=255)
                                                                                                       aGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATc                                        >  1:548881/1‑148 (MQ=255)
                                                                                                         aCATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATcgc                                      >  2:2433761/1‑148 (MQ=255)
                                                                                                          cATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATcgcg                                     <  2:3045076/148‑1 (MQ=255)
                                                                                                            tCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTCTTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGt                                   >  2:3303337/1‑148 (MQ=255)
                                                                                                            tCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGt                                   >  2:3558787/1‑148 (MQ=255)
                                                                                                            tCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGt                                   >  2:3002106/1‑148 (MQ=255)
                                                                                                             cAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTc                                  >  2:1413187/1‑148 (MQ=255)
                                                                                                              aGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCa                                 <  2:266934/148‑1 (MQ=255)
                                                                                                                cgcCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCAtt                               <  1:1158321/148‑1 (MQ=255)
                                                                                                                 gcCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTa                              <  2:1594295/148‑1 (MQ=255)
                                                                                                                  ccAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTaa                             <  2:884429/148‑1 (MQ=255)
                                                                                                                     gCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAAttt                          >  2:2260185/1‑148 (MQ=255)
                                                                                                                      cAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTg                         <  1:1376562/148‑1 (MQ=255)
                                                                                                                           gCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGAt                    >  1:1226574/1‑148 (MQ=255)
                                                                                                                           gCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGAt                    >  1:1903114/1‑148 (MQ=255)
                                                                                                                             atatCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAg                  <  1:1776369/148‑1 (MQ=255)
                                                                                                                                tCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACg               >  2:506710/1‑148 (MQ=255)
                                                                                                                                 cTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGc              >  2:3474434/1‑148 (MQ=255)
                                                                                                                                 cTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGc              <  1:1543905/148‑1 (MQ=255)
                                                                                                                                  tttACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCa             <  1:1921216/148‑1 (MQ=255)
                                                                                                                                   ttACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCaa            <  2:3249784/148‑1 (MQ=255)
                                                                                                                                    tACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAAt           <  2:1013308/148‑1 (MQ=255)
                                                                                                                                      cGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCg         >  1:1966960/1‑148 (MQ=255)
                                                                                                                                      cGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCg         <  2:3053773/148‑1 (MQ=255)
                                                                                                                                       ggCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCgg        <  1:891013/148‑1 (MQ=255)
                                                                                                                                        gCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGa       >  1:2004374/1‑148 (MQ=255)
                                                                                                                                        gCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGa       >  1:1684216/1‑148 (MQ=255)
                                                                                                                                         cTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGAc      >  2:3166188/1‑148 (MQ=255)
                                                                                                                                           gATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGc    <  1:3002107/148‑1 (MQ=255)
                                                                                                                                            aTTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCa   <  1:1953640/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCAC  >  CP009273/2217878‑2218166

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: