Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,218,022 C→A 100% 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,0220CA99.0% 975.2 / ‑1.3 296E110* (GAA→TAA) pbpGD‑alanyl‑D‑alanine endopeptidase
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base A (123/169);  total (127/169)
Fisher's exact test for biased strand distribution p-value = 3.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 296 total aligned reads displayed.

TGTTGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCA  >  CP009273/2217875‑2218168
                                                                                                                                                   |                                                                                                                                                  
  ttgttCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTa                                                                                                                                                >  1:3648273/1‑150 (MQ=255)
   tgttCATTCCGAGCGATGTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTaa                                                                                                                                               <  2:3062648/150‑2 (MQ=255)
       cATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCAAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAg                                                                                                                                           <  1:349693/150‑1 (MQ=255)
        aTTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGt                                                                                                                                          >  1:4239706/1‑150 (MQ=255)
          tCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCg                                                                                                                                        >  1:3615524/1‑150 (MQ=255)
            cGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTa                                                                                                                                      >  2:2501729/1‑150 (MQ=255)
             gAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAc                                                                                                                                     >  2:4677102/1‑150 (MQ=255)
                 gATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTAcgcgc                                                                                                                                 >  2:1826369/1‑150 (MQ=255)
                    ttcgctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATCTCTTTACGGCTGAGTTCGCTCTTCCGTCGTACGCGCGac                                                                                                                              >  2:4094552/1‑149 (MQ=255)
                       gctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAg                                                                                                                           <  2:1448091/150‑1 (MQ=255)
                        ctttcgcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGa                                                                                                                          >  1:3631580/1‑150 (MQ=255)
                             gcATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCt                                                                                                                     >  2:2270049/1‑150 (MQ=255)
                              cattcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCtt                                                                                                                    >  2:2310817/1‑150 (MQ=255)
                               attcattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCttt                                                                                                                   <  2:3598253/150‑1 (MQ=255)
                                  cattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAt                                                                                                                >  1:136919/1‑150 (MQ=255)
                                  cattGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAt                                                                                                                <  1:2828824/150‑1 (MQ=255)
                                   attGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtc                                                                                                               <  2:4545393/150‑1 (MQ=255)
                                   attGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCAtc                                                                                                               <  1:1865590/150‑1 (MQ=255)
                                      gCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTcc                                                                                                            <  2:3654254/150‑1 (MQ=255)
                                       ccTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCg                                                                                                           <  2:1222111/150‑1 (MQ=255)
                                       ccTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCg                                                                                                           <  1:3364148/150‑1 (MQ=255)
                                        cTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCgg                                                                                                          <  2:4532006/150‑1 (MQ=255)
                                           aataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGt                                                                                                       >  1:1715627/1‑150 (MQ=255)
                                            ataaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTc                                                                                                      <  2:3435584/150‑1 (MQ=255)
                                             taaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCt                                                                                                     <  1:4649913/150‑1 (MQ=255)
                                             taaAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCt                                                                                                     <  2:1575749/150‑1 (MQ=255)
                                              aaaGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTg                                                                                                    <  2:3341993/150‑1 (MQ=255)
                                              aaaGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTg                                                                                                    <  2:2233699/150‑1 (MQ=255)
                                               aaGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTgg                                                                                                   >  1:3750277/1‑150 (MQ=255)
                                                aGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGc                                                                                                  <  2:4039599/150‑1 (MQ=255)
                                                   cTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGa                                                                                               >  1:3250260/1‑150 (MQ=255)
                                                    tttGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGat                                                                                              <  1:2155385/150‑1 (MQ=255)
                                                      tGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGatat                                                                                            >  2:3909071/1‑150 (MQ=255)
                                                      tGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGatat                                                                                            >  1:585838/1‑150 (MQ=255)
                                                      tGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGatac                                                                                            <  1:1676895/150‑2 (MQ=255)
                                                         aaCCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCa                                                                                         <  1:1564237/150‑1 (MQ=255)
                                                             accGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTCTTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACttt                                                                                     >  2:1140447/1‑150 (MQ=255)
                                                               cGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTa                                                                                   <  2:84527/150‑1 (MQ=255)
                                                                  gATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGtt                                                                                <  2:3223017/150‑1 (MQ=255)
                                                                  gATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGtt                                                                                <  2:2863859/150‑1 (MQ=255)
                                                                    tAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGtttt                                                                              <  1:29538/150‑1 (MQ=255)
                                                                     aaTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttttt                                                                             >  2:3171158/1‑150 (MQ=255)
                                                                     aaTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGttttt                                                                             <  2:547739/150‑1 (MQ=255)
                                                                      aTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTc                                                                            <  1:1015332/150‑1 (MQ=255)
                                                                      aTGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTc                                                                            <  2:1635849/150‑1 (MQ=255)
                                                                        ggTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCAt                                                                          >  1:1720126/1‑150 (MQ=255)
                                                                          tGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACATAATATCTTTACGGCTGATTTAGATATTCAGTCGTACGCTCTTATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATcc                                                                        <  1:2904332/150‑1 (MQ=255)
                                                                          tGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATcc                                                                        >  1:2175922/1‑150 (MQ=255)
                                                                           gcgcAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACttt                                                                                     <  2:4784972/136‑1 (MQ=255)
                                                                           gcgcAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACttt                                                                                     >  1:4784972/1‑136 (MQ=255)
                                                                                 ggCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCa                                                                 >  2:2604745/1‑150 (MQ=255)
                                                                                 ggCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCa                                                                 >  2:2426823/1‑150 (MQ=255)
                                                                                 ggCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCa                                                                 <  1:2117090/150‑1 (MQ=255)
                                                                                    ttGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTc                                                              <  1:3418209/150‑1 (MQ=255)
                                                                                    ttGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTc                                                              <  1:4540791/150‑1 (MQ=255)
                                                                                    ttGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTc                                                              <  1:2610031/150‑1 (MQ=255)
                                                                                    ttGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTc                                                              >  2:2954431/1‑150 (MQ=255)
                                                                                       ccgccgCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTg                                                           <  2:162642/150‑1 (MQ=255)
                                                                                        cgccgcGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGc                                                          <  2:2807367/150‑1 (MQ=255)
                                                                                           cgcgcgGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATc                                                       >  1:3936482/1‑150 (MQ=255)
                                                                                           cgcgcgGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATc                                                       <  2:2592518/150‑1 (MQ=255)
                                                                                            gcgcgGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATCTCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATcc                                                      >  2:4329035/1‑150 (MQ=255)
                                                                                             cgcgGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCa                                                     >  2:753352/1‑150 (MQ=255)
                                                                                               cgGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGc                                                   >  1:2651623/1‑150 (MQ=255)
                                                                                                   tttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGACTAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACaa                                               <  2:1186271/150‑1 (MQ=255)
                                                                                                   tttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACaa                                               >  2:4690701/1‑150 (MQ=255)
                                                                                                   tttCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACaa                                               >  1:868865/1‑150 (MQ=255)
                                                                                                     tCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAAcc                                             <  1:1484585/150‑1 (MQ=255)
                                                                                                      cTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCa                                            >  2:3976569/1‑150 (MQ=255)
                                                                                                             cATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGt                                     <  1:4960543/150‑1 (MQ=255)
                                                                                                              aTCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTc                                    >  1:1285248/1‑150 (MQ=255)
                                                                                                                cAGCGCCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCAt                                  <  1:116768/150‑1 (MQ=255)
                                                                                                                 agagcCAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGAGTCTGGCTGATATCCACTTTTATTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCAGCGCGGTCAtt                                 >  1:526977/4‑150 (MQ=255)
                                                                                                                     ccAGCAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAAtt                             >  2:923587/1‑150 (MQ=255)
                                                                                                                         cAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCt                         >  1:5122139/1‑150 (MQ=255)
                                                                                                                         cAACAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCt                         >  1:3187026/1‑150 (MQ=255)
                                                                                                                          aaCAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTg                        >  2:2010585/1‑150 (MQ=255)
                                                                                                                           aCAGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGa                       <  2:3697047/150‑1 (MQ=255)
                                                                                                                             aGCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATa                     <  2:616803/150‑1 (MQ=255)
                                                                                                                              gCATATCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAg                    >  2:2657252/1‑150 (MQ=255)
                                                                                                                                 tatCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACg                 <  1:744901/150‑1 (MQ=255)
                                                                                                                                 tatCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACg                 >  2:4918904/1‑150 (MQ=255)
                                                                                                                                 tatCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACg                 <  1:3029171/150‑1 (MQ=255)
                                                                                                                                  atCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGc                >  2:3889185/1‑150 (MQ=255)
                                                                                                                                   tCTTTACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCa               <  2:4247046/150‑1 (MQ=255)
                                                                                                                                     tttACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAAt             <  2:1029576/150‑1 (MQ=255)
                                                                                                                                       tACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCg           <  2:1527496/150‑1 (MQ=255)
                                                                                                                                       tACGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCg           <  2:2167240/150‑1 (MQ=255)
                                                                                                                                        aCGGCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCgg          <  2:1484277/150‑1 (MQ=255)
                                                                                                                                          ggCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGAc        >  1:1446030/1‑150 (MQ=255)
                                                                                                                                           gCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGAca       <  2:597126/150‑2 (MQ=255)
                                                                                                                                           gCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACg       <  1:2501729/150‑1 (MQ=255)
                                                                                                                                           gCTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACg       >  2:1335196/1‑150 (MQ=255)
                                                                                                                                            cTGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGc      <  2:3513971/150‑1 (MQ=255)
                                                                                                                                             tGATTTAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCa     <  1:4287343/150‑1 (MQ=255)
                                                                                                                                                tttAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCa  <  1:4160914/150‑1 (MQ=255)
                                                                                                                                                tttAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCa  <  1:5001917/150‑1 (MQ=255)
                                                                                                                                                tttAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCa  <  2:4428317/150‑1 (MQ=255)
                                                                                                                                                tttAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCa  <  1:1252165/150‑1 (MQ=255)
                                                                                                                                                tttAGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCa  <  2:3750277/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGTTGTTCATTCCGAGCGATTTCGCTTTCGCATTCATTGCCTTAATAAAGGCTTTGTAACCACCGGGATAATGGTGCGCAAGGCTTGCCGCCGCGCGGTTTTCTGAAGACATCAGCGCCAGCAACAGCATATCTTTACGGCTGATTTCGCTATTCAGTCGTACGCGCGAATAGACCCCTTTCATCTCCGGCGTCTGGCTGATATCCACTTTTAGTTTTTCATCCAGCGGCAGTCGTGCATCCAGCACAACCATCGCGGTCATTAATTTGCTGATAGACGCAATCGGACGCACCA  >  CP009273/2217875‑2218168

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: