Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,815,180 T→C 8.7% E21G (GAA→GGA)  dnaE2 ← DNA‑directed DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,815,1800TC 8.7% 177.5 / 8.0 69E21G (GAA→GGA) dnaE2DNA‑directed DNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (35/28);  new base C (3/3);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

TGTCATACCAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGCTGTAGCAAGTTGCTCAAATTCATCCCAGCTGTCCA  >  CP000730/1815047‑1815320
                                                                                                                                     |                                                                                                                                             
tGTCATACCAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTAtttt                                                                                                                                        <  1:717231/141‑1 (MQ=255)
   cATACCAAAAATCGGTTTAATGTTATTTGCTATACGTGCATCATAAGATTTAGGAAAACCATACAATACATTGGTGTCAGTCAGGAGAAGTGCCTCAACATTTTCAGAAACAGCAAGTGTTACGGCATCTTCGATTTTTaa                                                                                                                                     >  1:447609/1‑141 (MQ=255)
   cATACCAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAACATTTAGGAAAACCATACAATACATTGCTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCACGTCTTCCGGCAccttccctccctac                                                                                                                                     >  1:875087/1‑132 (MQ=255)
      aCCAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCCAGTAATACGTCAGCTTCTATTTTTAAGCt                                                                                                                                  <  1:457237/141‑1 (MQ=255)
      aCCAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCACCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCt                                                                                                                                  <  1:815275/141‑1 (MQ=255)
       ccAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCtt                                                                                                                                 <  1:406274/141‑1 (MQ=255)
        cAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCCTCTTCTATTTTTAAGCTTg                                                                                                                                >  1:759368/1‑141 (MQ=255)
        cAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTg                                                                                                                                >  2:917267/1‑141 (MQ=255)
        cAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTg                                                                                                                                >  2:678945/1‑141 (MQ=255)
        cAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTg                                                                                                                                >  1:40676/1‑141 (MQ=255)
        cAAAAATCGGTTTAATGTTATTTGATATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAGGTCTTAAGGCATATTCTTTTTTTAAGCTTg                                                                                                                                >  1:609839/1‑141 (MQ=255)
           aaaTCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAAt                                                                                                                             <  1:29937/141‑1 (MQ=255)
                ggTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTACATTTata                                                                                                                        >  2:816596/1‑139 (MQ=255)
                 gTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGGATCTTCTATTTTTAAGCTTGAATTTaaaa                                                                                                                       >  2:378182/1‑139 (MQ=255)
                 gTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTCCTATTTTTAAGCTTGAATTTAACa                                                                                                                       >  2:693603/1‑141 (MQ=255)
                 gTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTCCTATTTTTAAGCTTGAATTTAACa                                                                                                                       >  2:600323/1‑141 (MQ=255)
                 gTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTCCTATTTTTAAGCTTGAATTTAACa                                                                                                                       >  2:263542/1‑141 (MQ=255)
                  tttAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACaa                                                                                                                      >  1:786676/1‑141 (MQ=255)
                   ttAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACaaa                                                                                                                     <  1:716587/141‑1 (MQ=255)
                    tAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCTGTCACTGATCATTTCAGGCCATCTTCTGTTTTTAAGCTTGATGTTAACAAAt                                                                                                                    <  1:690308/141‑1 (MQ=255)
                    tAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTACGCTTGAATTTAACAAAt                                                                                                                    <  1:501788/141‑1 (MQ=255)
                    tAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAAt                                                                                                                    >  1:335000/1‑141 (MQ=255)
                       tGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCAt                                                                                                                >  1:881668/1‑141 (MQ=255)
                        gttagggtcTATACATGCATCATGAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTTAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATa                                                                                                               >  2:225987/9‑141 (MQ=255)
                              tGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGt                                                                                                         <  1:170867/141‑1 (MQ=255)
                              tGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGt                                                                                                         <  1:163085/141‑1 (MQ=255)
                              tGCGATACATGCATCATAAAATTTAGGAAAACCATACCATACATGGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGt                                                                                                         <  2:5405/141‑1 (MQ=255)
                                    aCATGCTTCATAAAATTTAGGACACCCATACACTACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAAGAAATCATAAGCCGTATGAat                                                                                                   <  1:574487/141‑1 (MQ=255)
                                        gCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATAttt                                                                                               >  2:555082/1‑141 (MQ=255)
                                        gCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATAttt                                                                                               >  1:699943/1‑141 (MQ=255)
                                           tcatAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTaaa                                                                                            <  1:554946/141‑1 (MQ=255)
                                           tcatAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTaaa                                                                                            <  1:602767/141‑1 (MQ=255)
                                            catAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTATCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCATATGAATATTTAAat                                                                                           >  2:883892/1‑141 (MQ=255)
                                              tAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCGTCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAatat                                                                                         >  1:625025/1‑141 (MQ=255)
                                               aaaaTTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTGAACAAATCATAAGCCGTATGAATATTTAAATATg                                                                                        >  2:719064/1‑141 (MQ=255)
                                               aaaaTTTAGGAAAACCATAAAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATg                                                                                        >  1:564502/1‑141 (MQ=255)
                                                     tAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGccacca                                                                                  <  1:816596/141‑1 (MQ=255)
                                                       ggAAAACCGTACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTATCAGACACAGCAAGTCTTACGGCATCTTCTATTTTGAAGCTTGAATTTAACAAATCATATGCCGTATGAATATTTAAATATGCCACCATg                                                                                >  2:694136/1‑141 (MQ=255)
                                                       ggAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCAta                                                                                >  2:635213/1‑140 (MQ=255)
                                                       ggAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATg                                                                                >  1:160656/1‑141 (MQ=255)
                                                       ggAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAACTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATg                                                                                <  2:519157/141‑1 (MQ=255)
                                                            aCCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGa                                                                           <  2:33154/141‑1 (MQ=255)
                                                               atacaatacaTTGGTGTCAGTTATGGCACGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAATTCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATg                                                                        <  1:535427/141‑1 (MQ=255)
                                                                 acaatacaTTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGc                                                                      >  2:539210/1‑141 (MQ=255)
                                                                 acaatacaTTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGc                                                                      >  1:439488/1‑141 (MQ=255)
                                                                  caatacaTTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGACACCATGATTGAATGGcc                                                                     <  2:865297/141‑1 (MQ=255)
                                                                      acaTTGGTTTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCtt                                                                 <  2:786676/141‑1 (MQ=255)
                                                                         ttGGTGTAAGTTATGGCAAGTGCACAAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCt                                                              <  1:12433/141‑1 (MQ=255)
                                                                               tCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGt                                                        <  1:7370/141‑1 (MQ=255)
                                                                                cAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTGGtt                                                       >  1:273174/1‑141 (MQ=255)
                                                                                cAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTAGTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGtt                                                       <  2:87204/141‑1 (MQ=255)
                                                                                    tATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGCTTGAATGGCCCCTTCCTATTAGTTAAGt                                                   >  2:106060/1‑141 (MQ=255)
                                                                                     aTGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGAAATCATATATTATTCATCTTTAATTTAAAAAATGATTTTGCCATTGAAAAGTGAAATATGCCGGCATTATTTAGTGACCCCTTGCTATTAGtacatgt                                                  >  2:429268/1‑135 (MQ=255)
                                                                                          aaGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTtgtg                                             >  2:270567/1‑141 (MQ=255)
                                                                                          aaGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTCCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGGATGGCCCCGTTCTATTAGTTGAGTTTtgtg                                             <  1:693603/141‑1 (MQ=255)
                                                                                          aaGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTCCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTtgtg                                             <  1:263542/141‑1 (MQ=255)
                                                                                          aaGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTCCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTACTAGTTAAGTTTtgtg                                             <  1:600323/141‑1 (MQ=255)
                                                                                            gTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCg                                           >  2:123454/1‑141 (MQ=255)
                                                                                               cATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATGTTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTACGTTTTGTGCGTaa                                        >  2:288016/1‑141 (MQ=255)
                                                                                                  cAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCCCCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGc                                     <  1:735545/141‑1 (MQ=255)
                                                                                                  cAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGc                                     <  1:635507/141‑1 (MQ=255)
                                                                                                        tttCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGCTGTAGc                               >  1:117022/1‑141 (MQ=255)
                                                                                                         ttCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAGGCCTGAATATAACAAATCATAAGCCGTATGAATATTTAAATATTCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGCTGTAGCa                              >  2:619018/1‑141 (MQ=255)
                                                                                                         ttCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCCTGATTGAATGGCCCCTGTCTATTAGTTAAGTTTTGTGCGTAAAGCTGTAGCa                              >  1:480639/1‑141 (MQ=255)
                                                                                                         ttCAGACACAGCAAGTCTTACGGCATCTTCTATTTGTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCGACCATGATTGGAGGGCGCCTTTCTGTTAGTTAAGTTGTGTGCGTAAAGCTGTAGCa                              >  1:290373/1‑141 (MQ=255)
                                                                                                               cacaGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGCTTTGTGCGTAAAGCTGTAGCAAGATGc                        >  1:886815/1‑141 (MQ=255)
                                                                                                                  aggaaGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGCTGTAGCAAGTTGCTCa                     <  2:27803/138‑1 (MQ=255)
                                                                                                                           tACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAGGCTGTAAAAAGTTGCTCAAATTCAttc            >  2:234654/1‑139 (MQ=255)
                                                                                                                             cGGCATCTTCTATTTTTAAGCTTTAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGCTGTAGCAAGTTGCTCAAATTCATCCCa          <  1:191050/141‑1 (MQ=255)
                                                                                                                                gaTCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGCTGTAGCAAGTTGCTCAAATTCATCCCAGCt       <  1:120141/140‑1 (MQ=255)
                                                                                                                                   cttctATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGCTGTAGCAAGTTGCTCAAATTCATCCCAGCTGTc    >  2:872497/1‑141 (MQ=255)
                                                                                                                                   cttctATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTATGCGTAAAGCTGTAGCAGGTTGCGCAAATTCATCCCATCTGTc    >  1:122798/1‑141 (MQ=255)
                                                                                                                                     tctATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGCTGTAGCAAGTTGCTCAAATTCATCCCAGCTGTCCa  >  2:374010/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                             
TGTCATACCAAAAATCGGTTTAATGTTATTTGCTATACATGCATCATAAAATTTAGGAAAACCATACAATACATTGGTGTCAGTTATGGCAAGTGCATCAACATTTTCAGACACAGCAAGTCTTACGGCATCTTCTATTTTTAAGCTTGAATTTAACAAATCATAAGCCGTATGAATATTTAAATATGCCACCATGATTGAATGGCCCCTTTCTATTAGTTAAGTTTTGTGCGTAAAGCTGTAGCAAGTTGCTCAAATTCATCCCAGCTGTCCA  >  CP000730/1815047‑1815320

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: