Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 525,588 A→G 5.1% K776K (AAA→AAG rhsD → Rhs protein with DUF4329 family 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_000913525,5880AG 5.1% 189.4 / 2.7 78K776K (AAA→AAGrhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base G (2/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

ACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGC  >  NC_000913/525440‑525733
                                                                                                                                                    |                                                                                                                                                 
aCGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTCTGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGTAACTTCTGTGGCAGCATGAGACGaaa                                                                                                                                                   <  1:212032/149‑1 (MQ=255)
aCGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGaaa                                                                                                                                                   >  1:404657/1‑149 (MQ=255)
aCGGCGGGCTGACAGACAGCCGCCACCTGTGCGAAGGCCACCGTGGTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGGGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGaaa                                                                                                                                                   <  2:149615/149‑1 (MQ=18)
   gctggctgACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAacac                                                                                                                                                <  1:837331/149‑1 (MQ=255)
     tggctgACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGc                                                                                                                                              <  1:304648/149‑1 (MQ=255)
        ctgACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATa                                                                                                                                           >  1:189198/1‑149 (MQ=255)
         tgACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGCCGAAACACGCATAc                                                                                                                                          >  2:45586/1‑149 (MQ=255)
            cagacaTCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAAc                                                                                                                                       <  2:148578/149‑1 (MQ=255)
             agacaTCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGCACCCGGAGCCGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACg                                                                                                                                      >  1:1022918/1‑149 (MQ=255)
             agacaTCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACg                                                                                                                                      >  2:144575/1‑149 (MQ=255)
               acaTCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAg                                                                                                                                    >  1:727281/1‑149 (MQ=255)
                 aTCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAGCACGCATACAACGAGCa                                                                                                                                  >  1:92498/1‑149 (MQ=255)
                 aTCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAGCACGCATACAACGAGCa                                                                                                                                  >  1:96095/1‑149 (MQ=255)
                 aTCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCa                                                                                                                                  >  2:45020/1‑149 (MQ=255)
                     gCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGCCGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGCAACACGCATCCAACGAGCAgggg                                                                                                                              >  1:76212/1‑149 (MQ=255)
                     gCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAgggg                                                                                                                              >  1:505874/1‑149 (MQ=255)
                        aCCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTg                                                                                                                           >  1:481463/1‑149 (MQ=255)
                         ccTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTgg                                                                                                                          <  2:462391/149‑1 (MQ=255)
                         acTGAGCGAAGGCCACCGTGTTGCCTTCCTCTATGGCTATGACTATAAATGCCGCCTGACCGGCGAATGCCATACTGTGGAGAACCCGGAGACGGGGGATCTGCTGTGGCAGCATGAGACTAAACACGCATACAACGAGCAGGGGCTgg                                                                                                                          <  2:486436/148‑1 (MQ=18)
                                  aGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcgcg                                                                                                                 <  1:144575/149‑1 (MQ=255)
                                   ggCCACCCTGTGCGCGTCCACTATGGCTATGACGATAACGGCCGCCTGACCGGCGATTGCCAGACGGTGGAGAACCCGGAGTCGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGt                                                                                                                <  1:751610/149‑1 (MQ=17)
                                      cACCGTGTTGTCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCAc                                                                                                             >  1:376105/1‑149 (MQ=255)
                                      cACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGACACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCAc                                                                                                             <  2:230381/149‑1 (MQ=255)
                                      cACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCAc                                                                                                             >  2:32872/1‑149 (MQ=255)
                                      cACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCAc                                                                                                             >  2:994835/1‑149 (MQ=255)
                                        ccGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTTGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGc                                                                                                           <  2:289922/149‑1 (MQ=255)
                                           tgtTGCCGTCCTCTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGTCGGTTGAGAACCCGGAGACGGGTGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCgg                                                                                                        <  2:944735/149‑1 (MQ=255)
                                           gtttGCCGTCCACTATGTGTATGACGATAAAGGGCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGATCTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCgg                                                                                                        <  1:1005091/147‑1 (MQ=255)
                                            gtTGCCGTCCACTATGGCTATGACGATAAAGGCCTCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACTCCTGa                                                                                                       >  2:743491/1‑149 (MQ=255)
                                               gCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGCGCAGGGGCTGGCAAACCGCGTCACGCCGGACAg                                                                                                    >  1:796509/1‑149 (MQ=255)
                                               gCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAg                                                                                                    >  2:172842/1‑149 (MQ=255)
                                               gCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAg                                                                                                    <  2:798061/149‑1 (MQ=255)
                                                ccGGCCACTATGGCTATGACGATAAAGGCCGCCGGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGc                                                                                                   <  1:517501/149‑1 (MQ=255)
                                                  gTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCt                                                                                                 <  1:1010100/149‑1 (MQ=255)
                                                       cTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTgccgc                                                                                            <  1:32872/149‑1 (MQ=255)
                                                         aTGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGccgccg                                                                                          >  2:330122/1‑149 (MQ=255)
                                                          tGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCCCGCCGGACCGCCTGCCGCCgg                                                                                         >  2:79122/1‑149 (MQ=255)
                                                                 gACGATAAAGGCCGCCTGACCGGCGAATGCCCGACGGGGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGGCACGCCGGACAGCCTGCCGCCGGTGGAGTg                                                                                  <  1:915476/149‑1 (MQ=255)
                                                                   cGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTTGAGAACCCGGAGACGGTGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGc                                                                                <  2:61669/149‑1 (MQ=255)
                                                                    gATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACCGCCTGCCGCCGGTGGAGTGGCt                                                                               >  2:896944/1‑149 (MQ=255)
                                                                      aaaaGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGATCTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGa                                                                             <  2:511424/148‑1 (MQ=255)
                                                                           gccgccTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTAt                                                                        >  2:77204/1‑149 (MQ=255)
                                                                              gccTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGc                                                                     >  1:644143/1‑149 (MQ=255)
                                                                                cTGTCCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAg                                                                   <  1:696481/149‑1 (MQ=255)
                                                                                cTGACCGGCTAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAGCACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAg                                                                   <  2:92498/149‑1 (MQ=255)
                                                                                cTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAGCACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAg                                                                   <  2:96095/149‑1 (MQ=255)
                                                                                  gACCGGCGAACGCCCGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGAGGCAGCAGGAGCCGAAACACGCAAACCACGCGCCGGGGCGGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCg                                                                 <  1:559862/149‑1 (MQ=17)
                                                                                     cGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGTACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGtt                                                              >  1:764123/1‑149 (MQ=255)
                                                                                       gCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGCCGTAGGGCAGCGGTTAc                                                            >  1:569202/1‑149 (MQ=255)
                                                                                       gCGAATGCCAGACGCTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACACCGAGCAGGGGCTGGCACACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTAc                                                            <  1:93390/149‑1 (MQ=255)
                                                                                             gCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCg                                                      <  2:727281/149‑1 (MQ=255)
                                                                                              ccAGACGGTTTAGAACCCTGAGACGGGGTAACTGCTGTGGCAGCATGAGACGAAACACTCATACAACGAGCAGCGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCgg                                                     <  1:979438/149‑1 (MQ=18)
                                                                                                 gACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAACCACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGaa                                                  <  1:295810/149‑1 (MQ=255)
                                                                                                 gACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGaa                                                  <  1:418869/149‑1 (MQ=255)
                                                                                                    ggTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCTGCATGAGACGAAAAACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCCCGCCGGCCTGCCTCCCGCAGGTGGCGTGGATTACCTAGGCCACCGCTAACCTGGCGGGAAGGa                                               >  2:669708/1‑149 (MQ=17)
                                                                                                      tGGAGACCCCGGAGCCGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGGGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAg                                             <  2:322973/149‑1 (MQ=255)
                                                                                                        gagaACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCt                                           >  1:199439/1‑149 (MQ=255)
                                                                                                          gaACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTgg                                         <  2:29553/149‑1 (MQ=255)
                                                                                                             cccGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCg                                      <  1:596669/149‑1 (MQ=255)
                                                                                                             cccGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCg                                      >  1:243730/1‑149 (MQ=255)
                                                                                                              ccGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCgg                                     >  1:315605/1‑149 (MQ=255)
                                                                                                              ccGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCgg                                     <  2:189198/149‑1 (MQ=255)
                                                                                                              ccGGAGACGGGGGAACAGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCACGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCgg                                     <  1:94194/149‑1 (MQ=255)
                                                                                                                ggAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGa                                   >  1:231050/1‑149 (MQ=255)
                                                                                                                 gagaCGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGAc                                  <  2:376105/149‑1 (MQ=255)
                                                                                                                    aCGGGGGAACTGCTGGGGCAGCATGAGACGAAACCCGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGAcgcc                               <  2:404657/149‑1 (MQ=255)
                                                                                                                        gggAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTg                           <  2:34937/149‑1 (MQ=255)
                                                                                                                        gggAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTg                           >  2:936573/1‑149 (MQ=255)
                                                                                                                         ggAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTgg                          <  1:442273/149‑1 (MQ=255)
                                                                                                                          gAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGt                         <  1:330122/149‑1 (MQ=255)
                                                                                                                           aaCTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGAATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTc                        >  1:193046/1‑149 (MQ=255)
                                                                                                                           aaCTGCTGGGGCAGCATGAGACGCAACACGCATACACCGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCGGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTc                        <  1:115138/149‑1 (MQ=255)
                                                                                                                             ctgctgTGGCATCAGGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGa                      <  2:88713/149‑1 (MQ=255)
                                                                                                                             ctgctgTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGa                      >  1:543030/1‑149 (MQ=255)
                                                                                                                             ctgccgTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCGGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGa                      <  1:108948/149‑1 (MQ=255)
                                                                                                                              tgctgTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAg                     <  2:528397/149‑1 (MQ=255)
                                                                                                                               tctgTGGCTGCATGAGACGAAACACGCTTAGAACGAGCAGTGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCACTGGTCGatt                    <  2:190184/148‑3 (MQ=18)
                                                                                                                                     gcagcaTGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGCCGCCGCTCGGCGAGTATActc              >  2:869246/1‑147 (MQ=255)
                                                                                                                                       agcaTGTGCCGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCgg            <  1:694961/149‑1 (MQ=255)
                                                                                                                                            gagaCGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAg       >  1:265612/1‑149 (MQ=255)
                                                                                                                                             tgaCGAAACCCGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAgg      <  2:545493/148‑1 (MQ=255)
                                                                                                                                                 gAAACACGCATGCAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGATGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGc  >  1:474742/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGC  >  NC_000913/525440‑525733

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: