Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 625,165 G→A 8.5% G137E (GGA→GAA)  entH → enterobactin synthesis proofreading thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273625,1650GA 8.5% 978.8 / 58.5 282G137E (GGA→GAA) entHenterobactin synthesis proofreading thioesterase
Reads supporting (aligned to +/- strand):  ref base G (101/157);  new base A (8/16);  total (109/173)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

AACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGTTATAGAAACAAAAT  >  CP009273/625016‑625313
                                                                                                                                                     |                                                                                                                                                    
aaCACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTggg                                                                                                                                                      <  2:3400188/150‑1 (MQ=255)
  cacacCATCGCCCGGTGTCTGAGGGAAAGTTACGCGGCGTCTGCCAGCCGCTGCATCTTGTTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAAAAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGAAt                                                                                                                                                    <  1:1384595/150‑1 (MQ=255)
  cacacCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGAt                                                                                                                                                    >  2:1923375/1‑150 (MQ=255)
    cacCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGa                                                                                                                                                  <  1:3301484/150‑1 (MQ=255)
    cacCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGa                                                                                                                                                  <  2:167808/150‑1 (MQ=255)
     acCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGaa                                                                                                                                                 >  2:4563913/1‑150 (MQ=255)
     acCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGAATGaa                                                                                                                                                 <  1:933686/150‑1 (MQ=255)
     acCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGAATGaa                                                                                                                                                 <  1:4523798/150‑1 (MQ=255)
        aTCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAAccc                                                                                                                                              <  1:4510108/150‑1 (MQ=255)
        aTCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGAATGAAccc                                                                                                                                              <  2:5986/150‑1 (MQ=255)
                gtgtCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAAt                                                                                                                                      >  2:4865422/1‑150 (MQ=255)
                gtgtCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAAt                                                                                                                                      <  1:3266623/150‑1 (MQ=255)
                gtgtCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAAt                                                                                                                                      <  2:3576880/150‑1 (MQ=255)
                gtgtCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGAATGAACCCGACAGAAt                                                                                                                                      >  2:1197281/1‑150 (MQ=255)
                  gtCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTa                                                                                                                                    >  1:356084/1‑150 (MQ=255)
                   tCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAg                                                                                                                                   >  1:1684526/1‑150 (MQ=255)
                      gAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATg                                                                                                                                <  1:3779449/150‑1 (MQ=255)
                          gAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGAt                                                                                                                            <  2:3659044/150‑1 (MQ=255)
                           aaaGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGAtt                                                                                                                           <  1:1822739/150‑1 (MQ=255)
                           aaaGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGAtt                                                                                                                           <  2:3957828/150‑1 (MQ=255)
                            aaGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTg                                                                                                                          <  2:524823/150‑1 (MQ=255)
                              ggTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCa                                                                                                                        <  1:2332165/150‑1 (MQ=255)
                                 aCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGa                                                                                                                     >  1:3509210/1‑150 (MQ=255)
                                    cggcgTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGaaaa                                                                                                                  <  1:5141641/150‑1 (MQ=255)
                                        gTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCg                                                                                                              <  2:238688/150‑1 (MQ=255)
                                        gTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCg                                                                                                              >  2:1406468/1‑150 (MQ=255)
                                         tCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCgg                                                                                                             <  1:1019900/150‑1 (MQ=255)
                                          cTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGt                                                                                                            >  2:4727898/1‑150 (MQ=255)
                                           tGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGtt                                                                                                           <  2:24565/150‑1 (MQ=255)
                                           tGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGtt                                                                                                           <  2:4235564/150‑1 (MQ=255)
                                            gCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCTGGAAAACTCGGTTa                                                                                                          <  1:431230/150‑1 (MQ=255)
                                              cAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAAc                                                                                                        <  1:5069444/150‑1 (MQ=255)
                                                  cGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAg                                                                                                    <  1:3977646/150‑1 (MQ=255)
                                                  cGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGAATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAg                                                                                                    >  1:3131354/1‑150 (MQ=255)
                                                      gCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGAt                                                                                                <  2:2694183/150‑1 (MQ=255)
                                                       cATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGACTTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATc                                                                                               >  2:2816196/1‑150 (MQ=255)
                                                         tCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGa                                                                                             <  2:4628369/150‑1 (MQ=255)
                                                           ttGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGt                                                                                           <  2:2856029/150‑1 (MQ=255)
                                                           ttGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGt                                                                                           <  1:1252832/150‑1 (MQ=255)
                                                            tGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGtt                                                                                          >  1:4874933/1‑150 (MQ=255)
                                                               tCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAAc                                                                                       <  2:3768373/150‑1 (MQ=255)
                                                                cGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACa                                                                                      <  1:2854583/150‑1 (MQ=255)
                                                                cGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACa                                                                                      >  1:1520317/1‑150 (MQ=255)
                                                                cGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACa                                                                                      <  1:2078066/150‑1 (MQ=255)
                                                                 ggCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACAt                                                                                     >  2:1851313/1‑150 (MQ=255)
                                                                         cAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTgttaagt                                                                             <  2:1024482/150‑1 (MQ=255)
                                                                          agagCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTgttaagtt                                                                            <  1:77305/150‑1 (MQ=255)
                                                                           gagCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTgttaagtta                                                                           >  1:4111678/1‑150 (MQ=255)
                                                                                gggAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAatat                                                                      <  1:3202631/150‑1 (MQ=255)
                                                                                 ggAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATAtt                                                                     <  1:3433452/150‑1 (MQ=255)
                                                                                 ggAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATAtt                                                                     >  2:1108094/1‑150 (MQ=255)
                                                                                 ggAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGAATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATAtt                                                                     <  2:4387343/150‑1 (MQ=255)
                                                                                  gAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTg                                                                    <  2:3162805/150‑1 (MQ=255)
                                                                                  gAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTg                                                                    >  1:2924325/1‑150 (MQ=255)
                                                                                   ataTCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTgg                                                                   <  2:830600/148‑1 (MQ=255)
                                                                                    aaTCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGt                                                                  <  1:2166064/150‑1 (MQ=255)
                                                                                    aaTCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGt                                                                  <  2:4959373/150‑1 (MQ=255)
                                                                                     aTCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGtt                                                                 <  2:4512526/150‑1 (MQ=255)
                                                                                          cgtTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAAc                                                            <  1:3710904/150‑1 (MQ=255)
                                                                                          cgtTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAAc                                                            >  1:4333405/1‑150 (MQ=255)
                                                                                          cgtTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAAc                                                            >  2:2216876/1‑150 (MQ=255)
                                                                                             tttCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTcc                                                         >  2:3131947/1‑150 (MQ=255)
                                                                                              ttCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCg                                                        <  2:4845808/150‑1 (MQ=255)
                                                                                               tCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGTCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTGAAATATTGTTTTCAACTCCGa                                                       >  2:1656085/1‑150 (MQ=255)
                                                                                                cGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGAt                                                      >  1:4039070/1‑150 (MQ=255)
                                                                                                 gATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGAtt                                                     <  2:4186144/150‑1 (MQ=255)
                                                                                                 gATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGAtt                                                     >  2:2089033/1‑150 (MQ=255)
                                                                                                 gATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGAtt                                                     <  2:4740624/150‑1 (MQ=255)
                                                                                                     aaCAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACa                                                 >  1:1113723/1‑150 (MQ=255)
                                                                                                     aaCAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACa                                                 <  1:2272731/150‑1 (MQ=255)
                                                                                                       cAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATg                                               <  2:3664979/150‑1 (MQ=255)
                                                                                                         ggggCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGt                                             <  2:3006703/150‑1 (MQ=255)
                                                                                                         ggggCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTATGGAATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGt                                             <  1:3086185/150‑1 (MQ=255)
                                                                                                              ggcgTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCtg                                        >  1:2537450/1‑150 (MQ=255)
                                                                                                                 gTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTgtgt                                     >  1:4630495/1‑150 (MQ=255)
                                                                                                                  ttGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGtgtt                                    <  2:1586645/150‑1 (MQ=255)
                                                                                                                    gctgcACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGtgttgt                                  >  2:2778657/1‑150 (MQ=255)
                                                                                                                    gctgcACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGtgttgt                                  <  1:2414808/150‑1 (MQ=255)
                                                                                                                     ctgcACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTgttgtt                                 <  2:4051456/150‑1 (MQ=255)
                                                                                                                      tgcACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTa                                >  2:2886811/1‑150 (MQ=255)
                                                                                                                      tgcACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTa                                >  2:552270/1‑150 (MQ=255)
                                                                                                                        cACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTaaa                              <  1:106176/150‑1 (MQ=255)
                                                                                                                        cACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTaaa                              <  1:547240/150‑1 (MQ=255)
                                                                                                                         aCTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAAt                             >  1:4149570/1‑150 (MQ=255)
                                                                                                                              tCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATCTTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTAc                        <  2:3865009/150‑1 (MQ=255)
                                                                                                                              tCGGCTGGGTACGGCAGTTTTGGAATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTAc                        <  1:1005493/150‑1 (MQ=255)
                                                                                                                              tCGGCTGGGTACGGCAGTTTTGGAATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTAc                        <  1:4246291/150‑1 (MQ=255)
                                                                                                                                 gCTGGGTACGGCAGTTTTGGAATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACaaa                     <  1:738002/150‑1 (MQ=255)
                                                                                                                                     ggTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGTAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATg                 <  1:1621017/150‑1 (MQ=255)
                                                                                                                                      gTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACACTGATGt                <  2:1721997/150‑1 (MQ=255)
                                                                                                                                      gTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGt                <  1:3618974/150‑1 (MQ=255)
                                                                                                                                       tACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGtt               >  2:3220030/1‑150 (MQ=255)
                                                                                                                                         cGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGTtat             <  1:4783907/150‑1 (MQ=255)
                                                                                                                                          ggCAGTTTTGGGATGAACCCGACAGAATTAGATTAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGTtata            >  1:4978431/1‑150 (MQ=255)
                                                                                                                                          ggCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGTtata            >  1:1945782/1‑150 (MQ=255)
                                                                                                                                             aGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGTTATAGaa         >  1:388324/1‑150 (MQ=255)
                                                                                                                                             aGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGTTATAGaa         <  1:4665275/150‑1 (MQ=255)
                                                                                                                                               ttttGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGTTATAGAAAc       <  2:1492021/150‑1 (MQ=255)
                                                                                                                                                  tGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGTTATAGAAACaaa    <  1:2216876/150‑1 (MQ=255)
                                                                                                                                                   gggATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGTTATAGAAACaaaa   <  1:481328/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCAGTTTTGGGATGAACCCGACAGAATTAGATGAGATTGCAGGAAAACTCGGTTAACGGAGTGATCGAGTTAACATTGTTAAGTTAAATATTGGTTTCAACTCCGATTTACATGGTTGCTGTGTTGTTAAATTGTACAAAGATGTTATAGAAACAAAAT  >  CP009273/625016‑625313

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: