Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,429,155 A→G N36S (AAC→AGC)  stfR → putative prophage side tail fiber protein StfR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,1550AG98.0% 147.4 / ‑4.5 50N36S (AAC→AGC) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (32/16);  total (32/18)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

CCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCG  >  NC_000913/1428924‑1429390
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
ccccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAAt                                                                                                                                                                                                                            <  1:8274/251‑1 (MQ=255)
            cGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGcc                                                                                                                                                                                                               >  2:103470/1‑251 (MQ=255)
             gACTGGGGCACTGACGGTTTAATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCg                                                                                                                                                                                                              <  2:202779/251‑1 (MQ=255)
              aCTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCgg                                                                                                                                                                                                             >  1:304401/1‑251 (MQ=255)
                tGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGGGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCgggg                                                                                                                                                                                                           >  2:959/1‑250 (MQ=255)
                         gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCCCGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGcc                                                                                                                                                                                                  >  1:236747/1‑250 (MQ=255)
                           cGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATg                                                                                                                                                                                                <  1:350915/251‑1 (MQ=255)
                                                          tGATTAACCAGGTAACACAGTGTTATGAAAGCCCGCCGGTTCAGCCGGGCATTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGAGTTATTc                                                                                                                                                                 >  1:144584/1‑251 (MQ=255)
                                                                 ccAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGGGGGGTGAATATGGCAGTAAAGATTTCAGGTGTAATGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACCGCGTTAATCTGTtggt                                                                                                                                                          >  2:283406/1‑251 (MQ=255)
                                                                  cAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTtggtg                                                                                                                                                         >  2:133380/1‑251 (MQ=255)
                                                                          acaGTGATATGACAGCCCGTCGGTTCCGGCGGGCTTTTTTGTGGGGGGAATATGGCAGTAAAGATTTCAGGTGTACTGCAAGACGGCACAGGAAAAACGGTACAGCACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTTGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGGTGAGTACGGTCCGTACCGCGGTATTCTGTtggtggctggact                                                                                                                                                 >  2:288918/1‑244 (MQ=255)
                                                                                          ccgccgGTTCAGGCGGGCTTTTTAGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACCCTGGCCTCAGAAAAACCGGATGAAGCCGGGCGTTAAAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGcc                                                                                                                                >  2:19387/1‑251 (MQ=255)
                                                                                          ccgccgGTTCAGGAGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGcc                                                                                                                                >  1:66164/1‑251 (MQ=255)
                                                                                            gccgGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGCAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGGGGAAGGATTCCCGCCGTCACATGCCgg                                                                                                                              >  2:230745/1‑251 (MQ=255)
                                                                                             ccgGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAACGATTTCAGGTGTACTGAAAGTCGGCACAGGAAAACCGGGAAAGAACTGCACAATACCGATGAAAGCCAAACGTAACAGCACCACGGCGGTGGGGAGCACGCTGGCCCCAGCAAAACCGGAAGAAGCCGGGCGTTACAGCATGGAAGTTGAGCACGGTCAGTACAGCGGTATTCTGGTGGTGGCAGGATTCACGCCGGCTCATTCCggg                                                                                                                             >  2:126933/1‑251 (MQ=255)
                                                                                                   cAGGCGGGTTTTTTTGTGGGGTGAATAGGGCAGTAAAGATTGCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCAGGGACGTGGAGTACGTTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCAtt                                                                                                                       <  2:343880/251‑1 (MQ=255)
                                                                                                      gCGGGTTTTTTTGGGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAATTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAcc                                                                                                                    <  1:57813/251‑1 (MQ=255)
                                                                                                        gggCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCgt                                                                                                                  >  2:355160/1‑251 (MQ=255)
                                                                                                               tttGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGa                                                                                                            >  2:63570/1‑250 (MQ=255)
                                                                                                                       gTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATtctc                                                                                                    >  2:220769/1‑250 (MQ=255)
                                                                                                                               ggCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGc                                                                                                                                                                                                                <  2:216344/135‑1 (MQ=255)
                                                                                                                               ggCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGc                                                                                                                                                                                                                >  1:216344/1‑135 (MQ=255)
                                                                                                                               ggCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTa                                                                                           <  1:100681/251‑1 (MQ=255)
                                                                                                                                        gATTTCAGGTGTACTGAAAGACGGCACAGGAAACCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATg                                                                                  >  2:156994/1‑251 (MQ=255)
                                                                                                                                             cAGGTGTACTGAAAGACGGCACAGGCAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTAt                                                                                                                                                                   >  2:48670/1‑166 (MQ=255)
                                                                                                                                             cAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGGACAGCGTTAt                                                                                                                                                                   <  1:48670/166‑1 (MQ=255)
                                                                                                                                                 ggtACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCAAAACCCAGCAGAAAGCAAAACGAATCAGCTCCAGGGGGGGGGTGAACCCGCTGGCCTGAGAAAATCCGGATGAGGCCGGGCGTTCCAGCATGGAGGTTGAGTACGGTCAGTACAGTGTAGTTCTGTTGTGGTAAGGATCACCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAAGGATTGTCTCg                                                                         <  2:342319/250‑1 (MQ=255)
                                                                                                                                                  gtACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAAAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCg                                                                         <  1:294110/250‑1 (MQ=255)
                                                                                                                                                   tACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGt                                                                       >  2:343327/1‑251 (MQ=255)
                                                                                                                                                       gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCGACCCGGTACGCTGAATGATTTTCTCGGTGc                                                                     <  1:959/249‑1 (MQ=255)
                                                                                                                                                       gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGAGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGcc                                                                    >  2:292924/1‑250 (MQ=255)
                                                                                                                                                         aaGACGGCACAGGAAAACCGGTACAGAAATGCACAATCCAGCTGAAAGCACAACGTAACAGCCCCACGGAGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTAACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGCTTTTCTCGGTGCCATg                                                                 >  2:137765/1‑251 (MQ=255)
                                                                                                                                                                cacaGGAAACCCGGTACAGAACTGCACAATCCAGCTGAAAGAAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCg                                                                                              <  2:342395/215‑1 (MQ=255)
                                                                                                                                                                cacaGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCg                                                                                              >  1:342395/1‑215 (MQ=255)
                                                                                                                                                                   agaaaacccgggcaaGAACTGCCCAATCCAGCTGAAAGAAAATCGTAACAGCACCACGGTGGGGGTGAGCACTCTGGCCTCAAAAACCACGGATGAGGCGGGGCGTTACATCATGGACGTTGAGTCCGGTCAGTACAGCGTTATTCTGTGGGTGGAAGTATCCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGgatg                                                       <  2:284529/237‑1 (MQ=255)
                                                                                                                                                                      aaaaCCGGTACAGCACTGCATAATCCAGCTGACTGCAAAACGTAACAGCGCCCCGGAGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGCTGAAGCCGGGCGTTACAGCATGGCCGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCAGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGCTTTTCTCGGTGCCATGCCGGAGGctgatg                                                    >  2:355757/1‑251 (MQ=255)
                                                                                                                                                                         aCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGccc                                                 >  1:329341/1‑251 (MQ=255)
                                                                                                                                                                           cGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGt                                               >  1:360220/1‑251 (MQ=255)
                                                                                                                                                                               aCAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCAGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAACATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATTACGGAGGATGATGCCCGTCCgg                                           >  1:237743/1‑251 (MQ=255)
                                                                                                                                                                                 ggAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCCGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCgg                                           >  2:126462/2‑249 (MQ=255)
                                                                                                                                                                                           caTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGatgatg                                                    >  2:319590/2‑230 (MQ=255)
                                                                                                                                                                                           aaTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGatgatg                                                    <  1:319590/230‑1 (MQ=255)
                                                                                                                                                                                              ccAGCTGAAAGCAAACAGTAACAGCACCACGGGGGTGGTGAGCCCGCTGGCCACAGAAAATCCGGAGGAAGCCGGGGGTTACAGCATGGACGTTGGGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATGCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGta                            <  2:145234/251‑2 (MQ=255)
                                                                                                                                                                                                 gCTGAAAGCAAAACGTAACAGCCCCACGGTGGTGGTGAGCACGCGGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATtct                                                                                                     <  2:159625/175‑1 (MQ=255)
                                                                                                                                                                                                 gCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATtct                                                                                                     >  1:159625/1‑175 (MQ=255)
                                                                                                                                                                                                   gaaaagaaaaaagaaccagcacaccGGTGGTGTTGAGCACGCTGGCCGCAGAAAAGCGGGAAGAAGCCGGGCGTTACAGAATGTACGTTGAGTACGGTCAGTACAGCGTTATTCTTTTGGTGGAAGGATCCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGaa                       <  2:326626/227‑1 (MQ=255)
                                                                                                                                                                                                     aaaGCAACACGTAACCGCACCACGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACt                     >  2:35476/1‑251 (MQ=255)
                                                                                                                                                                                                         cgaaaCGTAAGAGCACCACGGTGGTGGTGAGCACGCAGGACTCAGAAAATCCGGATGAAGCCGGGCGTTAAAGCATGGACGTTGCGTACGGTCAGTACAGCGCTATTATGTTGGTGGAAGGATTCCCGCCGTCCCATGCCGGGACAATTACCGCGTATGAAGATTCTCCACCCGGTACGATGAATGATTTTCTCGGAGCCATGACGGCGGATGATGCACGTCCGGAGGCACTGCGCCGTTTTGAACTGAtg                 >  2:145911/3‑251 (MQ=255)
                                                                                                                                                                                                                    aGCACCAGGGTGGTGGTGACAACGTGGGCCTCAGAAAACCCGGAGGAAGCCGGGCGTTACAACATGGAGTTTGAGTACGGTCAGTACAGGGTTATTCTTTTGGTGGAAGGATCCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGt      <  2:249707/251‑1 (MQ=255)
                                                                                                                                                                                                                        ccacGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGgcg  >  1:290904/1‑251 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
CCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCG  >  NC_000913/1428924‑1429390

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