Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,152,819 T→G 6.0% K127T (AAG→ACG)  trmA ← tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,152,8190TG 6.0% 646.6 / 29.1 227K127T (AAG→ACG) trmAtRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (109/104);  new base G (3/11);  total (112/115)
Fisher's exact test for biased strand distribution p-value = 5.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

CTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATAAGTTCACTGGCGGCGG  >  CP009273/4152673‑4152966
                                                                                                                                                  |                                                                                                                                                   
cTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATAAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGcta                                                                                                                                                    <  1:2867026/148‑2 (MQ=255)
  gCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGcttct                                                                                                                                                  >  1:3554621/1‑148 (MQ=255)
  gCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGcttct                                                                                                                                                  >  1:3554117/1‑148 (MQ=255)
    gACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCGTCTTa                                                                                                                                                <  1:1015063/148‑1 (MQ=255)
     aCCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTAt                                                                                                                                               >  2:81141/1‑148 (MQ=255)
     aCCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTAt                                                                                                                                               <  2:3149915/148‑1 (MQ=255)
         gCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCGTCTTATGGta                                                                                                                                           <  1:527519/148‑1 (MQ=255)
         gCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCGTCTTATGGta                                                                                                                                           <  1:2302241/148‑1 (MQ=255)
          cAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGtat                                                                                                                                          >  2:1863278/1‑148 (MQ=255)
          cAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGtat                                                                                                                                          >  2:1198393/1‑148 (MQ=255)
          cAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGtat                                                                                                                                          <  2:2823938/148‑1 (MQ=255)
          cAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGtat                                                                                                                                          <  2:49367/148‑1 (MQ=255)
           aGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGtata                                                                                                                                         <  2:119293/148‑1 (MQ=255)
            gACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAg                                                                                                                                        >  2:225674/1‑148 (MQ=255)
            gACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAg                                                                                                                                        >  2:2156138/1‑148 (MQ=255)
            gACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAg                                                                                                                                        <  2:1774068/148‑1 (MQ=255)
             aCGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGc                                                                                                                                       >  1:1193220/1‑148 (MQ=255)
               gTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAg                                                                                                                                     >  2:760210/1‑148 (MQ=255)
               gTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAg                                                                                                                                     <  1:1627620/148‑1 (MQ=255)
                ttCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAgg                                                                                                                                    >  2:814136/1‑148 (MQ=255)
                 tCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAggg                                                                                                                                   >  2:1422181/1‑148 (MQ=255)
                 tCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCGTCTTATGGTATAGCAggg                                                                                                                                   <  2:1985303/148‑1 (MQ=255)
                 tCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCGTCTTATGGTATAGCAggg                                                                                                                                   <  2:165911/148‑1 (MQ=255)
                     cGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAac                                                                                                                               <  2:1333388/148‑1 (MQ=255)
                      gATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAacc                                                                                                                              <  2:2173430/148‑1 (MQ=255)
                      gATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGGTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAacc                                                                                                                              <  2:1386264/148‑1 (MQ=255)
                       aTGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAacca                                                                                                                             <  1:2014173/148‑1 (MQ=255)
                        tGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAaccac                                                                                                                            >  1:827699/1‑148 (MQ=255)
                           aaTCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGc                                                                                                                         >  2:2930566/1‑148 (MQ=255)
                            aTCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGcc                                                                                                                        >  2:3623714/1‑148 (MQ=255)
                                tGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGAt                                                                                                                    <  2:1156078/148‑1 (MQ=255)
                                tGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGAt                                                                                                                    <  2:2063733/148‑1 (MQ=255)
                                 gATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGAtt                                                                                                                   >  1:11120/1‑148 (MQ=255)
                                    ccAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACt                                                                                                                <  1:412527/148‑1 (MQ=255)
                                     cAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTc                                                                                                               <  2:357896/148‑1 (MQ=255)
                                      aGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCa                                                                                                              <  2:3688466/148‑1 (MQ=255)
                                       gCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAg                                                                                                             >  1:2095997/1‑148 (MQ=255)
                                        cTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAgt                                                                                                            <  1:3347087/148‑1 (MQ=255)
                                        cTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAgg                                                                                                            >  2:2483574/1‑147 (MQ=255)
                                         tCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAgtg                                                                                                           <  1:1570218/148‑1 (MQ=255)
                                         tCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCGTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAgtg                                                                                                           >  2:836229/1‑148 (MQ=255)
                                          cGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAgtgt                                                                                                          <  1:1729069/148‑1 (MQ=255)
                                           gATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTa                                                                                                         <  1:1029608/148‑1 (MQ=255)
                                             ttttGGTGTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGt                                                                                                       <  2:2707626/148‑1 (MQ=255)
                                               ttGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGa                                                                                                     <  2:3543927/148‑1 (MQ=255)
                                                 ggTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAgg                                                                                                   <  2:2292705/148‑1 (MQ=255)
                                                     ttCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTaat                                                                                               >  2:1688551/1‑148 (MQ=255)
                                                         ttGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTaatcaat                                                                                           >  1:2351214/1‑148 (MQ=255)
                                                         ttGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTaatcaat                                                                                           <  2:1261156/148‑1 (MQ=255)
                                                              cGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGa                                                                                      >  2:14299/1‑148 (MQ=255)
                                                                aCCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGaac                                                                                    >  2:442509/1‑148 (MQ=255)
                                                                 ccAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGaaca                                                                                   <  1:1440907/148‑1 (MQ=255)
                                                                  cAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGaacaa                                                                                  >  2:2313815/1‑148 (MQ=255)
                                                                  cAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGaacaa                                                                                  <  2:3242772/148‑1 (MQ=255)
                                                                   aaTCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGaacaac                                                                                 >  1:28181/1‑148 (MQ=255)
                                                                      cAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTtg                                                                              >  2:2989417/1‑148 (MQ=255)
                                                                         aTGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTgg                                                                           >  1:583474/1‑148 (MQ=255)
                                                                                ttCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGaa                                                                    >  1:1682526/1‑148 (MQ=255)
                                                                                 tCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAAc                                                                   >  1:3282049/1‑148 (MQ=255)
                                                                                  cAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACg                                                                  <  1:1043313/148‑1 (MQ=255)
                                                                                  cAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACg                                                                  >  1:1819783/1‑148 (MQ=255)
                                                                                     aTTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGa                                                               <  1:3795860/148‑1 (MQ=255)
                                                                                     aTTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGa                                                               >  1:1880969/1‑148 (MQ=255)
                                                                                     aTTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGa                                                               >  1:2694172/1‑148 (MQ=255)
                                                                                     aTTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGa                                                               <  1:760210/148‑1 (MQ=255)
                                                                                     aTTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGa                                                               >  1:3693464/1‑148 (MQ=255)
                                                                                      ttCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGat                                                              >  2:3463865/1‑148 (MQ=255)
                                                                                      ttCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGat                                                              <  2:466145/148‑1 (MQ=255)
                                                                                       tCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGatt                                                             >  1:3589230/1‑148 (MQ=255)
                                                                                          gcgcgcgcAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGattatt                                                          >  1:966636/1‑148 (MQ=255)
                                                                                           cgcgcgcAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGAttatta                                                         <  1:405042/148‑1 (MQ=255)
                                                                                                 cAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGcaca                                                   >  2:3472023/1‑148 (MQ=255)
                                                                                                   gTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACAcc                                                 <  1:1066823/148‑1 (MQ=255)
                                                                                                    tGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACAccc                                                >  1:685817/1‑148 (MQ=255)
                                                                                                    tGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACAccc                                                >  1:1954923/1‑148 (MQ=255)
                                                                                                     gCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCg                                               >  2:1679294/1‑148 (MQ=255)
                                                                                                       aTCGCGCAGGGCCTCCGCTTCCTGACGCCACTCCTCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCa                                             <  2:950976/148‑1 (MQ=255)
                                                                                                          gcgcAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAatc                                          <  1:1376699/148‑1 (MQ=255)
                                                                                                               gggCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGc                                     >  1:2580448/1‑148 (MQ=255)
                                                                                                                ggCCTCCGCTTCCTGACGCCACTCATCATCCAGCGTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGcc                                    >  2:2775967/1‑148 (MQ=255)
                                                                                                                  ccTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGt                                  <  1:193236/148‑1 (MQ=255)
                                                                                                                    tCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGtca                                <  1:509728/148‑1 (MQ=255)
                                                                                                                     ccGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGtcat                               >  1:67663/1‑148 (MQ=255)
                                                                                                                       gCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGtcatca                             <  2:3554621/148‑1 (MQ=255)
                                                                                                                       gCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGtcatca                             <  2:3554117/148‑1 (MQ=255)
                                                                                                                          tCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACt                          <  1:3819420/148‑1 (MQ=255)
                                                                                                                            cTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTgg                        <  1:3222027/148‑1 (MQ=255)
                                                                                                                                cGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGa                    >  1:3242108/1‑148 (MQ=255)
                                                                                                                                 gCCACTCATAATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGAt                   <  1:2260643/148‑1 (MQ=255)
                                                                                                                                   cACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATaa                 <  1:2490988/148‑1 (MQ=255)
                                                                                                                                    aCTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATAAg                >  1:207027/1‑148 (MQ=255)
                                                                                                                                     cTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCTCCGTCCTCAACTGGTTGATAAGt               <  2:1055282/148‑1 (MQ=255)
                                                                                                                                     cTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATAAGt               >  1:2772818/1‑148 (MQ=255)
                                                                                                                                     cTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATAAGt               >  1:1260258/1‑148 (MQ=255)
                                                                                                                                      tcatcaTCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATAAGtt              <  1:2639290/148‑1 (MQ=255)
                                                                                                                                      tcatcaTCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATAAGtt              <  2:2351214/148‑1 (MQ=255)
                                                                                                                                       catcatCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATAAGTTc             >  1:3145208/1‑148 (MQ=255)
                                                                                                                                             ccAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCACTCATCGCCGTCATCAACTGGTTGATAATTTCACTggg       <  2:2502787/148‑2 (MQ=255)
                                                                                                                                                gCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATAAGTTCACTggcggc    >  2:617097/1‑148 (MQ=255)
                                                                                                                                                 cttcttATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATAAGTTCACTGgcggcg   <  2:3098346/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGGGAAACCACCGCCTGATTACTCAGTGTAGTGAGGTAATCAATCTGGAACAACTTGTGGCGCAGAACGGGATTATTACGCACACCCGCAATCATCGCCGTCATCAACTGGTTGATAAGTTCACTGGCGGCGG  >  CP009273/4152673‑4152966

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: