Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,084,798 C→A 87.1% 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,7980CA87.1% 475.6 / 74.0 208L102L (CTC→CTArsmE16S rRNA m(3)U1498 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (12/15);  new base A (100/81);  total (112/96)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 208 total aligned reads displayed.

CCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATCC  >  CP009273/3084656‑3084943
                                                                                                                                              |                                                                                                                                                 
 cAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGACATCGATCGAACTAGGTGTa                                                                                                                                             >  1:84741/1‑148 (MQ=255)
    cgcCCGCACACAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGc                                                                                                                                          <  2:1822144/148‑1 (MQ=255)
    cgcCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGc                                                                                                                                          <  2:1822634/148‑1 (MQ=255)
    cgcCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGc                                                                                                                                          >  1:415257/1‑148 (MQ=255)
          cAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCAtt                                                                                                                                    <  2:2967631/148‑1 (MQ=255)
          cAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCAtt                                                                                                                                    >  2:46142/1‑148 (MQ=255)
            aaaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTAc                                                                                                                                  >  2:1494786/1‑148 (MQ=255)
             aaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACg                                                                                                                                 <  2:1592012/148‑1 (MQ=255)
              aaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGc                                                                                                                                >  2:861502/1‑148 (MQ=255)
                aaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCa                                                                                                                              <  2:2289172/148‑1 (MQ=255)
                  gCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACt                                                                                                                            <  2:1116839/148‑1 (MQ=255)
                     tGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACtttt                                                                                                                         >  2:1179980/1‑148 (MQ=255)
                       gAAGTGAAGGTGCTGGTAGGCCAGATCGACGATCGCGATTCTCGGCTGCATATTCTCCTCGGTCTGGTGATGTCGCGTGGTGTAAAAATGGAATTTACTATCCAGATATCGATCGAACTAGGTGTAAGCCTCATTACGCCACtttttt                                                                                                                       <  2:540526/148‑1 (MQ=255)
                       gAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACtttttt                                                                                                                       <  2:2946433/148‑1 (MQ=255)
                        aaGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTc                                                                                                                      <  1:3266265/148‑1 (MQ=255)
                          gTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTg                                                                                                                    >  2:23693/1‑148 (MQ=255)
                              aGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgcg                                                                                                                <  1:346115/148‑1 (MQ=255)
                               ggTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAACTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgcgc                                                                                                               >  2:2419761/1‑148 (MQ=255)
                               ggTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgcgc                                                                                                               <  1:3126553/148‑1 (MQ=255)
                                      gAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTgcggcg                                                                                                        >  1:882319/1‑148 (MQ=255)
                                      gAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTgcggcg                                                                                                        >  1:3041600/1‑148 (MQ=255)
                                      gAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTgcggcg                                                                                                        <  2:1839726/148‑1 (MQ=255)
                                       aaGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGt                                                                                                       <  1:1103522/148‑1 (MQ=255)
                                       aaGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGt                                                                                                       <  1:622525/148‑1 (MQ=255)
                                       aaGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGt                                                                                                       <  1:2561984/148‑1 (MQ=255)
                                       aaGGCCAGATCGACGATCGCGAATCTCCGCTGCATATGCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGt                                                                                                       <  1:2562065/148‑1 (MQ=255)
                                        aGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGtt                                                                                                      >  1:2417452/1‑148 (MQ=255)
                                         ggCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTa                                                                                                     >  1:277977/1‑148 (MQ=255)
                                          gCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGACCTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTaa                                                                                                    <  2:1262731/148‑1 (MQ=255)
                                              gATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTg                                                                                                >  1:1797062/1‑148 (MQ=255)
                                              gATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTg                                                                                                >  2:1904736/1‑148 (MQ=255)
                                              gATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTg                                                                                                >  2:234839/1‑148 (MQ=255)
                                              gATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACGAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTg                                                                                                <  1:609939/148‑1 (MQ=255)
                                               aTCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTgg                                                                                               >  1:1974718/1‑148 (MQ=255)
                                               aTCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTgg                                                                                               >  1:1809385/1‑148 (MQ=255)
                                                 cgacgaTCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGAt                                                                                             >  2:1361686/1‑148 (MQ=255)
                                                   acgaTCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAg                                                                                           >  1:1192357/1‑148 (MQ=255)
                                                    cgaTCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGt                                                                                          >  1:1008550/1‑148 (MQ=255)
                                                      aTCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGa                                                                                        <  2:2770969/148‑1 (MQ=255)
                                                       tCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGACAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGaa                                                                                       <  2:3099663/148‑1 (MQ=255)
                                                       tCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGaa                                                                                       >  1:3134606/1‑148 (MQ=255)
                                                       tCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGaa                                                                                       <  2:584822/148‑1 (MQ=255)
                                                       tCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGaa                                                                                       <  2:3168496/148‑1 (MQ=255)
                                                       tCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGaa                                                                                       >  1:3223987/1‑148 (MQ=255)
                                                       tCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATGTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGaa                                                                                       <  2:1457074/148‑1 (MQ=255)
                                                          cgAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGt                                                                                    >  2:3276244/1‑148 (MQ=255)
                                                          cgAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGt                                                                                    <  2:961406/148‑1 (MQ=255)
                                                             aTCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTg                                                                                 >  2:544649/1‑148 (MQ=255)
                                                              tctcCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGa                                                                                <  2:1480061/148‑1 (MQ=255)
                                                               ctcCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGaa                                                                               >  2:1919788/1‑148 (MQ=255)
                                                                 ccGCGGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACa                                                                             >  1:3113933/1‑148 (MQ=255)
                                                                   gCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaag                                                                           <  1:675498/148‑1 (MQ=255)
                                                                   gCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaag                                                                           <  1:3023667/148‑1 (MQ=255)
                                                                    cTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaaga                                                                          >  1:2889158/1‑148 (MQ=255)
                                                                       cATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGc                                                                       <  2:908496/148‑1 (MQ=255)
                                                                         tatTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCtt                                                                     >  1:3226389/1‑148 (MQ=255)
                                                                                 tCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTg                                                             <  1:2129555/148‑1 (MQ=255)
                                                                                 tCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTg                                                             <  1:977108/148‑1 (MQ=255)
                                                                                  cGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTgg                                                            <  1:1494786/148‑1 (MQ=255)
                                                                                       agttGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaa                                                       <  2:2607281/145‑1 (MQ=255)
                                                                                        ggTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaag                                                      >  2:2627045/1‑148 (MQ=255)
                                                                                         gTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaaga                                                     <  2:1097404/148‑1 (MQ=255)
                                                                                          tGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGAt                                                    <  2:2589344/148‑1 (MQ=255)
                                                                                              gTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCa                                                <  2:2087596/148‑1 (MQ=255)
                                                                                               tCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCaa                                               <  1:2926406/148‑1 (MQ=255)
                                                                                                  cgTGGTGAAAAAATGGAATTTACTCTCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtg                                            <  2:2719764/148‑1 (MQ=255)
                                                                                                  cgTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtg                                            <  1:1534683/148‑1 (MQ=255)
                                                                                                     ggtgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtgctg                                         <  1:2560768/148‑1 (MQ=255)
                                                                                                      gtgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtgctgc                                        <  1:3202281/148‑1 (MQ=255)
                                                                                                       tgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGcc                                       <  1:536244/148‑1 (MQ=255)
                                                                                                       tgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGcc                                       <  2:3229903/148‑1 (MQ=255)
                                                                                                          aaaaaTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCtgt                                    >  2:3295367/1‑148 (MQ=255)
                                                                                                            aaaTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGa                                  >  2:222277/1‑148 (MQ=255)
                                                                                                            aaaTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGGTAAACTGGATAGTGAACGGCGGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGa                                  >  2:221556/1‑148 (MQ=255)
                                                                                                              aTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTTGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGc                                >  1:3251079/1‑148 (MQ=255)
                                                                                                               tGGAATTTACTATCCAGAACTCGATCGAACTAGCTGTAAGCCTCATTACGCCACTTTTATCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCa                               >  1:2946482/1‑148 (MQ=255)
                                                                                                               tGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCa                               <  1:269465/148‑1 (MQ=255)
                                                                                                               tGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCa                               >  1:509267/1‑148 (MQ=255)
                                                                                                                ggAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAAAGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAg                              >  1:1742348/1‑148 (MQ=255)
                                                                                                                    tttACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGtgtg                          >  1:2977504/1‑148 (MQ=255)
                                                                                                                        cTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCg                      >  1:1610281/1‑148 (MQ=255)
                                                                                                                        cTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCg                      >  2:1216346/1‑148 (MQ=255)
                                                                                                                          ttCCAGAAAGCGATAGCAATCGGGGTAAGCCTCATTACGCCAGTTTTTTCTGAGCTCTGAGTCGTTAAACTGGATAGTGAACGTCTGAACAATAAGCTTCAGCTGTGGCAGATGATTTCAATTGCTGCCTGTGAGCAGTGTGGTCGTa                    <  1:595365/147‑1 (MQ=255)
                                                                                                                          aTCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTa                    <  1:46142/148‑1 (MQ=255)
                                                                                                                          aTCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTGAAACTGGATAGGGAACGTCTGAACAAGAAGCTTCAGAAGTGGCAGAAGATTGCAATAGCTGCCTGTGAGCAGTGTGGTCGTa                    >  1:1446048/1‑148 (MQ=255)
                                                                                                                            ccAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAAc                  <  2:1989603/148‑1 (MQ=255)
                                                                                                                               gAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCgg               >  1:835360/1‑148 (MQ=255)
                                                                                                                               gAAATCGATCGAACTAGGTGTAAGCCTCATTACGACACTTATTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGATGCCTGTGAGCAGTGTGGTCGTAACCgg               >  2:435481/1‑148 (MQ=255)
                                                                                                                                aaaTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCggg              <  1:1899059/148‑1 (MQ=255)
                                                                                                                                aaaTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCggg              >  2:2703392/1‑148 (MQ=255)
                                                                                                                                aaaTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCggg              >  2:385692/1‑148 (MQ=255)
                                                                                                                                aaaTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCggg              >  2:2696107/1‑148 (MQ=255)
                                                                                                                                    cgatcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGcc          >  1:2591468/1‑148 (MQ=255)
                                                                                                                                    caatcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGcc          >  2:2030153/3‑148 (MQ=255)
                                                                                                                                    caatcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGcc          >  2:2029980/3‑148 (MQ=255)
                                                                                                                                       tcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGa       >  2:659100/1‑148 (MQ=255)
                                                                                                                                        cgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGaa      <  1:861518/148‑1 (MQ=255)
                                                                                                                                         gaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGaaa     >  1:2468869/1‑148 (MQ=255)
                                                                                                                                           aCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATc   >  2:1925178/1‑148 (MQ=255)
                                                                                                                                            cTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATcc  >  1:2298310/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
CCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATCC  >  CP009273/3084656‑3084943

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: