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,7980CA95.7% 616.9 / 15.4 209L102L (CTC→CTArsmE16S rRNA m(3)U1498 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (6/4);  new base A (102/97);  total (108/101)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 209 total aligned reads displayed.

CCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATC  >  CP009273/3084656‑3084942
                                                                                                                                              |                                                                                                                                                
 cAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAACTCGATCGAACTAGGTGTa                                                                                                                                            >  1:1304894/1‑148 (MQ=255)
   gcgcCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAg                                                                                                                                          >  2:291456/1‑148 (MQ=255)
       catcAAAAAAATCGTGGAATTGAAGGTGCTGGAATGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCTCCTGGTTCAGGTGATTTCGCGTGGTGAAATAATTGATTTTACTATCAAGTAATCGATCGAACTAGGTGTAAGCCTc                                                                                                                                      <  2:3271747/145‑1 (MQ=255)
       cAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTc                                                                                                                                      <  2:770606/148‑1 (MQ=255)
        aGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCa                                                                                                                                     <  1:1132702/148‑1 (MQ=255)
          cAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCAtt                                                                                                                                   <  2:3020893/148‑1 (MQ=255)
          cAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCAtt                                                                                                                                   <  1:2884051/148‑1 (MQ=255)
           aaaaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTa                                                                                                                                  <  1:1279165/148‑1 (MQ=255)
            aaaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTAc                                                                                                                                 >  2:6588/1‑148 (MQ=255)
            aaaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTAc                                                                                                                                 >  2:1448879/1‑148 (MQ=255)
              aaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGc                                                                                                                               >  2:401262/1‑148 (MQ=255)
              aaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGc                                                                                                                               >  2:2962670/1‑148 (MQ=255)
                 aGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCAc                                                                                                                            >  2:2396579/1‑148 (MQ=255)
                 aGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCAc                                                                                                                            >  1:2037908/1‑148 (MQ=255)
                 aGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCAc                                                                                                                            >  1:67647/1‑148 (MQ=255)
                 aGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCAc                                                                                                                            >  1:107960/1‑148 (MQ=255)
                  gCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACt                                                                                                                           >  2:2510253/1‑148 (MQ=255)
                  gCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACt                                                                                                                           >  1:610135/1‑148 (MQ=255)
                   cGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACtt                                                                                                                          >  1:770541/1‑148 (MQ=255)
                    gTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCAACTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACttt                                                                                                                         <  2:3665792/148‑1 (MQ=255)
                     tGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACtttt                                                                                                                        <  1:1623026/148‑1 (MQ=255)
                      ggAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTCTAAGCCTCATTACGCCACttttt                                                                                                                       >  2:1974307/1‑148 (MQ=255)
                       gAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACtttttt                                                                                                                      >  2:2644770/1‑148 (MQ=255)
                       gAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACtttttt                                                                                                                      >  2:727731/1‑148 (MQ=255)
                        aaGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTc                                                                                                                     <  2:230599/148‑1 (MQ=255)
                          gTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTg                                                                                                                   <  1:1977888/148‑1 (MQ=255)
                          gTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTg                                                                                                                   <  1:2262814/148‑1 (MQ=255)
                           tGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGa                                                                                                                  >  2:3141702/1‑148 (MQ=255)
                             aaGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgc                                                                                                                >  2:2469340/1‑148 (MQ=255)
                               ggTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgcgc                                                                                                              <  2:3424861/148‑1 (MQ=255)
                                 tGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTg                                                                                                            <  1:2793796/148‑1 (MQ=255)
                                   cTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTgcg                                                                                                          <  1:1642402/148‑1 (MQ=255)
                                      gAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTgcggcg                                                                                                       >  1:375590/1‑148 (MQ=255)
                                      gAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTgcggcg                                                                                                       <  1:3750361/148‑1 (MQ=255)
                                         ggCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTa                                                                                                    >  1:679352/1‑148 (MQ=255)
                                          gCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTaa                                                                                                   >  1:2727419/1‑148 (MQ=255)
                                           ccAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTaaa                                                                                                  <  2:1395362/148‑1 (MQ=255)
                                           ccAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTaaa                                                                                                  <  2:1642307/148‑1 (MQ=255)
                                           ccAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTaaa                                                                                                  >  1:129678/1‑148 (MQ=255)
                                           ccAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTaaa                                                                                                  <  1:2302514/148‑1 (MQ=255)
                                           ccAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTaaa                                                                                                  <  1:378856/148‑1 (MQ=255)
                                            cAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAAc                                                                                                 >  1:821141/1‑148 (MQ=255)
                                            cAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAAc                                                                                                 <  2:3479111/148‑1 (MQ=255)
                                              gATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTg                                                                                               <  1:3708712/148‑1 (MQ=255)
                                                  gacgaTCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATa                                                                                           <  2:1195895/148‑1 (MQ=255)
                                                       tCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGaa                                                                                      <  2:1561318/148‑1 (MQ=255)
                                                         gcgAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACg                                                                                    <  1:3115622/148‑1 (MQ=255)
                                                                    cttcATATTCACCTCGGTCAGGTGATGTCGCGTGTTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaaga                                                                         <  1:3574958/145‑1 (MQ=255)
                                                                    cTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaaga                                                                         <  1:3814112/148‑1 (MQ=255)
                                                                    cTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaaga                                                                         <  2:560872/148‑1 (MQ=255)
                                                                      gCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGACCGTCTGAACaagaag                                                                       >  1:1046018/1‑148 (MQ=255)
                                                                         tatTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAACTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCtt                                                                    >  1:2973696/1‑148 (MQ=255)
                                                                           ttCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTca                                                                  <  2:744519/148‑1 (MQ=255)
                                                                            tCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTcag                                                                 >  2:2915329/1‑148 (MQ=255)
                                                                               ccTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTcagcag                                                              <  1:3098972/148‑1 (MQ=255)
                                                                                cTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGt                                                             >  1:2704262/1‑148 (MQ=255)
                                                                                  cGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTgg                                                           <  1:363910/148‑1 (MQ=255)
                                                                                    gTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCa                                                         <  1:926108/148‑1 (MQ=255)
                                                                                        ggTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaag                                                     >  1:2353940/1‑148 (MQ=255)
                                                                                        ggTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaag                                                     >  1:634914/1‑148 (MQ=255)
                                                                                          tGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGAt                                                   >  1:3154099/1‑148 (MQ=255)
                                                                                              gTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCa                                               <  2:435661/148‑1 (MQ=255)
                                                                                                cgcgTGGTGAAAAAATGGAATTTACTATCCAGAAATCTATCGAACTAGGTGTAAGCCTCATTACTCCACTTTTTTCTGAGCGCTGCGGCGTTAAACAGTATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAAt                                             >  2:3526546/1‑148 (MQ=255)
                                                                                                cgcgTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAAt                                             >  2:3691609/1‑148 (MQ=255)
                                                                                                cgcgTGGTGAAAAAATGGAATTTACTATACATATATCGATCGACCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGTCGTTAAACTGGATAGTGAACGTCTGAACTATAAGCTTCAGCAGTGGCAGAAGATTGCAAt                                             <  1:3644583/148‑1 (MQ=255)
                                                                                                  cgTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtg                                           <  1:1448879/148‑1 (MQ=255)
                                                                                                      gtgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtgctgc                                       <  2:130014/148‑1 (MQ=255)
                                                                                                       tgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGcc                                      <  1:3597190/148‑1 (MQ=255)
                                                                                                           aaaaTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCtgtg                                  >  1:3540643/1‑148 (MQ=255)
                                                                                                            aaaTGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGa                                 >  2:3123132/1‑148 (MQ=255)
                                                                                                            aaaTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGa                                 >  1:946239/1‑148 (MQ=255)
                                                                                                                ggAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAg                             <  1:6588/148‑1 (MQ=255)
                                                                                                                  aaTTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAgtg                           <  1:1202367/148‑1 (MQ=255)
                                                                                                                    tttACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGtgtg                         >  1:2963910/1‑148 (MQ=255)
                                                                                                                     ttACTCTCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTgg                        <  2:2652595/148‑1 (MQ=255)
                                                                                                                     ttACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGAAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTgg                        >  1:2916132/1‑148 (MQ=255)
                                                                                                                        cTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCg                     <  1:2621032/148‑1 (MQ=255)
                                                                                                                          aTCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTa                   >  2:2154882/1‑148 (MQ=255)
                                                                                                                           tCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGACGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTaa                  >  1:200793/1‑148 (MQ=255)
                                                                                                                           tCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTaa                  >  2:3385457/1‑148 (MQ=255)
                                                                                                                            ccAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAAc                 >  1:1616548/1‑148 (MQ=255)
                                                                                                                            ccAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAAc                 >  1:1896815/1‑148 (MQ=255)
                                                                                                                                 aaTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGt            <  2:2727419/148‑1 (MQ=255)
                                                                                                                                  atcgatcgAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTCCCCGGGTg           >  1:2299171/1‑148 (MQ=255)
                                                                                                                                   tcgatcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGc          >  2:597593/1‑148 (MQ=255)
                                                                                                                                    cgatcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGAGTCAGAAGATAGCAATTGCTGCCTGTGAGCAGAGAGGTCGTAACCGGGTGcc         >  2:1743703/1‑148 (MQ=255)
                                                                                                                                      atcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCgg       >  2:1770380/1‑148 (MQ=255)
                                                                                                                                      atcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCgg       >  2:1768960/1‑148 (MQ=255)
                                                                                                                                       tcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGa      >  2:2334548/1‑148 (MQ=255)
                                                                                                                                       tcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGa      >  2:194669/1‑148 (MQ=255)
                                                                                                                                        cgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCGGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGaa     <  1:2666515/148‑1 (MQ=255)
                                                                                                                                          aaCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAAt   <  1:2443230/148‑1 (MQ=255)
                                                                                                                                          aaCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAAt   <  1:1485721/148‑1 (MQ=255)
                                                                                                                                          aaCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAAt   <  1:544884/148‑1 (MQ=255)
                                                                                                                                          aaCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAAt   <  1:3392164/148‑1 (MQ=255)
                                                                                                                                          aaCTAGGTGTAAGCCTCAATACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAAt   <  2:2954478/148‑1 (MQ=255)
                                                                                                                                           aCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATc  >  2:3249663/1‑148 (MQ=255)
                                                                                                                                           aCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATc  >  2:1881688/1‑148 (MQ=255)
                                                                                                                                           aCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATc  >  2:1130984/1‑148 (MQ=255)
                                                                                                                                           aCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATc  >  2:857870/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATC  >  CP009273/3084656‑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: