Predicted mutation
evidence seq id position mutation annotation gene description
RA CP028101 875,381 A→G K85E (AAA→GAA)  C7Y58_04335 → autotransporter domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP028101875,3810AG98.0% 155.4 / ‑3.6 51K85E (AAA→GAA) C7Y58_04335autotransporter domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (19/31);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01

GCAAATGGAGTACCAACAAGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAAAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGCCTACAAATAGTGAAAATTATTCATTATTGCCTAAAT  >  CP028101/875273‑875627
                                                                                                             |                                                                                                                                                                                                                                                      
gCAAATGGAGTACCAACAAGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGAAAAATCACCTTGGATGTCATGGCAATTTGGAATATACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGgaga                                                                                                          >  1:576470/1‑251 (MQ=14)
 cAAATGGAGGACCAACAAGAGGAGGAAAT‑GCGACAACAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTGAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCCCCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATaaaggaaaagg                                                                                                             <  2:504248/247‑1 (MQ=11)
    aTGGAGTACCAACAAGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAg                                                                                                      >  2:237283/1‑251 (MQ=14)
       gAGTACCAACAAGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAGGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCa                                                                                                   >  2:27188/1‑251 (MQ=14)
       gAGTACCAACAAGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCa                                                                                                   >  1:177107/1‑251 (MQ=14)
       gAGTACCAACAAGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATAGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGc                                                                                                    <  2:190706/250‑1 (MQ=14)
           aCCAACAAGAG‑AAGAAATCGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCATAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACTTAGGATGTCATGGCAATTTGGAATAAACTATATGTAGAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGaaa                                                                                               >  1:569004/1‑251 (MQ=255)
           aCCAACAAGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGaaa                                                                                               >  1:459791/1‑251 (MQ=14)
           aCCAACAAGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGaaa                                                                                               <  2:33953/251‑1 (MQ=14)
            ccAACAAGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGaaaa                                                                                              <  1:15251/251‑1 (MQ=14)
                 aaGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATc                                                                                         <  1:117951/251‑1 (MQ=14)
                  agag‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATc                                                                                         >  1:147733/1‑250 (MQ=14)
                                   cATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATAttt                                                                        <  1:439333/251‑1 (MQ=18)
                                          agaAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGTCGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAag                                                                   <  1:616211/249‑1 (MQ=18)
                                          agaAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAGGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAag                                                                   <  1:27188/249‑1 (MQ=18)
                                             aaaTTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCAtt                                                                                     <  2:301204/228‑1 (MQ=14)
                                             aaaTTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCAtt                                                                                     >  1:301204/1‑228 (MQ=14)
                                             aaaTTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATa                                                              >  1:405730/1‑251 (MQ=21)
                                                 ttAAGAATCTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACaaa                                                          <  1:664291/251‑1 (MQ=37)
                                                  tAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAAt                                                         <  2:625340/251‑1 (MQ=37)
                                                   aaGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGTAGCAAGAGCAGAAAACTCAGCAAGTTTAGCAGGCTTCAGATTCGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAAGTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGTGATACAAAtt                                                        <  1:517956/251‑1 (MQ=11)
                                                        aCTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTAtta                                                   >  2:695711/1‑250 (MQ=38)
                                                        aCTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTGTAGCAGGCTTAAGATTAGAATTAATACAATTAATTGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGTAAAAGGGGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTc                                                   <  2:25505/251‑1 (MQ=255)
                                                             gTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAg                                              >  2:2024/1‑251 (MQ=255)
                                                             gTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAg                                              >  2:159024/1‑251 (MQ=255)
                                                              tAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAg                                              <  2:156196/250‑1 (MQ=255)
                                                              tAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGACCTCGCGGCCTTGTTTTCTTCTCTCTATTTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAg                                              <  1:578150/250‑1 (MQ=14)
                                                                 gAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGTAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAACTAGTAAAATCACGTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGata                                           <  1:620280/250‑1 (MQ=255)
                                                                 gAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGata                                           <  2:576470/250‑1 (MQ=255)
                                                                  aaGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGatata                                         <  2:624492/251‑1 (MQ=255)
                                                                  aaGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGatata                                         <  2:128490/251‑1 (MQ=255)
                                                                  aaGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGatat                                          <  2:106315/250‑1 (MQ=255)
                                                                     tttGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATTCAAATTTATTCAATAGATATATGc                                      <  2:577919/251‑1 (MQ=255)
                                                                     tttGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATg                                       <  1:347765/250‑1 (MQ=255)
                                                                     tttGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGc                                      <  1:241786/251‑1 (MQ=255)
                                                                     tttGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGc                                      <  2:147733/251‑1 (MQ=255)
                                                                     tttGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGc                                      <  2:44323/251‑1 (MQ=255)
                                                                      ttGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGcc                                     <  2:523440/251‑1 (MQ=255)
                                                                      ttGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGcc                                     <  1:337326/251‑1 (MQ=255)
                                                                           aTCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAg                                              <  2:83794/237‑1 (MQ=255)
                                                                           aTCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAg                                              >  1:83794/1‑237 (MQ=255)
                                                                           aTCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGata                                           <  1:153657/240‑1 (MQ=255)
                                                                           aTCTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGata                                           >  2:153657/1‑240 (MQ=255)
                                                                             cTAAGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGCCTACAAAt                              >  1:432845/1‑251 (MQ=255)
                                                                               aaGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGCCTACAAATAg                            >  1:322946/1‑251 (MQ=255)
                                                                               aaGATAGACGAAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGCCTACAAATAg                            <  1:184522/251‑1 (MQ=255)
                                                                                 gATAGACGCAGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCGGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGCCTACAAATAGta                          >  2:219714/1‑250 (MQ=255)
                                                                                         aaGCAAGAGCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGCCTACAAATAGTGCAAattat                  >  1:653826/1‑251 (MQ=255)
                                                                                                agCAGAAAACTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGCCTACAAATAGTGAAAATTATTCatt             <  2:322946/249‑1 (MQ=255)
                                                                                                         cTCAGAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGCCTACAAATAGTGAAAATTATTCATTATTGCCTAAAt  >  2:519279/1‑251 (MQ=255)
                                                                                                         cTCAAAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAGAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGCCTACAAATAGTGAAAATTATTCATTATTGCCTAAAt  >  2:216897/1‑251 (MQ=255)
                                                                                                             |                                                                                                                                                                                                                                                      
GCAAATGGAGTACCAACAAGAG‑AAGAAATAGCGACATCAAGAGAAAATTTAAGAAACTCAGTAGGAAGTTTGCAATCTAAGATAGACGAAGCAAGAGCAGAAAACTCAAAAAGTTTAGCAGGCTTAAGATTAGAATTAATACAATTAATGGAACAAGGGGATCAAGTAGTAAAATCACCTTGGATGTCATGGCAATTTGGAATAAACTATATGTACAGTAAATGGAATGGTACATATAAAGGAAAAGGAGATAAGGCAGAAAAATATCCATTTGAAGGAATATTTACAAGAGATACAAATTTATTCAATAGATATATGCCTACAAATAGTGAAAATTATTCATTATTGCCTAAAT  >  CP028101/875273‑875627

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑