Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,424,363 T→A 100% T138S (ACC→TCC)  tsaC ← tRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein, threonine‑dependent ADP‑forming ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,424,3630TA98.9% 607.6 / ‑2.4 184T138S (ACC→TCC) tsaCtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein, threonine‑dependent ADP‑forming ATPase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (100/82);  total (102/82)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 187 total aligned reads displayed.

TTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGG  >  CP009273/3424220‑3424507
                                                                                                                                               |                                                                                                                                                
ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGgaagaa                                                                                                                                              <  2:669247/148‑1 (MQ=255)
ttCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGgaagaa                                                                                                                                              >  1:3253000/1‑148 (MQ=255)
 tCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGaaa                                                                                                                                             >  2:2201261/1‑148 (MQ=255)
 tCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGaaa                                                                                                                                             >  2:1580533/1‑148 (MQ=255)
    cccGTCAGGTCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCa                                                                                                                                          <  2:1230812/148‑1 (MQ=255)
    cccGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCa                                                                                                                                          >  2:2503930/1‑148 (MQ=255)
    cccGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCa                                                                                                                                          <  2:2577842/148‑1 (MQ=255)
     ccGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAg                                                                                                                                         <  1:1371856/148‑1 (MQ=255)
     ccGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAg                                                                                                                                         <  1:3119467/148‑1 (MQ=255)
      cGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGc                                                                                                                                        >  1:1333097/1‑148 (MQ=255)
        tCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCgg                                                                                                                                      >  2:1138413/1‑148 (MQ=255)
        tCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCgg                                                                                                                                      <  1:204237/148‑1 (MQ=255)
         cAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGt                                                                                                                                     >  1:3215147/1‑148 (MQ=255)
          aGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGtt                                                                                                                                    <  1:3468650/148‑1 (MQ=255)
              cagcgcgGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTAcc                                                                                                                                <  2:181924/145‑1 (MQ=255)
                 cgcgGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATa                                                                                                                             <  1:2143250/148‑1 (MQ=255)
                 cgcgGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATa                                                                                                                             >  1:42560/1‑148 (MQ=255)
                    ggATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGc                                                                                                                          <  1:3498938/148‑1 (MQ=255)
                     gATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGcc                                                                                                                         <  1:558907/148‑1 (MQ=255)
                     gATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGcc                                                                                                                         >  1:2939128/1‑148 (MQ=255)
                      aTTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCt                                                                                                                        >  1:1515465/1‑148 (MQ=255)
                      aTTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTAGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCt                                                                                                                        >  1:1515636/1‑148 (MQ=255)
                         tCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGc                                                                                                                     >  2:2493811/1‑148 (MQ=255)
                         tCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGc                                                                                                                     <  2:1485448/148‑1 (MQ=255)
                         tCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGc                                                                                                                     <  1:1166354/148‑1 (MQ=255)
                            gAATGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAATTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGcaca                                                                                                                  <  1:3389424/148‑1 (MQ=255)
                            gAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGcaca                                                                                                                  <  1:2460227/148‑1 (MQ=255)
                             acggATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                 <  1:642073/146‑1 (MQ=255)
                             aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTTTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                 >  2:2826326/1‑148 (MQ=255)
                             aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                 <  2:2726321/148‑1 (MQ=255)
                             aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                 >  2:667912/1‑148 (MQ=255)
                             aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                 <  2:499994/148‑1 (MQ=255)
                             aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACaa                                                                                                                 <  1:2265159/148‑1 (MQ=255)
                             aaGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGATGGCACTGGAAGAAACCAGCGGTTTACCATAAGACTGGCACaa                                                                                                                 >  2:1075139/1‑148 (MQ=255)
                                 aTTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCa                                                                                                             >  1:3005703/1‑148 (MQ=255)
                                  tttAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCaa                                                                                                            <  1:2334956/148‑1 (MQ=255)
                                     aaaCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAacca                                                                                                         >  1:2936265/1‑148 (MQ=255)
                                       aCGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAaccacc                                                                                                       >  1:1209314/1‑148 (MQ=255)
                                          cccccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAAc                                                                                                    >  2:185197/1‑148 (MQ=255)
                                           ccccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACg                                                                                                   >  2:3395703/1‑148 (MQ=255)
                                            cccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCACCgg                                                                                                  >  2:447826/1‑148 (MQ=255)
                                            cccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACgg                                                                                                  >  2:3097690/1‑148 (MQ=255)
                                            cccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACgg                                                                                                  <  2:2331725/148‑1 (MQ=255)
                                             ccccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGa                                                                                                 <  2:2989574/148‑1 (MQ=255)
                                             ccacGTTTCACCAGGCAAAACCGGGAACGCCGCTCCAAATTGTGCGCGAACTTCGTCTAATGTTCGAAAAGGTGGCAATCCACTCAAGTTTGCACTGGAAGAAACCAGCGGTTTACCATAAGCCAGTCACAAAGAAACCACCAACGGa                                                                                                 >  2:124122/4‑148 (MQ=255)
                                               ccGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATg                                                                                               <  2:943227/148‑1 (MQ=255)
                                                cGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATgg                                                                                              <  2:2222869/148‑1 (MQ=255)
                                                 gTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGt                                                                                             >  1:313474/1‑148 (MQ=255)
                                                   ttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCg                                                                                           >  1:1177537/1‑148 (MQ=255)
                                                   ttCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCg                                                                                           >  1:3136211/1‑148 (MQ=255)
                                                     cACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGt                                                                                         >  1:591450/1‑148 (MQ=255)
                                                      aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                                                                                        >  2:2678985/1‑148 (MQ=255)
                                                      aCCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTg                                                                                        >  2:1398808/1‑148 (MQ=255)
                                                       ccAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGa                                                                                       >  2:862856/1‑148 (MQ=255)
                                                       ccAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGAGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGa                                                                                       >  1:286202/1‑148 (MQ=255)
                                                           gCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCg                                                                                   >  2:1634025/1‑148 (MQ=255)
                                                                 ccGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGc                                                                             >  2:681375/1‑148 (MQ=255)
                                                                   gggAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCaa                                                                           >  1:1249419/1‑148 (MQ=255)
                                                                     gAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGc                                                                         >  1:402487/1‑148 (MQ=255)
                                                                        cgccgcGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGaa                                                                      <  2:3041273/148‑1 (MQ=255)
                                                                        cgccgcGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGaa                                                                      <  2:1008350/148‑1 (MQ=255)
                                                                           cgcgCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCa                                                                   >  2:1454216/1‑148 (MQ=255)
                                                                             cgcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCaaa                                                                 >  1:833769/1‑148 (MQ=255)
                                                                              gcCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAg                                                                <  1:300484/148‑1 (MQ=255)
                                                                               ccAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGc                                                               <  2:1406637/148‑1 (MQ=255)
                                                                               ccAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGc                                                               >  2:3231683/1‑148 (MQ=255)
                                                                                   aTTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGcc                                                           >  1:1685906/1‑148 (MQ=255)
                                                                                         cgcgAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCaa                                                     >  2:2313059/1‑148 (MQ=255)
                                                                                           cgAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAAcc                                                   >  1:2335927/1‑148 (MQ=255)
                                                                                              aCTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAgc                                                >  2:2905578/1‑148 (MQ=255)
                                                                                               cTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAgcg                                               >  1:3428985/1‑148 (MQ=255)
                                                                                                    tCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtg                                          >  1:143122/1‑148 (MQ=255)
                                                                                                    tCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtg                                          >  2:1305480/1‑148 (MQ=255)
                                                                                                     cTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtgt                                         >  2:625286/1‑148 (MQ=255)
                                                                                                     cTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGgtgt                                         >  2:3137605/1‑148 (MQ=255)
                                                                                                        cTGTTCGACATGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttgt                                      <  1:2202123/148‑1 (MQ=255)
                                                                                                        cTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGtgttgt                                      >  2:1256032/1‑148 (MQ=255)
                                                                                                           ttCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGc                                   >  1:442640/1‑148 (MQ=255)
                                                                                                           ttCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGc                                   >  2:2607054/1‑148 (MQ=255)
                                                                                                                 aaGGTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGgcgc                             <  2:3452641/148‑1 (MQ=255)
                                                                                                                    gTGGCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCggg                          >  2:1982826/1‑148 (MQ=255)
                                                                                                                       gCAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGaaa                       <  2:1690278/148‑1 (MQ=255)
                                                                                                                       gCAATCCAATAAAGTTGGCACTGGAAGAAACCTGCGGTTTACAATAAGCCTGGCAAAAAGAAACCACCACGGGATGGTCGGTTACTCGTACACCAAGCGAATAAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGaaa                       <  1:1861568/148‑1 (MQ=255)
                                                                                                                        cAATCCACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAg                      >  1:2420668/1‑148 (MQ=255)
                                                                                                                            ccACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACaa                  >  1:472766/1‑148 (MQ=255)
                                                                                                                             cACTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACaaa                 <  2:1333097/148‑1 (MQ=255)
                                                                                                                               cTCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAgg               >  1:2713362/1‑148 (MQ=255)
                                                                                                                                tCAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGt              >  2:813058/1‑148 (MQ=255)
                                                                                                                                 cAAGTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTg             >  2:1177725/1‑148 (MQ=255)
                                                                                                                                    gTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACa          >  2:769699/1‑148 (MQ=255)
                                                                                                                                    gTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACa          >  2:2156400/1‑148 (MQ=255)
                                                                                                                                    gTTGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACa          >  2:387039/1‑148 (MQ=255)
                                                                                                                                      tGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAgg        <  1:2547319/148‑1 (MQ=255)
                                                                                                                                      tGGCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAgg        <  1:3090088/148‑1 (MQ=255)
                                                                                                                                        gCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCACCGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAc      <  1:447826/148‑1 (MQ=255)
                                                                                                                                        gCACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAc      <  2:3215095/148‑1 (MQ=255)
                                                                                                                                         cACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGACAGACAAAGGTGACAGGAcc     >  2:2272313/1‑148 (MQ=255)
                                                                                                                                         cACTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGAcc     <  2:560146/148‑1 (MQ=255)
                                                                                                                                          aCTGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCt    <  2:563015/148‑1 (MQ=255)
                                                                                                                                            tGGAAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTgg  <  2:2335927/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TTCACCCGTCAGGGCATCGCGGATTTCTGAAGGATTTAAACGCCCCCCCGTTTCACCAGGCACAACCGGGAACGCCGCGCCAAATTGTGCGCGAACTTCGTCTACTGTTCGACAAGGTGGCAATCCACTCAAGTTGGCACTGGTAGAAACCAGCGGTTTACCATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGG  >  CP009273/3424220‑3424507

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: