Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,818,010 C→A 100% H46N (CAT→AAT)  trmH → tRNA mG18‑2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,818,0100CA97.8% 617.0 / 3.5 197H46N (CAT→AAT) trmHtRNA mG18‑2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base A (114/78);  total (116/81)
Fisher's exact test for biased strand distribution p-value = 4.04e-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 197 total aligned reads displayed.

TAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTC  >  CP009273/3817866‑3818154
                                                                                                                                                |                                                                                                                                                
  aTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGaa                                                                                                                                             <  2:1268437/148‑1 (MQ=255)
  aTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGaa                                                                                                                                             >  2:1840247/1‑148 (MQ=255)
  aTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGAGTATTCGTACCGCAGATGCCGTTGGCGTAAATGaa                                                                                                                                             >  2:1837679/1‑148 (MQ=255)
   tGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAg                                                                                                                                            <  2:1316113/148‑1 (MQ=255)
   tGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAg                                                                                                                                            >  2:1879655/1‑148 (MQ=255)
    gTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGt                                                                                                                                           >  1:2596708/1‑148 (MQ=255)
    gTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGt                                                                                                                                           >  2:2840539/1‑148 (MQ=255)
      tttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTc                                                                                                                                         <  2:2739648/148‑1 (MQ=255)
      tttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTc                                                                                                                                         >  2:1664166/1‑148 (MQ=255)
      tttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTc                                                                                                                                         >  1:915524/1‑148 (MQ=255)
      tttATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTc                                                                                                                                         >  1:191578/1‑148 (MQ=255)
            aaCCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCg                                                                                                                                   >  2:3150781/1‑148 (MQ=255)
             aCCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGt                                                                                                                                  <  2:2607043/148‑1 (MQ=255)
              cccAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTc                                                                                                                                 <  1:2247916/148‑1 (MQ=255)
              cccAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTc                                                                                                                                 <  1:3006498/148‑1 (MQ=255)
              cccAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTc                                                                                                                                 <  2:102302/148‑1 (MQ=255)
                cAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTg                                                                                                                               <  1:1715613/148‑1 (MQ=255)
                 aaCACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTgg                                                                                                                              >  1:298160/1‑148 (MQ=255)
                  acacGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGc                                                                                                                             <  1:2958192/148‑1 (MQ=255)
                  acacGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGc                                                                                                                             <  1:2717196/148‑1 (MQ=255)
                   cacGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGcc                                                                                                                            <  2:968652/148‑1 (MQ=255)
                    acGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCt                                                                                                                           >  1:1460207/1‑148 (MQ=255)
                    acGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCt                                                                                                                           >  1:3509367/1‑148 (MQ=255)
                     cGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTg                                                                                                                          >  1:1892704/1‑148 (MQ=255)
                      gTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTgg                                                                                                                         <  1:994153/148‑1 (MQ=255)
                      gTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTgg                                                                                                                         <  1:1001176/148‑1 (MQ=255)
                      gTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTgg                                                                                                                         <  1:1001116/148‑1 (MQ=255)
                      gTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTgg                                                                                                                         <  1:1472128/148‑1 (MQ=255)
                       ttATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGt                                                                                                                        >  2:1616814/1‑148 (MQ=255)
                       ttATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGt                                                                                                                        >  1:2312860/1‑148 (MQ=255)
                        tATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTa                                                                                                                       >  2:1946536/1‑148 (MQ=255)
                         aTGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAg                                                                                                                      <  2:2566254/148‑1 (MQ=255)
                          tGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGc                                                                                                                     >  2:1028456/1‑148 (MQ=255)
                             aCGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGc                                                                                                                  >  1:3048081/1‑148 (MQ=255)
                              cGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCa                                                                                                                 >  1:1565858/1‑148 (MQ=255)
                                cATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATg                                                                                                               <  1:6942/148‑1 (MQ=255)
                                    tGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCa                                                                                                           <  2:404/148‑1 (MQ=255)
                                     gCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCAc                                                                                                          >  2:486502/1‑148 (MQ=255)
                                        aaaTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCAt                                                                                                       >  1:2478996/1‑148 (MQ=255)
                                          aTGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATgg                                                                                                     <  2:1180638/148‑1 (MQ=255)
                                            gCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCt                                                                                                   >  2:303177/1‑148 (MQ=255)
                                            gCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCt                                                                                                   >  2:2467/1‑148 (MQ=255)
                                             cTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGTAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCtt                                                                                                  <  2:3074353/148‑1 (MQ=255)
                                              tCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTc                                                                                                 <  1:2169963/148‑1 (MQ=255)
                                               cGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCg                                                                                                >  1:3042016/1‑148 (MQ=255)
                                               cGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCg                                                                                                >  1:6485/1‑148 (MQ=255)
                                               cGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCg                                                                                                >  1:2772941/1‑148 (MQ=255)
                                               cGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGACTGGCCTGGAAGCCGCATGCGCACCATGTCTTCg                                                                                                >  1:6516/1‑148 (MQ=255)
                                                gCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCgg                                                                                               >  2:560278/1‑148 (MQ=255)
                                                   aGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAg                                                                                            >  1:940921/1‑148 (MQ=255)
                                                   aGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAg                                                                                            >  1:1735539/1‑148 (MQ=255)
                                                    ggcggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgc                                                                                           >  1:1291564/1‑148 (MQ=255)
                                                    ggcggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAgc                                                                                           >  2:660395/1‑148 (MQ=255)
                                                     gcggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGc                                                                                              <  1:706589/144‑1 (MQ=255)
                                                     gcggcAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGc                                                                                              >  2:706589/1‑144 (MQ=255)
                                                            ccTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgta                                                                                   >  2:814900/1‑148 (MQ=255)
                                                            ccTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgta                                                                                   <  2:312275/148‑1 (MQ=255)
                                                             cTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgtag                                                                                  >  2:3523366/1‑148 (MQ=255)
                                                               gATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGgtagta                                                                                <  1:3112294/148‑1 (MQ=255)
                                                                aTCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTaa                                                                               <  1:1008263/148‑1 (MQ=255)
                                                                  cTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACa                                                                             >  1:1501087/1‑148 (MQ=255)
                                                                         tCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGt                                                                      >  2:3266086/1‑148 (MQ=255)
                                                                         tCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGt                                                                      <  1:560278/148‑1 (MQ=255)
                                                                                ggAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTg                                                               <  1:1639848/148‑1 (MQ=255)
                                                                                 gAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGa                                                              >  1:2063266/1‑148 (MQ=255)
                                                                                    cAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGaaaa                                                           <  1:2840539/148‑1 (MQ=255)
                                                                                    cAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGaaaa                                                           >  2:634696/1‑148 (MQ=255)
                                                                                     aGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAac                                                          <  2:2215555/148‑1 (MQ=255)
                                                                                      ggTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAaca                                                         <  1:707834/148‑1 (MQ=255)
                                                                                       gTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAacac                                                        <  2:1714545/148‑1 (MQ=255)
                                                                                       gTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAacac                                                        <  1:318300/148‑1 (MQ=255)
                                                                                         ccACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAacacac                                                      >  1:3084235/1‑148 (MQ=255)
                                                                                           acaAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCg                                                    >  1:2640228/1‑148 (MQ=255)
                                                                                           acaAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCg                                                    >  1:2644883/1‑148 (MQ=255)
                                                                                                 cTCATAACGGTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCAt                                              <  2:3237900/148‑1 (MQ=255)
                                                                                                   cATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTg                                            >  1:3054813/1‑148 (MQ=255)
                                                                                                   cATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTg                                            <  2:3341987/148‑1 (MQ=255)
                                                                                                     tAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTTGc                                          >  1:2344985/1‑148 (MQ=255)
                                                                                                     tAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGc                                          >  1:2393034/1‑148 (MQ=255)
                                                                                                     tAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGc                                          >  1:607223/1‑148 (MQ=255)
                                                                                                      aaCGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCg                                         <  2:865205/148‑1 (MQ=255)
                                                                                                      aaCGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTCACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCg                                         >  1:2941404/1‑148 (MQ=255)
                                                                                                          tttCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGc                                     >  1:2013139/1‑148 (MQ=255)
                                                                                                           ttCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGcc                                    >  2:2410158/1‑148 (MQ=255)
                                                                                                                  attattCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTc                             <  1:1669809/148‑1 (MQ=255)
                                                                                                                    tattCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAt                           >  1:1112250/1‑148 (MQ=255)
                                                                                                                      ttCGTACCGCAGATGCCGTTGGCTTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAtct                         >  1:745130/1‑148 (MQ=255)
                                                                                                                      ttCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAtct                         >  1:804637/1‑148 (MQ=255)
                                                                                                                       tCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTACCAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCAtctc                        >  2:169124/1‑148 (MQ=255)
                                                                                                                       tCGTACCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCAGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCTGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACAATTGGCGATGCCGTCGCTCAtcta                        >  2:169366/1‑147 (MQ=255)
                                                                                                                           aCCGCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTCCAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAg                    >  2:2319842/1‑148 (MQ=255)
                                                                                                                              gCAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGcc                 >  1:2203190/1‑148 (MQ=255)
                                                                                                                               cAGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCa                <  1:1263396/148‑1 (MQ=255)
                                                                                                                                aGATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAg               <  2:1818075/148‑1 (MQ=255)
                                                                                                                                 gATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCCgg              >  2:2549385/1‑148 (MQ=255)
                                                                                                                                 gATGCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAgg              >  2:766150/1‑148 (MQ=255)
                                                                                                                                    gCCGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCa           <  1:1156896/148‑1 (MQ=255)
                                                                                                                                      cGTTGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATg         <  2:1565858/148‑1 (MQ=255)
                                                                                                                                        ttGGCGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCa       >  1:3396112/1‑148 (MQ=255)
                                                                                                                                            cGTAAATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGAtt   <  1:3250897/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACCATGGCTTCGGCAGCGGCGGGTAGTAACAGCTGGGTACAGGTGAAAACACACCGCACCATTGGCGATGCCGTCGCTCATCTCAAAGGCCAGGGCATGCAGATTC  >  CP009273/3817866‑3818154

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