Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,958,599 A→T 29.8% L131Q (CTG→CAG)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,958,5990AT29.8% 406.7 / 231.7 252L131Q (CTG→CAG) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (84/93);  new base T (39/36);  total (123/129)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 255 total aligned reads displayed.

GCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCAC  >  CP009273/3958451‑3958747
                                                                                                                                                    |                                                                                                                                                    
  ggATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGc                                                                                                                                                    >  1:2353162/1‑149 (MQ=255)
   gATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCa                                                                                                                                                   >  1:3037784/1‑149 (MQ=255)
   gATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCAACCTCCAGAACCACCACCAGAAAGGCACCCAGGGTCAGCTGGTAATGCAAGCCGTCGGCAATAAAACCCCCCGAGGAGCCAAGCAGAATGTCAACGCCGCTTTCCAGCa                                                                                                                                                   >  1:3375593/1‑149 (MQ=255)
    atatCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTATTGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCAc                                                                                                                                                  <  1:3843422/149‑1 (MQ=255)
        cTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTc                                                                                                                                              <  2:3940547/149‑1 (MQ=255)
         tttAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCa                                                                                                                                             <  1:3584630/149‑1 (MQ=255)
         tttAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCa                                                                                                                                             <  2:1592566/149‑1 (MQ=255)
             ataaAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTATTGTGGTTCTGCTTGGCGTAGTAAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTg                                                                                                                                         <  1:2857659/149‑1 (MQ=255)
             ataaAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAAGGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTg                                                                                                                                         <  1:316455/149‑1 (MQ=255)
              taaAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGt                                                                                                                                        >  2:1788228/1‑149 (MQ=255)
               aaaGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGtt                                                                                                                                       <  2:2256587/149‑1 (MQ=255)
                caGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGttt                                                                                                                                      <  2:3313495/148‑1 (MQ=255)
                    ccAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCg                                                                                                                                  >  1:1910636/1‑149 (MQ=255)
                      aGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGTGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTa                                                                                                                                <  1:4194649/149‑1 (MQ=255)
                       gATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAg                                                                                                                               <  2:2054032/149‑1 (MQ=255)
                       gATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAg                                                                                                                               <  2:2964131/149‑1 (MQ=255)
                        aTCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGc                                                                                                                              >  1:4487982/1‑149 (MQ=255)
                            tACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATc                                                                                                                          >  1:1286383/1‑149 (MQ=255)
                              cATGCGATCGTCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATcac                                                                                                                        <  1:190115/149‑1 (MQ=255)
                              cATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTATTGTGGTTCTGGTTGGCGTATTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATcac                                                                                                                        <  1:2093750/149‑1 (MQ=255)
                                tGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATcacca                                                                                                                      <  2:2670951/149‑1 (MQ=255)
                                tGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATcacca                                                                                                                      >  1:596095/1‑149 (MQ=255)
                                   gATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCg                                                                                                                   >  1:1063839/1‑149 (MQ=255)
                                   gATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCg                                                                                                                   <  1:4200897/149‑1 (MQ=255)
                                        gCTTCGTCCTGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGc                                                                                                              >  2:657916/1‑149 (MQ=255)
                                           tCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAg                                                                                                           >  1:1117981/1‑149 (MQ=255)
                                              tCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGAcc                                                                                                        >  2:2108872/1‑149 (MQ=255)
                                                 aGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAg                                                                                                     <  1:2503448/149‑1 (MQ=255)
                                                  gTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGc                                                                                                    >  2:4573778/1‑149 (MQ=255)
                                                  gTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGc                                                                                                    >  2:2175653/1‑149 (MQ=255)
                                                    accaccaccTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCtt                                                                                                  <  2:1374170/149‑1 (MQ=255)
                                                       accaccTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAg                                                                                               >  1:3122935/1‑149 (MQ=255)
                                                       accaccTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAg                                                                                               <  2:909311/149‑1 (MQ=255)
                                                           ccTGAAGGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCa                                                                                           <  2:1469130/149‑1 (MQ=255)
                                                            cTGAATGGCACCGAGGTTATTGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAg                                                                                          <  1:2209177/149‑1 (MQ=255)
                                                             tGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGt                                                                                         <  1:2909490/149‑1 (MQ=255)
                                                               aaTGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAg                                                                                       <  2:4390625/149‑1 (MQ=255)
                                                                aTGGCCCCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGc                                                                                      >  2:3392226/1‑149 (MQ=255)
                                                                aTGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGc                                                                                      >  2:1003539/1‑149 (MQ=255)
                                                                 tGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTTGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCt                                                                                     <  1:501033/149‑1 (MQ=255)
                                                                 tGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCt                                                                                     <  1:1058578/149‑1 (MQ=255)
                                                                  ggCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCtt                                                                                    >  1:2992059/1‑149 (MQ=255)
                                                                   gCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTc                                                                                   >  1:736714/1‑149 (MQ=255)
                                                                       cGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGcc                                                                               >  2:2316353/1‑149 (MQ=255)
                                                                       cGAGGTTAATGTGGTTCTGCTTGGCGTAGTAAATTAAACGCCCCGTGGTGCCAATAAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGcc                                                                               <  1:1691457/149‑1 (MQ=255)
                                                                        gAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCa                                                                              >  1:4549643/1‑149 (MQ=255)
                                                                         aGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAg                                                                             >  1:1000012/1‑149 (MQ=255)
                                                                           gTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCg                                                                           <  1:1812171/149‑1 (MQ=255)
                                                                            ttAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCgg                                                                          <  1:3957083/149‑1 (MQ=255)
                                                                             tAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGt                                                                         <  2:2379877/149‑1 (MQ=255)
                                                                             tAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGt                                                                         <  2:3588584/149‑1 (MQ=255)
                                                                              aaTGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGtt                                                                        >  1:527294/1‑149 (MQ=255)
                                                                               aTGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTc                                                                       <  2:430523/149‑1 (MQ=255)
                                                                                     ttCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTc                                                                 <  2:4771939/149‑1 (MQ=255)
                                                                                     ttCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTc                                                                 >  1:1454685/1‑149 (MQ=255)
                                                                                     ttCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTc                                                                 >  2:4726964/1‑149 (MQ=255)
                                                                                        tGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATCTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGc                                                              >  2:115327/1‑149 (MQ=255)
                                                                                                  aGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCGGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGc                                                    >  1:369941/1‑149 (MQ=255)
                                                                                                    tCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCAc                                                  >  1:68400/1‑149 (MQ=255)
                                                                                                    tCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCAc                                                  >  1:2754062/1‑149 (MQ=255)
                                                                                                    tCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCAc                                                  <  2:1620717/149‑1 (MQ=255)
                                                                                                    tCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCAc                                                  <  2:4008046/149‑1 (MQ=255)
                                                                                                      aaTTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACgg                                                >  2:3192388/1‑149 (MQ=255)
                                                                                                       aTTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGc                                               >  1:177155/1‑149 (MQ=255)
                                                                                                       aTTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGc                                               >  1:4053035/1‑149 (MQ=255)
                                                                                                       aTTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGc                                               <  2:596095/149‑1 (MQ=255)
                                                                                                          aaaCGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAg                                            <  2:1278452/149‑1 (MQ=255)
                                                                                                           aaCGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGt                                           <  1:2654240/149‑1 (MQ=255)
                                                                                                            aCGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGtt                                          >  1:2269142/1‑149 (MQ=255)
                                                                                                             cGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTc                                         <  2:1006244/149‑1 (MQ=255)
                                                                                                              gCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCa                                        >  1:3775477/1‑149 (MQ=255)
                                                                                                              gCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCa                                        >  1:3311164/1‑149 (MQ=255)
                                                                                                              gCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCa                                        >  2:1855881/1‑149 (MQ=255)
                                                                                                                   gtggtgCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGt                                   <  2:3866690/149‑1 (MQ=255)
                                                                                                                     ggtgCCAATCAGAATGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCg                                 >  2:4615337/1‑149 (MQ=255)
                                                                                                                         ccAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGc                             <  1:3740270/149‑1 (MQ=255)
                                                                                                                         ccAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGc                             <  1:2803919/149‑1 (MQ=255)
                                                                                                                          cAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGcc                            <  1:3974695/149‑1 (MQ=255)
                                                                                                                          cAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGcc                            <  2:165235/149‑1 (MQ=255)
                                                                                                                             tCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATa                         >  2:1774811/1‑149 (MQ=255)
                                                                                                                               aGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAAt                       <  1:480628/149‑1 (MQ=255)
                                                                                                                                gAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAAtt                      >  2:4164995/1‑149 (MQ=255)
                                                                                                                                  aTGTCAACGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTaa                    <  2:1931281/149‑1 (MQ=255)
                                                                                                                                  aTGTCAACGCCGCGTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTaa                    <  2:3670095/149‑1 (MQ=255)
                                                                                                                                      cAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCa                <  2:4599813/149‑1 (MQ=255)
                                                                                                                                       aaCGCCGCTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCAc               >  2:3025521/1‑149 (MQ=255)
                                                                                                                                       aaCGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCAc               <  1:2186695/149‑1 (MQ=255)
                                                                                                                                          gccgcTTTCCTGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCAcgcg            <  1:1388513/149‑1 (MQ=255)
                                                                                                                                          gccgcTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCAcgcg            >  2:915943/1‑149 (MQ=255)
                                                                                                                                          gccgcTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCAcgcg            >  2:239878/1‑149 (MQ=255)
                                                                                                                                            cgcTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGc          >  1:4493946/1‑149 (MQ=255)
                                                                                                                                            cgcTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGc          <  1:93049/149‑1 (MQ=255)
                                                                                                                                            cgcTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGc          <  1:1279119/149‑1 (MQ=255)
                                                                                                                                             gcTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCt         >  1:4254596/1‑149 (MQ=255)
                                                                                                                                               tttCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGa       >  1:1477609/1‑149 (MQ=255)
                                                                                                                                               tttCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGa       <  2:2584630/149‑1 (MQ=255)
                                                                                                                                                ttCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGAt      >  1:4263166/1‑149 (MQ=255)
                                                                                                                                                  ccAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTc    >  1:4656213/1‑149 (MQ=255)
                                                                                                                                                  ccAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTc    <  2:909718/149‑1 (MQ=255)
                                                                                                                                                    aGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCAc  >  1:1596529/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGCCAATCAGAATGTCAACGCCGCTTTCCAGCACTTTCAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCAC  >  CP009273/3958451‑3958747

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: