Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,835,011 T→A 7.3% I90F (ATC→TTC)  U712_19780 ← ATP‑dependent RNA helicase dbpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,835,0110TA 7.3% 195.8 / 6.8 69I90F (ATC→TTC) U712_19780ATP‑dependent RNA helicase dbpA
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base A (2/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

GATATGATCCAGAACACGTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGCTTCCTGTCTGTGATTT  >  CP006881/3834874‑3835146
                                                                                                                                         |                                                                                                                                       
gATATGATCCAGAACACGTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGat                                                                                                                                         <  5:65748/138‑1 (MQ=255)
 atatGATCCCGAACACGGCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGatat                                                                                                                                       <  8:258704/139‑1 (MQ=255)
 atatGATCCAGAACACGTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGatat                                                                                                                                       <  7:180301/139‑1 (MQ=255)
 atatGATCCAGAACACGTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGatat                                                                                                                                       <  5:87168/139‑1 (MQ=255)
  tatGATCCAGAACACGTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATc                                                                                                                                      <  8:78176/139‑1 (MQ=255)
   atGATCCAGATCAGGTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATc                                                                                                                                      <  6:259407/138‑1 (MQ=255)
      aTCCAGAACACGTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTcc                                                                                                                                  <  5:240222/139‑1 (MQ=255)
                 gTCCCGGCGTGCCGGCAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCAc                                                                                                                       >  2:286254/1‑139 (MQ=255)
                 gTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTTATCCGTTTAAATCGTCCGATGTTCGTGATATGTTACTTTACTTGCAc                                                                                                                       >  6:376504/1‑139 (MQ=255)
                 gTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCAc                                                                                                                       >  8:96483/1‑139 (MQ=255)
                 gTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCAc                                                                                                                       >  8:79108/1‑139 (MQ=255)
                 gTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCAc                                                                                                                       >  1:33375/1‑139 (MQ=255)
                 gTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCATTACATGCAc                                                                                                                       >  1:100624/1‑139 (MQ=255)
                 gTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATAAAATGAGGATTTTCCAAATACGGCTGTCGCTTTTTTACGTTTAAATCGTCCGATGTTAGTGTTATGTTACTTTAAtttgac                                                                                                                       >  5:357655/1‑135 (MQ=255)
                    ccGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTTATCCCTTTAAATCGTCCGATTTTAGTGATATATTCCTTTACTTTCAc                                                                                                                       >  4:140889/1‑136 (MQ=255)
                    ccGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTGGTGTTATGTTCCTTTACTTGCAc                                                                                                                       >  3:34878/1‑136 (MQ=255)
                    ccGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCAc                                                                                                                       >  4:303107/1‑136 (MQ=255)
                    ccGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCAc                                                                                                                       <  3:23485/136‑1 (MQ=255)
                    ccGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATACGTTTAAATAGTCCGATGTTAGTGATATCTTACTTTACTTTCAc                                                                                                                       >  4:23485/1‑136 (MQ=255)
                      ggCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGTATGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCTTCCGATGTTAGTGATATATTCCTTTAATTTACCCGCaa                                                                                                                  >  8:184612/1‑139 (MQ=255)
                      ggCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGc                                                                                                                         >  7:19250/1‑132 (MQ=255)
                      ggCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGc                                                                                                                         >  7:12953/1‑132 (MQ=255)
                      ggCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGc                                                                                                                         <  8:12953/132‑1 (MQ=255)
                      ggCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGc                                                                                                                         <  8:19250/132‑1 (MQ=255)
                         gTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCaaccg                                                                                                               >  7:20114/1‑136 (MQ=255)
                                caacacTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        <  2:162551/139‑1 (MQ=255)
                                caacaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTATATAGTACGATGTTAGTGATATCTTCCTTTACTTGCACCGCAAGCTCACgtgt                                                                                                        <  8:27385/139‑1 (MQ=255)
                                caacaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACg                                                                                                           <  3:303107/136‑1 (MQ=255)
                                caacaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACg                                                                                                           <  4:258068/136‑1 (MQ=255)
                                caacaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        <  5:376504/139‑1 (MQ=255)
                                caacaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        <  1:286254/139‑1 (MQ=255)
                                caacaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        <  6:167892/139‑1 (MQ=255)
                                caacaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        <  8:343457/139‑1 (MQ=255)
                                caacaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        <  7:373922/139‑1 (MQ=255)
                                caacaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        <  1:129285/139‑1 (MQ=255)
                                caacaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        <  1:70461/139‑1 (MQ=255)
                                  acaaTGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCg                                                                                                      <  6:129603/139‑1 (MQ=255)
                                      tgtgGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCTAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGgtgt                                                                                                  <  1:75130/139‑1 (MQ=255)
                                      tgtgGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCCCCGCAAGCTCACGGGTCGgtgt                                                                                                  >  5:176574/1‑139 (MQ=255)
                                      tgtgGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        <  7:29130/133‑1 (MQ=255)
                                      tgtgGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        >  8:29130/1‑133 (MQ=255)
                                              ttttGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTACTTTACTTGCACCGCAAGCTCACGGGGGGGTGTTAAAATa                                                                                           >  1:100606/1‑138 (MQ=255)
                                                ttGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAgc                                                                                        >  7:43687/1‑139 (MQ=255)
                                                            aGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCg                                                                                 <  3:140889/134‑1 (MQ=255)
                                                             gCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        <  7:23378/110‑1 (MQ=255)
                                                             gCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGt                                                                                                        >  8:23378/1‑110 (MQ=255)
                                                                 tttGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGt                                                                           <  3:212359/135‑1 (MQ=255)
                                                                  ttGTTTATCAAAGGAGGATTTTCCAAATACGTCTGTCGCTTTGATCCGTTTAAATAGTACGCTGTTCGTTATATAGTACTTTACATTCACCCCAAGCTAACGGGTCGGTTTTAAAATAATCGCCTGAGGCTAGttttgc                                                                      >  6:41489/1‑137 (MQ=255)
                                                                   tGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCAt                                                                     >  1:269511/1‑139 (MQ=255)
                                                                        aTCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCaa                                                                >  8:88526/1‑139 (MQ=255)
                                                                            aaGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCg                                                            <  7:96483/139‑1 (MQ=255)
                                                                                cggaTTTTCCAAATACGGAGGGCGCTTTGATCCGTTTAAAGCGGCCGATGTTCGTGATCACTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGGGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAg                                                        <  6:225661/138‑1 (MQ=255)
                                                                                  gaTTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCAGCCCAATTCGCCAGcg                                                      >  5:190419/1‑138 (MQ=255)
                                                                                    ttttCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTca                                                    <  7:79108/139‑1 (MQ=255)
                                                                                       tCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTca                                                    >  1:34598/1‑136 (MQ=255)
                                                                                       tCCAAATACGGCTGGCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCATTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTcaca                                                  <  7:184612/138‑1 (MQ=255)
                                                                                       cccAAAAAGGGCGGCCGCTTTTACCCGTTAAAATCGTCCGAGTTGCGTTAACTCTACCTTAACTTCCCCCGCATGCGCCGGGGGCGGTGTAAAAATACGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTcaca                                                  <  2:216332/137‑1 (MQ=255)
                                                                                           aaTACGGCTGTCGCTTTGAACCGTTTAAAACGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTAAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACagagag                                             <  6:329888/139‑1 (MQ=255)
                                                                                            aTACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACagagag                                             >  8:170865/1‑138 (MQ=255)
                                                                                            aTACCGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCATACATTACTTTCACCGCAAGATCAAGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTAATACCAAATCGCCAGCTCACagagag                                             >  2:81583/1‑138 (MQ=255)
                                                                                                                  tttAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTc                      >  7:8238/1‑139 (MQ=255)
                                                                                                                  tttAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTc                      >  1:345915/1‑139 (MQ=255)
                                                                                                                     aaaTCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTc                      >  3:428697/1‑136 (MQ=255)
                                                                                                                       aTCGTCCGAGGTTGGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGCt                 <  8:414163/139‑1 (MQ=255)
                                                                                                                         cGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGc                  <  4:428697/136‑1 (MQ=255)
                                                                                                                           tCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGTATTGTTTTCATCCCAATTCGCCATCTCACAGAGAGGAATACCGAACGAAGCCGGTTTCCCGCTTc               >  7:297317/1‑137 (MQ=255)
                                                                                                                           tCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGCTTCct             >  6:59151/1‑139 (MQ=255)
                                                                                                                              gATGTTCGTGATATCTTCATTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGACCGAAGCCGTGTTCCCGCTTCct             >  4:66774/1‑136 (MQ=255)
                                                                                                                               aTGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGCTTCctg            >  3:345176/1‑136 (MQ=255)
                                                                                                                                    cGTGAAATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGCTTCctgtctgt       <  1:3519/136‑1 (MQ=255)
                                                                                                                                    cGTGAAATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGCTTCctgtctgt       >  2:3519/1‑136 (MQ=255)
                                                                                                                                    cGTGAAATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGCTTCctgtctgt       <  1:10268/136‑1 (MQ=255)
                                                                                                                                    cGTGAAATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGCTTCctgtctgt       >  2:10268/1‑136 (MQ=255)
                                                                                                                                         tatCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGCTTCCTGTCTGTGAttt  <  3:98972/136‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
GATATGATCCAGAACACGTCCCGGCGTGCCGACAACAATGTGGCTTTTTTGCTTAAGCTCAGCTTTTTGTTTATCAAAGGAGGATTTTCCAAATACGGCTGTCGCTTTGATCCGTTTAAATCGTCCGATGTTCGTGATATCTTCCTTTACTTGCACCGCAAGCTCACGGGTCGGTGTTAAAATAAGCGCCTGCGGCTTGTTTTCATCCCAATTCGCCAGCTCACAGAGAGGAATCCCGAACGAAGCCGTTTTCCCGCTTCCTGTCTGTGATTT  >  CP006881/3834874‑3835146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑