Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,923,700 G→T 19.7% R252S (CGC→AGC)  ravA ← hexameric AAA+ MoxR family ATPase, putative molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,923,7000GT19.7% 672.5 / 127.9 238R252S (CGC→AGC) ravAhexameric AAA+ MoxR family ATPase, putative molecular chaperone
Reads supporting (aligned to +/- strand):  ref base G (86/105);  new base T (25/22);  total (111/127)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

CAATACATCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAGGCATCAG  >  CP009273/3923551‑3923849
                                                                                                                                                     |                                                                                                                                                     
caATACATCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGcg                                                                                                                                                       <  1:722602/150‑1 (MQ=255)
 aaTACATCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGcga                                                                                                                                                      >  2:3357274/1‑150 (MQ=255)
  aTACATCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGcgac                                                                                                                                                     >  1:4459496/1‑150 (MQ=255)
  aTACATCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGcgac                                                                                                                                                     >  1:762573/1‑150 (MQ=255)
  aTACATCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGcgac                                                                                                                                                     >  2:3494933/1‑150 (MQ=255)
    aCATCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGctacga                                                                                                                                                   >  2:552735/1‑150 (MQ=255)
     cATCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGAt                                                                                                                                                  <  1:907524/150‑1 (MQ=255)
     cATCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGAt                                                                                                                                                  <  2:2809648/150‑1 (MQ=255)
      aTCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATc                                                                                                                                                 <  1:4356197/150‑1 (MQ=255)
      aTCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATc                                                                                                                                                 >  2:331735/1‑150 (MQ=255)
       tCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATcc                                                                                                                                                >  2:822256/1‑150 (MQ=255)
             tgttgttgtATCAAATTCAGGCTTTGCGCTTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACa                                                                                                                                          <  1:4155560/150‑1 (MQ=255)
                  ttgtATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAgg                                                                                                                                     <  2:3084845/150‑1 (MQ=255)
                  ttgtATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAgg                                                                                                                                     <  2:1121920/150‑1 (MQ=255)
                         aaaTTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATcc                                                                                                                              <  1:2543903/150‑1 (MQ=255)
                          aaTTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCg                                                                                                                             <  2:1734479/150‑1 (MQ=255)
                           aTTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCgg                                                                                                                            <  1:3259287/150‑1 (MQ=255)
                           aTTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCgg                                                                                                                            <  1:1693218/150‑1 (MQ=255)
                            ttCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGt                                                                                                                           <  1:197814/150‑1 (MQ=255)
                            ttCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGt                                                                                                                           <  2:1474586/150‑1 (MQ=255)
                               aGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAATGCGCATCCGGTAAt                                                                                                                        >  2:361268/1‑150 (MQ=255)
                                ggCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGACCCGAGACATAAGGCGCATCCGGTAAtt                                                                                                                       >  1:4446989/1‑150 (MQ=255)
                                 gCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAAttt                                                                                                                      >  2:455911/1‑150 (MQ=255)
                                 gCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAAttt                                                                                                                      >  1:506547/1‑150 (MQ=255)
                                  cTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTa                                                                                                                     >  1:3018949/1‑150 (MQ=255)
                                   tttGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTAt                                                                                                                    >  2:3336809/1‑150 (MQ=255)
                                   tttGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTAt                                                                                                                    <  2:2048663/150‑1 (MQ=255)
                                          gTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTg                                                                                                             <  2:3578198/150‑1 (MQ=255)
                                            cATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCt                                                                                                           <  2:626261/150‑1 (MQ=255)
                                             aTACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTg                                                                                                          <  1:2742218/150‑1 (MQ=255)
                                                 cacaGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGgcgc                                                                                                      >  2:2163251/1‑150 (MQ=255)
                                                   caGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCgcag                                                                                                    >  1:4564580/1‑150 (MQ=255)
                                                    aGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCgcagc                                                                                                   <  1:638099/150‑1 (MQ=255)
                                                     ggCACTCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCgcagca                                                                                                  <  2:666719/150‑1 (MQ=255)
                                                       cAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATa                                                                                                >  1:1105432/1‑150 (MQ=255)
                                                       cAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATa                                                                                                <  1:1839391/150‑1 (MQ=255)
                                                        aaTCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTAATCCAGTTGCTGGCGCAGCATaa                                                                                               >  1:2427462/1‑150 (MQ=255)
                                                          tCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATaaaa                                                                                             <  2:1374176/150‑1 (MQ=255)
                                                          tCTTTCAGCAAAATGAGATCGACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATACGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATaaaa                                                                                             <  2:2717238/150‑1 (MQ=255)
                                                           cTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATaaaaa                                                                                            <  2:4043932/150‑1 (MQ=255)
                                                           cTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATaaaaa                                                                                            <  2:4318948/150‑1 (MQ=255)
                                                           cTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATaaaaa                                                                                            <  1:2787270/150‑1 (MQ=255)
                                                              tCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCa                                                                                         >  2:2670506/1‑150 (MQ=255)
                                                               cAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAg                                                                                        <  2:3759307/150‑1 (MQ=255)
                                                               cAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAg                                                                                        >  2:3290576/1‑150 (MQ=255)
                                                                aGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGc                                                                                       >  1:4575385/1‑150 (MQ=255)
                                                                   aaaaTGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCa                                                                                    >  2:3126212/1‑150 (MQ=255)
                                                                            tCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGa                                                                           <  1:2126110/150‑1 (MQ=255)
                                                                             cTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCCGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGAt                                                                          >  2:2886476/1‑150 (MQ=255)
                                                                              tACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATc                                                                         <  2:2051891/150‑1 (MQ=255)
                                                                              tACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATc                                                                         >  1:3320196/1‑150 (MQ=255)
                                                                              tACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATc                                                                         <  1:4760628/150‑1 (MQ=255)
                                                                               aCCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCg                                                                        <  2:4795892/150‑1 (MQ=255)
                                                                                 cGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCggg                                                                      <  2:2975660/150‑1 (MQ=255)
                                                                                  ggggCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGc                                                                     >  1:4873389/1‑150 (MQ=255)
                                                                                    ggCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAg                                                                   <  2:2056841/150‑1 (MQ=255)
                                                                                      cAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCg                                                                 >  1:3909192/1‑150 (MQ=255)
                                                                                        aCAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTa                                                               >  1:3496715/1‑150 (MQ=255)
                                                                                          aGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAAt                                                             >  1:3167009/1‑150 (MQ=255)
                                                                                            cACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAAttt                                                           >  2:1279719/1‑150 (MQ=255)
                                                                                                 cGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCa                                                      <  2:2144308/150‑1 (MQ=255)
                                                                                                 cGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCa                                                      <  2:2144411/150‑1 (MQ=255)
                                                                                                 cGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCa                                                      <  2:2144572/150‑1 (MQ=255)
                                                                                                   aCCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAAt                                                    <  2:3177569/150‑1 (MQ=255)
                                                                                                   aCCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAAt                                                    <  2:2360398/150‑1 (MQ=255)
                                                                                                   aCCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAAt                                                    <  2:3651966/150‑1 (MQ=255)
                                                                                                     cGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAAtct                                                  <  1:681884/150‑1 (MQ=255)
                                                                                                     cGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAAtct                                                  <  2:4587190/150‑1 (MQ=255)
                                                                                                       cTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATctct                                                >  2:1971303/1‑150 (MQ=255)
                                                                                                       cTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATctct                                                >  1:202994/1‑150 (MQ=255)
                                                                                                        tAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCtt                                               <  1:4116523/150‑1 (MQ=255)
                                                                                                        tAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCtt                                               >  1:1321171/1‑150 (MQ=255)
                                                                                                        tAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTACAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCtt                                               >  1:566574/1‑150 (MQ=255)
                                                                                                         aaaaaaGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCttt                                              >  1:1578616/1‑150 (MQ=255)
                                                                                                         aaaaaaGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCttt                                              <  1:4059437/150‑1 (MQ=255)
                                                                                                           aaaaGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCt                                            <  1:2394621/150‑1 (MQ=255)
                                                                                                           aaaaGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCt                                            <  1:1444865/150‑1 (MQ=255)
                                                                                                             aaGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGc                                          >  2:4637237/1‑150 (MQ=255)
                                                                                                              aGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGcc                                         <  2:2823151/150‑1 (MQ=255)
                                                                                                               ggCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCa                                        <  1:1358639/150‑1 (MQ=255)
                                                                                                                  gcTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCg                                     <  1:3357435/150‑1 (MQ=255)
                                                                                                                     ggCCTGCAATAATCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTc                                  <  1:882983/150‑1 (MQ=255)
                                                                                                                        cTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCata                               <  1:3660298/150‑1 (MQ=255)
                                                                                                                            aataatCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATAttct                           >  1:4022494/1‑150 (MQ=255)
                                                                                                                             ataatCGGATCGCTTTTTTCCAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTGTCTGCCAGCGTTCATAttctt                          >  2:2558302/1‑150 (MQ=255)
                                                                                                                             ataatCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATAttctt                          >  1:1075060/1‑150 (MQ=255)
                                                                                                                                  cGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCt                     >  1:2265987/1‑150 (MQ=255)
                                                                                                                                  cGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCt                     <  1:3494933/150‑1 (MQ=255)
                                                                                                                                     aTCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCtgtg                  <  1:3126212/150‑1 (MQ=255)
                                                                                                                                      tCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGa                 <  1:2050575/150‑1 (MQ=255)
                                                                                                                                       cGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGAc                >  1:854826/1‑150 (MQ=255)
                                                                                                                                              ttCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAg         >  2:2550647/1‑150 (MQ=255)
                                                                                                                                                ccAGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAGGc       >  1:2147901/1‑150 (MQ=255)
                                                                                                                                                 cAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAGGCa      <  1:2910122/150‑1 (MQ=255)
                                                                                                                                                 cAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAGGCa      >  1:3641223/1‑150 (MQ=255)
                                                                                                                                                  aGCTACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAGGCAt     <  1:3459331/150‑1 (MQ=255)
                                                                                                                                                  aGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAGGCAt     >  1:4341591/1‑150 (MQ=255)
                                                                                                                                                    cgacgaTCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAGGCATCa   <  2:2885030/150‑1 (MQ=255)
                                                                                                                                                     gacgaTCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAGGCATCAg  >  2:4062005/1‑150 (MQ=255)
                                                                                                                                                     gacgaTCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAGGCATCAg  >  2:1791094/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAATACATCAATTTGTTGTTGTATCAAATTCAGGCTTTGCGCGTCATACCACAGGCAATCTTTCAGCAAAATGAGATCTACCGGGGCAACAGCACTGCGACCGCTAAAAAAGGCGCTGGCCTGCAATAATCGGATCGCTTTTTTCCAGCGACGATCCGAGACATAAGGCGCATCCGGTAATTTATCCAGTTGCTGGCGCAGCATAAAAATCAGCTCAAATACATGATCGGGCAGCGTAATTTCACCAATCTCTTTCTGCCAGCGTTCATATTCTTCATCTGTGACCTGCAAGGCATCAG  >  CP009273/3923551‑3923849

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: