Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,084,798 C→A 100% L102L (CTC→CTA)  rsmE → 16S rRNA m(3)U1498 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,084,7980CA98.4% 565.1 / 0.6 182L102L (CTC→CTA) rsmE16S rRNA m(3)U1498 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (98/81);  total (100/82)
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.22e-01
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 182 total aligned reads displayed.

TTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATC  >  CP009273/3084653‑3084942
                                                                                                                                                 |                                                                                                                                                
ttACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAgg                                                                                                                                                <  2:1422991/148‑1 (MQ=255)
ttACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAgg                                                                                                                                                <  2:640770/148‑1 (MQ=255)
 tACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGgt                                                                                                                                               <  1:3056989/148‑1 (MQ=255)
  aCCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGACATCGATCGAACTAGgtg                                                                                                                                              >  2:1323576/1‑148 (MQ=255)
  aCCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGgtg                                                                                                                                              <  1:396485/148‑1 (MQ=255)
  aCCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGgtg                                                                                                                                              >  2:3066811/1‑148 (MQ=255)
   ccAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGgtgt                                                                                                                                             <  1:3053956/148‑1 (MQ=255)
    cAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAACATGGAATTTACTATCCAGAACTCGATCGAACTAGGTGTa                                                                                                                                            >  2:2583475/1‑148 (MQ=255)
       cgcCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGc                                                                                                                                         >  2:2250640/1‑148 (MQ=255)
       cgcCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGc                                                                                                                                         >  2:1312655/1‑148 (MQ=255)
         ccAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCt                                                                                                                                       <  2:2717499/148‑1 (MQ=255)
          cAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTc                                                                                                                                      <  2:404396/148‑1 (MQ=255)
           aGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCa                                                                                                                                     <  1:3255100/148‑1 (MQ=255)
           aGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCa                                                                                                                                     <  1:1513678/148‑1 (MQ=255)
            gCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCAt                                                                                                                                    >  2:507083/1‑148 (MQ=255)
             cAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCAtt                                                                                                                                   <  1:1802583/148‑1 (MQ=255)
             cAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCAtt                                                                                                                                   >  2:283051/1‑148 (MQ=255)
              aaaaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTa                                                                                                                                  <  2:3138633/148‑1 (MQ=255)
                aaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGACATCGATCGAACTAGGTGTAAGCCTCATTACg                                                                                                                                >  2:323315/1‑148 (MQ=255)
                aaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGGTCGAACTAGGTGTAAGCCTCATTACg                                                                                                                                <  1:2686145/148‑1 (MQ=255)
                aaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACg                                                                                                                                >  1:1010773/1‑148 (MQ=255)
                   aaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCa                                                                                                                             >  1:1928593/1‑148 (MQ=255)
                    aGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCAc                                                                                                                            >  2:1116539/1‑148 (MQ=255)
                     gCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGNTGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACt                                                                                                                           <  1:1781021/148‑1 (MQ=255)
                     gCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACt                                                                                                                           <  1:4556/148‑1 (MQ=255)
                       gTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACttt                                                                                                                         <  1:2375315/148‑1 (MQ=255)
                       gTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACttt                                                                                                                         <  1:2375205/148‑1 (MQ=255)
                        tGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACtttt                                                                                                                        >  1:1132516/1‑148 (MQ=255)
                           aaGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTc                                                                                                                     >  2:3136837/1‑148 (MQ=255)
                                 aGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgcg                                                                                                               <  2:1196897/148‑1 (MQ=255)
                                 aGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgcg                                                                                                               <  1:1568429/148‑1 (MQ=255)
                                  ggTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgcgc                                                                                                              <  1:3176461/148‑1 (MQ=255)
                                         gAAGGCCTGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTgcggcg                                                                                                       <  1:1828296/148‑1 (MQ=255)
                                         gAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTgcggcg                                                                                                       >  2:50098/1‑148 (MQ=255)
                                          aaGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGt                                                                                                      <  2:684972/148‑1 (MQ=255)
                                           aGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGtt                                                                                                     >  1:1275903/1‑148 (MQ=255)
                                            ggCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTa                                                                                                    >  2:1616530/1‑148 (MQ=255)
                                             gCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTaa                                                                                                   <  2:2744141/148‑1 (MQ=255)
                                              ccAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATCTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTaaa                                                                                                  <  1:1912159/148‑1 (MQ=255)
                                                 gATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTg                                                                                               >  1:2090459/1‑148 (MQ=255)
                                                  aTCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTgg                                                                                              >  2:335280/1‑148 (MQ=255)
                                                  aTCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGGGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTgg                                                                                              >  1:2154417/1‑148 (MQ=255)
                                                   tCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGa                                                                                             >  2:2004884/1‑148 (MQ=255)
                                                    cgacgaTCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGAt                                                                                            >  2:1303612/1‑148 (MQ=255)
                                                      acgaTCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAg                                                                                          >  1:990862/1‑148 (MQ=255)
                                                         aTCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGa                                                                                       >  1:3363179/1‑148 (MQ=255)
                                                         aTCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGa                                                                                       >  2:33797/1‑148 (MQ=255)
                                                          tCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGaa                                                                                      >  2:1476394/1‑148 (MQ=255)
                                                              gAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTc                                                                                  <  2:1000392/148‑1 (MQ=255)
                                                                 tctcCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGa                                                                               >  2:803997/1‑148 (MQ=255)
                                                                 tctcCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGa                                                                               >  1:3118286/1‑148 (MQ=255)
                                                                    ccGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACa                                                                            >  2:2763666/1‑148 (MQ=255)
                                                                      gCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaag                                                                          <  2:1928593/148‑1 (MQ=255)
                                                                         gCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaagaag                                                                       >  2:2933401/1‑148 (MQ=255)
                                                                         gCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaagaag                                                                       >  1:293625/1‑148 (MQ=255)
                                                                                  ccTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTcagcag                                                              >  1:2783579/1‑148 (MQ=255)
                                                                                  ccTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTcagcag                                                              <  1:585157/148‑1 (MQ=255)
                                                                                     cGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTgg                                                           >  1:3386001/1‑148 (MQ=255)
                                                                                        tCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCag                                                        >  1:847086/1‑148 (MQ=255)
                                                                                        tCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCag                                                        >  2:1829441/1‑148 (MQ=255)
                                                                                         cAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCaga                                                       <  2:1010773/148‑1 (MQ=255)
                                                                                            gTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaaga                                                    <  1:3385859/148‑1 (MQ=255)
                                                                                            gTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaaga                                                    >  2:2361413/1‑148 (MQ=255)
                                                                                             tGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGAt                                                   >  1:1434398/1‑148 (MQ=255)
                                                                                              gATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGAtt                                                  >  1:1034366/1‑148 (MQ=255)
                                                                                                 gTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCa                                               >  1:2431390/1‑148 (MQ=255)
                                                                                                  tCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCaa                                              <  2:220581/148‑1 (MQ=255)
                                                                                                  tCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCaa                                              <  1:750735/148‑1 (MQ=255)
                                                                                                     cgTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtg                                           >  2:760462/1‑148 (MQ=255)
                                                                                                        ggtgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtgctg                                        <  1:2763666/148‑1 (MQ=255)
                                                                                                         gtgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtgctgc                                       >  1:1923045/1‑148 (MQ=255)
                                                                                                         gtgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtgctgc                                       >  2:2684013/1‑148 (MQ=255)
                                                                                                         gtgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtgctgc                                       >  1:826319/1‑148 (MQ=255)
                                                                                                          tgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGcc                                      >  1:2755146/1‑148 (MQ=255)
                                                                                                          tgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGcc                                      <  2:3118286/148‑1 (MQ=255)
                                                                                                           gAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCt                                     <  2:2320958/148‑1 (MQ=255)
                                                                                                           gAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCt                                     <  1:2266537/148‑1 (MQ=255)
                                                                                                            aaaaaaTGGAATTTACTATCCAGAAATCGATCGAACGAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCtg                                    <  1:168400/148‑1 (MQ=255)
                                                                                                             aaaaaTTGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCtgt                                   <  2:2632841/148‑1 (MQ=255)
                                                                                                             aaaaaTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCtgt                                   <  1:948369/148‑1 (MQ=255)
                                                                                                                  tGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCa                              <  2:2386449/148‑1 (MQ=255)
                                                                                                                   ggAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAg                             >  1:902002/1‑148 (MQ=255)
                                                                                                                     aaTTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAgtg                           <  1:374588/148‑1 (MQ=255)
                                                                                                                     aaTTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAgtg                           <  1:1957378/148‑1 (MQ=255)
                                                                                                                      aTTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAgtgt                          <  2:280433/148‑1 (MQ=255)
                                                                                                                      aTTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAgtgt                          >  1:2157148/1‑148 (MQ=255)
                                                                                                                       tttACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGtgtg                         >  1:3262473/1‑148 (MQ=255)
                                                                                                                       tttACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGtgtg                         >  2:2699404/1‑148 (MQ=255)
                                                                                                                         tACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGt                       >  2:2362274/1‑148 (MQ=255)
                                                                                                                             aTCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTa                   <  2:842499/148‑1 (MQ=255)
                                                                                                                              tCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTaa                  <  1:1669922/148‑1 (MQ=255)
                                                                                                                              tCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATAGCAATTGCTGCCGGTGAGCAGTGTGGTCGTaa                  >  2:2672192/1‑148 (MQ=255)
                                                                                                                               ccAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTaa                  >  2:922418/1‑147 (MQ=255)
                                                                                                                               ccAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAAc                 >  2:360895/1‑148 (MQ=255)
                                                                                                                               ccAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAAc                 >  2:3476986/1‑148 (MQ=255)
                                                                                                                                   aaaTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCggg             >  1:1638169/1‑148 (MQ=255)
                                                                                                                                   aaaTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCggg             <  1:1503430/148‑1 (MQ=255)
                                                                                                                                       cgatcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGcc         >  2:2130799/1‑148 (MQ=255)
                                                                                                                                          tcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGa      >  2:872436/1‑148 (MQ=255)
                                                                                                                                              aCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATc  <  1:2945660/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATC  >  CP009273/3084653‑3084942

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