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,5530CA98.0% 799.1 / 6.0 252R12L (CGG→CTG) bsmAbioflm peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base A (87/159);  total (91/161)
Fisher's exact test for biased strand distribution p-value = 5.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01
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 252 total aligned reads displayed.

GCTTACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTAT  >  CP009273/4406408‑4406697
                                                                                                                                                 |                                                                                                                                                
gCTTACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAgg                                                                                                                                                <  1:3665941/148‑1 (MQ=255)
      gCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATa                                                                                                                                          <  2:1384599/148‑1 (MQ=255)
       cTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATaa                                                                                                                                         <  1:415473/148‑1 (MQ=255)
            tATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTa                                                                                                                                    <  1:3016014/148‑1 (MQ=255)
            tATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTa                                                                                                                                    >  2:328016/1‑148 (MQ=255)
              ttCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTAtt                                                                                                                                  <  2:1378752/148‑1 (MQ=255)
               tCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTa                                                                                                                                 >  2:3595959/1‑148 (MQ=255)
               tCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTa                                                                                                                                 >  2:3548374/1‑148 (MQ=255)
                  cTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGt                                                                                                                              <  2:278763/148‑1 (MQ=255)
                   tGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGtt                                                                                                                             <  1:106529/148‑1 (MQ=255)
                    gTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGttt                                                                                                                            <  1:3432308/148‑1 (MQ=255)
                    gTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGttt                                                                                                                            <  2:2798817/148‑1 (MQ=255)
                    gTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGttt                                                                                                                            <  1:1403583/148‑1 (MQ=255)
                     tAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTc                                                                                                                           <  2:3373906/148‑1 (MQ=255)
                      aaCCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTcc                                                                                                                          >  1:920235/1‑148 (MQ=255)
                      aaCCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTcc                                                                                                                          >  1:2102686/1‑148 (MQ=255)
                         ccTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGc                                                                                                                       >  1:1353246/1‑148 (MQ=255)
                           ttGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTa                                                                                                                     >  1:2310792/1‑148 (MQ=255)
                                tCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAt                                                                                                                <  1:3599020/148‑1 (MQ=255)
                                 cTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAta                                                                                                               <  2:1030429/148‑1 (MQ=255)
                                 cTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAta                                                                                                               <  1:2564858/148‑1 (MQ=255)
                                      cgcgGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgta                                                                                                          <  1:574659/148‑1 (MQ=255)
                                      cgcgGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgta                                                                                                          <  1:2495761/148‑1 (MQ=255)
                                       gcggcgAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagtgtag                                                                                                         <  2:3338547/148‑1 (MQ=255)
                                           cgAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGcc                                                                                                     <  1:3047032/148‑1 (MQ=255)
                                            gAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCt                                                                                                    <  1:3159948/148‑1 (MQ=255)
                                             aaTCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtt                                                                                                   >  1:667610/1‑148 (MQ=255)
                                             aaTCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtt                                                                                                   <  2:1349035/148‑1 (MQ=255)
                                             aaTCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtt                                                                                                   <  2:847651/148‑1 (MQ=255)
                                              aTCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtta                                                                                                  >  2:3037699/1‑148 (MQ=255)
                                               tctcTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCttat                                                                                                 >  1:3325529/1‑148 (MQ=255)
                                               tctcTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCttat                                                                                                 <  2:3683096/148‑1 (MQ=255)
                                                ctctTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCttatt                                                                                                >  2:570501/1‑148 (MQ=255)
                                                 tctTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCttatta                                                                                               >  1:1079039/1‑148 (MQ=255)
                                                  ctTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTaa                                                                                              >  2:3623608/1‑148 (MQ=255)
                                                  ctTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGGCGTGCCTTCTTaa                                                                                              >  2:1961863/1‑148 (MQ=255)
                                                   ttGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAAt                                                                                             <  2:1168490/148‑1 (MQ=255)
                                                    tGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAAtt                                                                                            >  2:1124331/1‑148 (MQ=255)
                                                        gATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAAc                                                                                        >  2:750339/1‑148 (MQ=255)
                                                        gATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAAc                                                                                        >  2:1881525/1‑148 (MQ=255)
                                                        gATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAAc                                                                                        >  2:1895188/1‑148 (MQ=255)
                                                        gATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAAc                                                                                        <  2:2032178/148‑1 (MQ=255)
                                                           ggTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAACCTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACttt                                                                                     <  1:3221688/148‑1 (MQ=255)
                                                            gTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTg                                                                                    <  2:2958810/148‑1 (MQ=255)
                                                             tCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTgg                                                                                   >  1:1184258/1‑148 (MQ=255)
                                                                   aCTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAAt                                                                             >  2:2298407/1‑148 (MQ=255)
                                                                    cTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATg                                                                            <  2:379125/148‑1 (MQ=255)
                                                                    cTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATg                                                                            <  1:495149/148‑1 (MQ=255)
                                                                    cTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATg                                                                            <  2:2734493/148‑1 (MQ=255)
                                                                    cTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATg                                                                            <  2:3234820/148‑1 (MQ=255)
                                                                     tGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGt                                                                           <  1:826747/148‑1 (MQ=255)
                                                                      ggCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGtt                                                                          <  1:2212893/148‑1 (MQ=255)
                                                                       gCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTc                                                                         <  2:1745726/148‑1 (MQ=255)
                                                                           tGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGa                                                                     >  1:198028/1‑148 (MQ=255)
                                                                              tggctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAAtt                                                                  <  2:1418191/148‑1 (MQ=255)
                                                                                gctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGt                                                                <  1:590179/148‑1 (MQ=255)
                                                                                 ctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGtt                                                               <  1:8139/148‑1 (MQ=255)
                                                                                   gtgGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTAc                                                             >  2:3460326/1‑148 (MQ=255)
                                                                                   gtgGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTAc                                                             <  1:1512663/148‑1 (MQ=255)
                                                                                   gtgGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTAc                                                             <  2:1531589/148‑1 (MQ=255)
                                                                                    tgGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTAcc                                                            <  2:2819799/148‑1 (MQ=255)
                                                                                     ggCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCt                                                           <  1:961532/148‑1 (MQ=255)
                                                                                      gagTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTg                                                          <  1:2258663/146‑1 (MQ=255)
                                                                                      gCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTg                                                          <  1:2878870/148‑1 (MQ=255)
                                                                                      gCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTg                                                          <  1:338712/148‑1 (MQ=255)
                                                                                        gTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGt                                                        <  1:1569360/148‑1 (MQ=255)
                                                                                        gTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGt                                                        <  2:1804448/148‑1 (MQ=255)
                                                                                         tACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTa                                                       <  2:145273/148‑1 (MQ=255)
                                                                                         tACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTa                                                       <  2:145162/148‑1 (MQ=255)
                                                                                          aCCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTAc                                                      >  2:2193924/1‑148 (MQ=255)
                                                                                           ccTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTAcg                                                     >  1:3530898/1‑148 (MQ=255)
                                                                                           ccTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTAcg                                                     <  2:807197/148‑1 (MQ=255)
                                                                                             ttGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTAcgcg                                                   <  2:1401100/148‑1 (MQ=255)
                                                                                                cAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGAtt                                                >  1:3270119/1‑148 (MQ=255)
                                                                                                 aGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGAttt                                               <  2:2060099/148‑1 (MQ=255)
                                                                                                           aGGCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATANTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGaaa                                     <  1:1932166/148‑1 (MQ=255)
                                                                                                             gCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAAc                                   >  1:1678987/1‑148 (MQ=255)
                                                                                                             gCACTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTAACTGGTACGCGATTTTCGAAAGaaaaa                                   >  1:2295324/1‑147 (MQ=255)
                                                                                                                cTTAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTg                                <  1:1290416/148‑1 (MQ=255)
                                                                                                                  tAACATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGAc                              >  2:3263376/1‑148 (MQ=255)
                                                                                                                    aCATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGa                            >  2:2584049/1‑148 (MQ=255)
                                                                                                                    aCATCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGa                            >  1:1761439/1‑148 (MQ=255)
                                                                                                                       tCAACACCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAg                         >  1:1530492/1‑148 (MQ=255)
                                                                                                                         aaCACCAATAATAAACTGGCAAACAGGTAAATAACTCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCg                       >  2:1058894/1‑148 (MQ=255)
                                                                                                                           cacCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTa                     >  2:1742863/1‑148 (MQ=255)
                                                                                                                            acCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTaa                    >  1:1245364/1‑148 (MQ=255)
                                                                                                                            acCAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTaa                    <  2:1025342/148‑1 (MQ=255)
                                                                                                                             ccAATAATAAACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAAc                   >  1:3502777/1‑148 (MQ=255)
                                                                                                                              cAATAATAAACTCGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAAcc                  <  1:1445114/148‑1 (MQ=255)
                                                                                                                                 taataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCcac               <  2:235164/148‑1 (MQ=255)
                                                                                                                                  aataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCcaca              <  1:3447996/148‑1 (MQ=255)
                                                                                                                                   ataaCCTGGCACACAGGTAAATATCGATATTACGTTTCATGCTAACCATAGTGTAGTGCCTTATTATTTTAACTTTGGGGGTATGTTCCCGAATTGTTACCTGGTACGCGATTTTCTAATGTAAACGTGTCGAAAGCGTAACACACAt             <  2:1224945/148‑1 (MQ=255)
                                                                                                                                   ataaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACAt             <  2:2102966/148‑1 (MQ=255)
                                                                                                                                    taaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACAtc            <  2:1790745/148‑1 (MQ=255)
                                                                                                                                    taaACTGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACAtc            >  2:1464907/1‑148 (MQ=255)
                                                                                                                                     aaaCTGGCACACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACAtct           <  1:2126570/148‑1 (MQ=255)
                                                                                                                                         tGGCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAAc       <  1:648615/148‑1 (MQ=255)
                                                                                                                                          ggCAAACAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACa      >  1:1563296/1‑148 (MQ=255)
                                                                                                                                              aaCAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTAt  >  2:2658914/1‑148 (MQ=255)
                                                                                                                                              aaCAGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTAt  <  1:328016/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCTTACGCTACCTATTCGCTGTAACCCTTGCGTCTGGTCGCGGCGAATCTCTTGCGGATGGTCCGTTACTGGCGGTGCTGGCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTATTAATTTAACTTTGGGGGAATGTTCCCGAATTGTTACCTGGTACGCGATTTTCGAATGAAAACGTGACGAAAGCGTAACCCACATCTCAACATTAT  >  CP009273/4406408‑4406697

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: