Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,128,392 C→T P185L (CCG→CTG)  murJ → lipid II flippase MurJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,128,3920CT98.6% 229.4 / ‑5.7 70P185L (CCG→CTG) murJlipid II flippase MurJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/39);  minor base G (0/1);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

GCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACCGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTAT  >  NC_000913/1128151‑1128636
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
gccgcACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTgctggc                                                                                                                                                                                                                                             <  2:306803/251‑1 (MQ=255)
    cACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACCAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGAACTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGGTTCCGGCGTTTGCTCCAACACTGCTTAACATCAGTATGATTGTTTTCGCGCTGTTTGCCGCACCGTACTTTATCCCACGGGTgctggcgctg                                                                                                                                                                                                                                         <  2:375737/251‑1 (MQ=255)
                gcatagTGACCGCGCAAGGCTTCACTGACACAGCTGACAAATTTGCCCTGGCCAGCCAGCTACTAAAGATTGCCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTTGTGGGAGCGAGTCTGAATAGGTGGAGCCGCTTCTCGATGACGGCGTGTGCTCCAACACTGCGTAACATCAGCATGATTGGTTTCGCGCGGTTTTCAACGCCGTTCTTTAACCCGCTGGTgcttgcgatgccgttgtctgtt                                                                                                                                                                                                                             <  2:406415/246‑16 (MQ=255)
                         ccGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTcggc                                                                                                                                                                                                                    <  1:352517/251‑1 (MQ=255)
                             gcCAGGCTTCGCTGACACAGCTGACAAATTCGCCCTGATCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGTTGCTGGCGCTGGTGTGGGCTGTTACGGTCGgcggcg                                                                                                                                                                                                                >  1:474914/1‑251 (MQ=255)
                                  gCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTg                                                                                                                                                                                                           >  2:478559/1‑251 (MQ=255)
                                   cTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGc                                                                                                                                                                                                          >  1:469369/1‑251 (MQ=255)
                                       gCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAgtt                                                                                                                                                                                                      >  2:564061/1‑249 (MQ=255)
                                        cTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTg                                                                                                                                                                                                     >  2:261485/1‑251 (MQ=255)
                                          gACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGgt                                                                                                                                                                                                   <  1:399112/251‑1 (MQ=255)
                                          gACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGgt                                                                                                                                                                                                   <  1:478559/251‑1 (MQ=255)
                                          gACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGgt                                                                                                                                                                                                   <  2:125670/251‑1 (MQ=255)
                                           acacAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGgtg                                                                                                                                                                                                  >  1:385364/1‑251 (MQ=255)
                                              caGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTcggcg                                                                                                                                                                                                                   <  2:528184/231‑1 (MQ=255)
                                              caGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTcggcg                                                                                                                                                                                                                   >  1:528184/1‑231 (MQ=255)
                                               aGCTGACAAATTTGACCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGCATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTACCATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTCCGGTCGGCGGCGGCATGCAGCTGGTGNATc                                                                                                                                                                                              >  2:607720/1‑251 (MQ=255)
                                                            gCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCGTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGGTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCAcc                                                                                                                                                                                  <  1:34519/250‑1 (MQ=255)
                                                             ccgtGACCAGCCAGCGGCTATGGATTACCTTAGCCTACTTCTTGCTGATCTACCTGGCGTCGCTGGGGGGATCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCATCACTGCTTAACATCAGCATGATTTGTTTCGAGCTGTTCGCCGCACCGTTCTTTAACCCACTGCTGCTGGCTCTGCCGTTGGTTGTTTCGGTTGGCTGCGTCCTGCAGCTGGTGTTTCATGTACCGCACCtt                                                                                                                                                                                <  2:375048/248‑2 (MQ=255)
                                                                 gACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTgaa                                                                                                                                                                              <  2:172705/249‑1 (MQ=255)
                                                                 gACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTga                                                                                                                                                                               >  2:309481/1‑248 (MQ=255)
                                                                  aCCAGCCAGCTACTAACGATTACCTTACCCTATATCTTGCTGATCTCCCTTGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGCTTCCGGCGTTTGCTCCAACACTGCTTAACCTCACCATGATTGGTTTCGCGCTGTTTGCCGCACCGCACTTTAACCCACTCGTGATGGCGCTGGCGTGGGCTGTTCCGGTCGGCGGCGTCCTGCAGCTGGTGTATCCGCTACCGCACCTgaata                                                                                                                                                                            >  2:64215/1‑248 (MQ=255)
                                                                   ccagccagCTACTAAAGATTACCTTTCCCTAGATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGaagaag                                                                                                                                                                          <  1:452121/251‑1 (MQ=255)
                                                                       ccagcTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTGAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCg                                                                                                                                                                      <  2:307781/251‑1 (MQ=255)
                                                                              ctaAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCt                                                                                                                                                               >  1:131655/1‑251 (MQ=255)
                                                                                 aaGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGt                                                                                                                                                            <  1:353469/251‑1 (MQ=255)
                                                                                    aTTACCTTTCACTATATCTTGCTGATCTACCTGGCGCCGCTGGTGCGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTCCGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGGCCt                                                                                                                                                         >  2:113176/1‑251 (MQ=255)
                                                                                       aCCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACTTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGcc                                                                                                                                                      <  1:30678/251‑1 (MQ=255)
                                                                                        ccTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCcg                                                                                                                                                     >  1:641297/1‑251 (MQ=255)
                                                                                           ttCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCgcgc                                                                                                                                                  >  1:374973/1‑251 (MQ=255)
                                                                                           ttCCCGATATCTTGCTGATCTCCCTGGCGTTGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGTTTTCGCGCTGTTTGCCGCACCGTTCTTAAACCCACTGGTGCTGGCGATGGCGTGGGCTGTTACGATCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCcgc                                                                                                                                                    <  2:433247/249‑1 (MQ=255)
                                                                                            tCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTcc                                                                                                                                                          <  1:594942/242‑1 (MQ=255)
                                                                                            tCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTcc                                                                                                                                                          >  2:594942/1‑242 (MQ=255)
                                                                                              ccTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCAtt                                                                                                                                               <  1:524948/251‑1 (MQ=255)
                                                                                              ccTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCAtt                                                                                                                                               <  2:457326/251‑1 (MQ=255)
                                                                                                tataTCTTGCTGATCTCCCTGGCGTAGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGCTTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGACCTGCCGCGCATTa                                                                                                                                              >  2:307673/1‑250 (MQ=255)
                                                                                                  tatCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACt                                                                                                                                           >  1:218369/1‑251 (MQ=255)
                                                                                                              ctcCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCg                                                                                                                               >  1:553562/1‑251 (MQ=255)
                                                                                                                cccTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATg                                                                                                                                  <  2:469369/246‑1 (MQ=255)
                                                                                                                cccTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATg                                                                                                                                  <  1:657079/246‑1 (MQ=255)
                                                                                                                   tGGCGTCGCTGGTGGGTGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCCGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGc                                                                                                                           >  2:463100/1‑250 (MQ=255)
                                                                                                                    ggCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGGGTTTGCTTCATCACTTCTTAACATCAGCATGATTTGTTGCGCGCTGTTTGCCGCACCGTACTTTACCCCACTGGTGCTGGCGCTGGCGTGTGCTTTTACGGTCGGCGGCTCCCTGCAGCTGGTGTATCAGCTACCTCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCa                                                                                                                          <  1:325261/250‑1 (MQ=255)
                                                                                                                             ggtggGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGAGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCgtgg                                                                                                                >  1:578824/1‑251 (MQ=255)
                                                                                                                              gtggGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCgtggt                                                                                                               <  2:218369/251‑1 (MQ=255)
                                                                                                                                 ggAGCGATTCTTAATTCGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCTCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGaa                                                                                                            <  1:99167/251‑1 (MQ=255)
                                                                                                                                  gAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCCCCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGaaa                                                                                                           <  2:387389/251‑1 (MQ=255)
                                                                                                                                    gCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGGGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGGGGTGAAAc                                                                                                          >  2:637406/1‑250 (MQ=255)
                                                                                                                                          cTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCgt                                                                                                                  >  1:611093/1‑236 (MQ=255)
                                                                                                                                          cTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCgt                                                                                                                  <  2:611093/236‑1 (MQ=255)
                                                                                                                                            gAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGAc                                                                                                 <  1:660703/251‑1 (MQ=255)
                                                                                                                                             aaTACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGAcc                                                                                                >  2:324739/1‑251 (MQ=255)
                                                                                                                                                      aaCCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCAGGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGGGCTGGCGCTGGCGTTGGCTGTTGTGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACGTCCACGATGCCGGATCAATGCGCGTGGTGAAACATATGGGACCGGCGATCCt                                                                                       <  2:623719/251‑1 (MQ=255)
                                                                                                                                                       aCCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCtt                                                                                      <  1:492585/251‑1 (MQ=255)
                                                                                                                                                          gCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGAGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGCTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGc                                                                                   >  1:614100/1‑251 (MQ=255)
                                                                                                                                                                cGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCTCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGtct                                                                               <  2:353885/249‑1 (MQ=255)
                                                                                                                                                                       ggCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGCTGGGACCGGCGATCCTTGGCGTCTCTGTGAGcc                                                                      >  2:78244/1‑251 (MQ=255)
                                                                                                                                                                             tGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGAtct                                                                >  2:498545/1‑251 (MQ=255)
                                                                                                                                                                              gCNCCAACACTGCTTACCATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGGCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCTCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGAtct                                                                <  2:596672/250‑1 (MQ=255)
                                                                                                                                                                                 ccAACACTGCTTAACATCAGCATGAGTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCtt                                                            <  1:604865/251‑1 (MQ=255)
                                                                                                                                                                                   aaCACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTaa                                                          <  1:261485/251‑1 (MQ=255)
                                                                                                                                                                                            tAACATCAGCAGGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCTCCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGTCCGGCGATCCTTGGCGTCTCTGTGATCCAGATCTCCTTAATCATCaa                                                   <  2:317638/249‑1 (MQ=255)
                                                                                                                                                                                             aaCAGCAGCACGATTTGTTTCGCGCGGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCTCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCAGCAGCGGGTGTATCGGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCGTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAacac                                                <  2:544904/251‑1 (MQ=255)
                                                                                                                                                                                                  cAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttt                                           >  1:151086/1‑251 (MQ=255)
                                                                                                                                                                                                        gATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCt                                     <  1:58657/251‑1 (MQ=255)
                                                                                                                                                                                                            ggTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGtt                                 <  1:309481/251‑1 (MQ=255)
                                                                                                                                                                                                               ttCGCGCTGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGAACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCt                              >  1:78735/1‑251 (MQ=255)
                                                                                                                                                                                                                      tGTTTGCCGCACCGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATACACACCATTTTTGACTCGTTTCTTGCTTcc                       >  1:562609/1‑251 (MQ=255)
                                                                                                                                                                                                                           gCCGCACCGTACTTTAACCCACTGGTGCTTGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTc                  <  2:319285/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTc            >  2:424261/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gTACTTTAACCCACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCtt          <  1:651548/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ccccACTGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTAt  <  2:502881/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
GCCGCACCGTGGGTGATCATGGTGACCGCGCCAGGCTTCGCTGACACAGCTGACAAATTTGCCCTGACCAGCCAGCTACTAAAGATTACCTTTCCCTATATCTTGCTGATCTCCCTGGCGTCGCTGGTGGGAGCGATTCTGAATACGTGGAACCGCTTCTCGATTCCGGCGTTTGCTCCAACACTGCTTAACATCAGCATGATTGGTTTCGCGCTGTTTGCCGCACCGTACTTTAACCCACCGGTGCTGGCGCTGGCGTGGGCTGTTACGGTCGGCGGCGTCCTGCAGCTGGTGTATCAGCTACCGCACCTGAAGAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTAT  >  NC_000913/1128151‑1128636

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