Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 733,589 T→G L1336* (TTA→TGA)  rhsC → Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913733,5890TG98.7% 281.0 / ‑4.0 77L1336* (TTA→TGA) rhsCRhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (35/41);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

GCGGAACTGGTGAAAGAAGCGGTCCAGATATATGGGGGAATCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTAT  >  NC_000913/733361‑733807
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                          
gCGGAACTGGTGAAAGAAGCGGTCCAGATATAGGGGGGAATCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAggg                                                                                                                                                                                                                     <  1:389291/236‑1 (MQ=255)
     aCTGGTGAAAGAAGCGGTCCAGATATATGGGGGAATCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGtt                                                                                                                                                                                                                >  1:493676/1‑236 (MQ=255)
                     gTCCAGATATATGGGGGAATCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCcg                                                                                                                                                                                                >  2:633066/1‑236 (MQ=255)
                       ccAGATATATGGGTGACTCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCcgcg                                                                                                                                                                                              <  1:119571/236‑1 (MQ=255)
                       ccAGATATATGGGGGAATCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCcgcg                                                                                                                                                                                              <  1:633066/236‑1 (MQ=255)
                           atatatGGGGGAATCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTACAGAGGTAAGGCCGCGt                                                                                                                                                                                             >  2:425079/1‑233 (MQ=255)
                                tGGGGGAATCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGt                                                                                                                                                                                     >  1:636050/1‑236 (MQ=255)
                                    ggAATCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTaa                                                                                                                                                                                                             <  2:416261/208‑1 (MQ=255)
                                    ggAATCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTaa                                                                                                                                                                                                             >  1:416261/1‑208 (MQ=255)
                                      aaTCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGc                                                                                                                                                                                <  1:560544/235‑1 (MQ=255)
                                        tCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTg                                                                                                                                                                             >  2:442228/1‑236 (MQ=255)
                                        tCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTg                                                                                                                                                                             >  1:96358/1‑236 (MQ=255)
                                          cGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGat                                                                                                                                                                           <  2:22166/236‑1 (MQ=255)
                                          cGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGat                                                                                                                                                                           <  1:530292/236‑1 (MQ=255)
                                               tATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATAtt                                                                                                                                                                      <  2:96358/236‑1 (MQ=255)
                                                       atatCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCta                                                                                                                                                              <  1:116171/236‑1 (MQ=255)
                                                             ttGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAAtt                                                                                                                                                        >  1:660988/1‑236 (MQ=255)
                                                             ttGTGTCCCAGACGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAAtt                                                                                                                                                        >  1:125338/1‑236 (MQ=255)
                                                               gtgtCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTa                                                                                                                                                      >  1:506686/1‑236 (MQ=255)
                                                                    ccAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGat                                                                                                                                                                           <  1:102512/210‑1 (MQ=255)
                                                                    ccAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGat                                                                                                                                                                           >  2:102512/1‑210 (MQ=255)
                                                                            tAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGttt                                                                                                                                                                                                               >  2:687311/1‑166 (MQ=255)
                                                                            tAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGttt                                                                                                                                                                                                               <  1:687311/166‑1 (MQ=255)
                                                                                      tataCTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGc                                                                                                                                                                                <  2:586460/187‑1 (MQ=255)
                                                                                      tataCTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGc                                                                                                                                                                                >  1:586460/1‑187 (MQ=255)
                                                                                      tataCTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGc                                                                                                                                                  <  2:766334/217‑1 (MQ=255)
                                                                                      tataCTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGc                                                                                                                                                  >  1:766334/1‑217 (MQ=255)
                                                                                       ataCTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGa                                                                                                                               <  1:157627/235‑1 (MQ=255)
                                                                                        taCTTGTGGGGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGAc                                                                                                                              <  2:527462/235‑1 (MQ=255)
                                                                                          cTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTgg                                                                                                                           <  1:442228/236‑1 (MQ=255)
                                                                                             gtggtgGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCcg                                                                                                                                                                                                >  2:338823/1‑164 (MQ=255)
                                                                                             gtggtgGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCcg                                                                                                                                                                                                <  1:338823/164‑1 (MQ=255)
                                                                                             gtggtgGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCt                                                                                                                        >  1:603385/1‑236 (MQ=255)
                                                                                              tggtggCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTa                                                                                                                       >  1:457338/1‑236 (MQ=255)
                                                                                                gtggCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTaaaa                                                                                                                                         <  2:391161/216‑1 (MQ=255)
                                                                                                gtggCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTaaaa                                                                                                                                         >  1:391161/1‑216 (MQ=255)
                                                                                                  ggCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGc                                                                                                                                                  >  2:461479/1‑205 (MQ=255)
                                                                                                  ggCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGc                                                                                                                                                  <  1:461479/205‑1 (MQ=255)
                                                                                                  ggCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGa                                                                                                                               <  1:86999/224‑1 (MQ=255)
                                                                                                  ggCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGa                                                                                                                               >  2:86999/1‑224 (MQ=255)
                                                                                                  ggCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGa                                                                                                                   <  1:328624/236‑1 (MQ=255)
                                                                                                  ggCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGa                                                                                                                   <  1:510757/236‑1 (MQ=255)
                                                                                                     cAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGc                                                                                                                                                                                <  1:425651/172‑1 (MQ=255)
                                                                                                     cAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGc                                                                                                                                                                                >  2:425651/1‑172 (MQ=255)
                                                                                                            aaCGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGaa                                                                                                         <  2:660988/236‑1 (MQ=255)
                                                                                                             aCGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGaa                                                                                                         <  2:125338/235‑1 (MQ=255)
                                                                                                                gggAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGGTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACAt                                                                                                     >  1:353733/1‑236 (MQ=255)
                                                                                                                     aaTCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTAtt                                                                                                                                                                                         <  1:125667/147‑1 (MQ=255)
                                                                                                                     aaTCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTAtt                                                                                                                                                                                         >  2:125667/1‑147 (MQ=255)
                                                                                                                          aaaGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGaa                                                                                           <  1:425079/236‑1 (MQ=255)
                                                                                                                                 aTGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTcc                                                                                    >  1:645740/1‑236 (MQ=255)
                                                                                                                                     tCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATCTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGaa                                                                                >  2:658515/1‑236 (MQ=255)
                                                                                                                                      cTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGaaa                                                                               <  2:550839/236‑1 (MQ=255)
                                                                                                                                               ccAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGa                                                                                                                               <  1:268869/179‑1 (MQ=255)
                                                                                                                                               ccAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGa                                                                                                                               >  2:268869/1‑179 (MQ=255)
                                                                                                                                               ccAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAAc                                                                                                       <  2:333216/203‑1 (MQ=255)
                                                                                                                                               ccAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAAc                                                                                                       >  1:333216/1‑203 (MQ=255)
                                                                                                                                               ccAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTgg                                                                      >  1:362870/1‑236 (MQ=255)
                                                                                                                                               ccAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTgg                                                                      >  2:367311/1‑236 (MQ=255)
                                                                                                                                                        gCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTttat                                                             <  1:658515/236‑1 (MQ=255)
                                                                                                                                                             tAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTttatt                                                            >  2:328615/1‑232 (MQ=255)
                                                                                                                                                                 gATACAAAATAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGAttt                                                      <  2:570091/234‑1 (MQ=255)
                                                                                                                                                                  aTACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGt                                                    <  2:362870/235‑1 (MQ=255)
                                                                                                                                                                   tACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTcc                                                                                    >  1:371995/1‑202 (MQ=255)
                                                                                                                                                                   tACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTcc                                                                                    <  2:371995/202‑1 (MQ=255)
                                                                                                                                                                       aaaGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCAc                                              >  1:766936/1‑236 (MQ=255)
                                                                                                                                                                               tgctgGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGa                                          <  2:353733/232‑1 (MQ=255)
                                                                                                                                                                                   ggCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTa                                  <  2:506686/236‑1 (MQ=255)
                                                                                                                                                                                    gCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTaa                                 >  1:280109/1‑236 (MQ=255)
                                                                                                                                                                                     ccGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATg                                    <  1:328615/232‑1 (MQ=255)
                                                                                                                                                                                          atcagcTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAat                           <  1:408687/231‑1 (MQ=255)
                                                                                                                                                                                             aCCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTaa                                                                                                                      <  1:354869/142‑1 (MQ=255)
                                                                                                                                                                                             aCCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTaa                                                                                                                      >  2:354869/1‑142 (MQ=255)
                                                                                                                                                                                                            ataataGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCa         >  2:237282/1‑236 (MQ=255)
                                                                                                                                                                                                             taataGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCag        <  2:766936/236‑1 (MQ=255)
                                                                                                                                                                                                                taGTTGTGTGGAGAACTGTTGAAAAGGGAAGTTTAAAGAGGTAAGGCCGAGTTATTCTGTATTGCCTGATATTTTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCa         <  2:636050/232‑1 (MQ=255)
                                                                                                                                                                                                                           agaACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCtat  >  1:410575/1‑228 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                          
GCGGAACTGGTGAAAGAAGCGGTCCAGATATATGGGGGAATCCGTTCTATCATCAATATCTTTGTGTCCCAGATGGTAAAGGGGACTATACTTGTGGTGGCCAAGACCAACGGGGAGAATCAAAAGGAGATGGTCTATGGGGGCCAGGTAAAGCAAGTAATGATACAAAAGAAGCTGCTGGCCGTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTAT  >  NC_000913/733361‑733807

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