Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,406,553 C→A 100% R12L (CGG→CTG)  bsmA ← bioflm peroxide resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,406,5530CA97.0% 613.5 / 9.0 199R12L (CGG→CTG) bsmAbioflm peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (93/100);  total (96/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 199 total aligned reads displayed.

GCTTACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATG  >  CP009273/4406408‑4406698
                                                                                                                                                 |                                                                                                                                                 
gCTGACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAgg                                                                                                                                                 <  2:651620/148‑1 (MQ=255)
   tACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTaa                                                                                                                                              >  1:2919899/1‑148 (MQ=255)
      gCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATa                                                                                                                                           >  2:1396649/1‑148 (MQ=255)
       cTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATaa                                                                                                                                          <  2:418815/148‑1 (MQ=255)
        tACCTATTCGCTGTAACCCGTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCGGGCTGTGGCGTACCTTGCCGTGCGCTACAGGCACTTATCATCAACACCAATAATAAACTGGCAAACAGGtaactaac                                                                                                                                         >  2:3234105/1‑148 (MQ=255)
         aCCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACg                                                                                                                                        >  1:2391232/1‑148 (MQ=255)
         aCCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACg                                                                                                                                        >  2:1124450/1‑148 (MQ=255)
           cTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCt                                                                                                                                      <  1:1448160/148‑1 (MQ=255)
           cTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCt                                                                                                                                      >  1:1221822/1‑148 (MQ=255)
           cTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCt                                                                                                                                      <  2:844938/148‑1 (MQ=255)
            tATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTa                                                                                                                                     <  2:714534/148‑1 (MQ=255)
            tATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTa                                                                                                                                     <  2:1498211/148‑1 (MQ=255)
             aTTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTAt                                                                                                                                    >  2:2867619/1‑148 (MQ=255)
                cGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTAc                                                                                                                                 >  2:658331/1‑148 (MQ=255)
                   tGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGtt                                                                                                                              >  1:1828231/1‑148 (MQ=255)
                    gTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGttt                                                                                                                             <  1:3046826/148‑1 (MQ=255)
                     tAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTc                                                                                                                            >  2:1338330/1‑148 (MQ=255)
                      aaCCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTcc                                                                                                                           <  2:2738487/148‑1 (MQ=255)
                        cccTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTg                                                                                                                         >  1:116187/1‑148 (MQ=255)
                        cccTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTg                                                                                                                         >  1:2719525/1‑148 (MQ=255)
                         ccTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGc                                                                                                                        <  1:1085151/148‑1 (MQ=255)
                           ttGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTa                                                                                                                      <  2:630201/148‑1 (MQ=255)
                           ttGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTa                                                                                                                      <  2:2384479/148‑1 (MQ=255)
                           ttGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTa                                                                                                                      <  2:425807/148‑1 (MQ=255)
                               gTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCa                                                                                                                  <  2:3037191/148‑1 (MQ=255)
                                tCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAt                                                                                                                 <  2:2565396/148‑1 (MQ=255)
                                tCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAt                                                                                                                 >  2:2704785/1‑148 (MQ=255)
                                 cTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAta                                                                                                                <  1:807241/148‑1 (MQ=255)
                                 cTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAta                                                                                                                <  2:466206/148‑1 (MQ=255)
                                 cTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAta                                                                                                                >  1:263905/1‑148 (MQ=255)
                                 cTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAta                                                                                                                <  2:2260774/148‑1 (MQ=255)
                                     tCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgt                                                                                                            <  1:738659/148‑1 (MQ=255)
                                       gcggcgAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgtag                                                                                                          <  2:946018/148‑1 (MQ=255)
                                          gcgAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGc                                                                                                       <  1:1410045/148‑1 (MQ=255)
                                           cgAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGcc                                                                                                      >  1:2220902/1‑148 (MQ=255)
                                              aTCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtta                                                                                                   >  1:1256718/1‑148 (MQ=255)
                                              aTCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtta                                                                                                   <  2:2044583/148‑1 (MQ=255)
                                               tctcTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCttat                                                                                                  <  2:1225127/148‑1 (MQ=255)
                                                 tctTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCttatta                                                                                                <  1:2892772/148‑1 (MQ=255)
                                                  ctTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTaa                                                                                               <  1:3553093/148‑1 (MQ=255)
                                                   ttGCGGATGGTCCGTTACTGTCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACTATATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAAt                                                                                              >  2:1431137/1‑148 (MQ=255)
                                                   ttGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAAt                                                                                              <  1:472933/148‑1 (MQ=255)
                                                    tGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAAtt                                                                                             >  2:1169733/1‑148 (MQ=255)
                                                      cGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTa                                                                                           >  2:399884/1‑148 (MQ=255)
                                                       ggATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTaa                                                                                          <  2:2885465/148‑1 (MQ=255)
                                                       ggATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTaa                                                                                          <  1:1842273/148‑1 (MQ=255)
                                                        gATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAAc                                                                                         <  1:1020399/148‑1 (MQ=255)
                                                         aTGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACt                                                                                        >  1:1344694/1‑148 (MQ=255)
                                                            gTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTg                                                                                     <  1:2008783/148‑1 (MQ=255)
                                                             tCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTgg                                                                                    <  2:2125897/148‑1 (MQ=255)
                                                             tCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTgg                                                                                    >  2:2173307/1‑148 (MQ=255)
                                                             tCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTAGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTTTAGTGCCTTATTAATTTAACTTTgg                                                                                    >  2:2058970/1‑148 (MQ=255)
                                                                       gCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTc                                                                          <  2:1348353/148‑1 (MQ=255)
                                                                       gCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTc                                                                          >  1:508028/1‑148 (MQ=255)
                                                                       gCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTc                                                                          <  1:1783762/148‑1 (MQ=255)
                                                                         ggTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTccc                                                                        <  1:952354/148‑1 (MQ=255)
                                                                              tggctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAAtt                                                                   <  1:164272/148‑1 (MQ=255)
                                                                                gctgTGGCGTACCTTGCAGTTCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGt                                                                 >  1:1952281/1‑148 (MQ=255)
                                                                                gctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGt                                                                 >  1:528297/1‑148 (MQ=255)
                                                                                gctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGt                                                                 >  2:89768/1‑148 (MQ=255)
                                                                                gctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGt                                                                 >  1:2782667/1‑148 (MQ=255)
                                                                                 ctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGtt                                                                <  1:1922306/148‑1 (MQ=255)
                                                                                  tgtgGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTa                                                               <  2:2068579/148‑1 (MQ=255)
                                                                                   gtgGCGTACCTTGCAGTGCGCTACAGGCACTAAACATCAACACCAATAATAAACTTTCAAACAGGTAAATAACGCTATTACTTTTCCTGCTAACCATAGTTTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTAc                                                              <  2:3444019/148‑1 (MQ=255)
                                                                                       cGTACCTTGCATTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACATGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTgg                                                          <  2:533520/148‑1 (MQ=255)
                                                                                        gTACCTTGCAGTGCGCTACAGGCACATAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGt                                                         <  2:2977826/148‑1 (MQ=255)
                                                                                         tACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTa                                                        <  2:1221527/148‑1 (MQ=255)
                                                                                          aCCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTAc                                                       >  1:787992/1‑148 (MQ=255)
                                                                                            cTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTAcgc                                                     >  1:2174421/1‑148 (MQ=255)
                                                                                               gCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGAt                                                  <  1:2271705/148‑1 (MQ=255)
                                                                                                cAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGAtt                                                 >  1:707551/1‑148 (MQ=255)
                                                                                                  gTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGAtttt                                               >  2:1565714/1‑148 (MQ=255)
                                                                                                  gTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGAtttt                                               >  2:841074/1‑148 (MQ=255)
                                                                                                    gcgcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCg                                             <  1:1705820/148‑1 (MQ=255)
                                                                                                    gcgcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCg                                             <  2:444601/148‑1 (MQ=255)
                                                                                                          cAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGaa                                       >  2:377394/1‑148 (MQ=255)
                                                                                                            ggCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGaaaa                                     >  1:2799643/1‑148 (MQ=255)
                                                                                                              cACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACg                                   >  2:3135172/1‑148 (MQ=255)
                                                                                                                cTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTg                                 <  1:3142611/148‑1 (MQ=255)
                                                                                                                 ttAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGa                                <  1:1124450/148‑1 (MQ=255)
                                                                                                                  tAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGAc                               >  2:2438258/1‑148 (MQ=255)
                                                                                                                    aCATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGa                             >  2:2656541/1‑148 (MQ=255)
                                                                                                                     cATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGaa                            >  1:1744430/1‑148 (MQ=255)
                                                                                                                      aTCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGaaa                           >  2:1456075/1‑148 (MQ=255)
                                                                                                                       tCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAg                          >  2:3327119/1‑148 (MQ=255)
                                                                                                                       tCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAg                          >  2:550452/1‑148 (MQ=255)
                                                                                                                        cAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGc                         >  1:596506/1‑148 (MQ=255)
                                                                                                                        cAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGc                         <  1:123519/148‑1 (MQ=255)
                                                                                                                        cAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGc                         >  1:3261510/1‑148 (MQ=255)
                                                                                                                         aaCACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCg                        <  1:1831089/148‑1 (MQ=255)
                                                                                                                             ccAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAAc                    >  2:1044813/1‑148 (MQ=255)
                                                                                                                             ccAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAAc                    >  1:847281/1‑148 (MQ=255)
                                                                                                                              cAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAAcc                   <  1:1888264/148‑1 (MQ=255)
                                                                                                                                 taataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCcac                >  1:2659920/1‑148 (MQ=255)
                                                                                                                                  aataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCcaca               >  2:1801470/1‑148 (MQ=255)
                                                                                                                                    taaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACAtc             >  2:3546211/1‑148 (MQ=255)
                                                                                                                                        cTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTATCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCaa         >  1:2889842/1‑148 (MQ=255)
                                                                                                                                           gCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACAt      >  1:1070012/1‑148 (MQ=255)
                                                                                                                                             aaaCAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTa    <  1:420578/148‑1 (MQ=255)
                                                                                                                                               aCCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATg  >  1:1366538/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GCTTACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTATG  >  CP009273/4406408‑4406698

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