Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 92,966 T→A 5.6% V518V (GTT→GTA) ‡ ftsI → transpeptidase involved in septal peptidoglycan synthesis; penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091392,9660TA 5.6% 304.3 / 3.8 94V518V (GTT→GTA) ‡ftsItranspeptidase involved in septal peptidoglycan synthesis; penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base T (48/40);  new base A (4/2);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCG  >  NC_000913/92819‑93113
                                                                                                                                                   |                                                                                                                                                   
tatGATGGAAAGCGTGGCGCTACCCGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTg                                                                                                                                                    >  1:242946/1‑149 (MQ=255)
   gatgGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGAATAAAGGCTATCGTATCGCCATTAAAAACGGTACCGCGAAAAAGGTAGGCCCGGACGGTCGCTACATCAAAAAATATATAGCTTATACCGCAGGCGTTgcgc                                                                                                                                                 >  2:261442/1‑149 (MQ=255)
         aaGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGc                                                                                                                                             >  6:48293/1‑147 (MQ=255)
         aaGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGa                                                                                                                                           <  8:293077/149‑1 (MQ=255)
         aaGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGGCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGa                                                                                                                                           <  7:361793/149‑1 (MQ=255)
            cGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTc                                                                                                                                        <  8:65199/149‑1 (MQ=255)
                  gcTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCcgc                                                                                                                                  >  4:169502/1‑149 (MQ=255)
                            ggcggcggcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCt                                                                                                                        <  2:242946/149‑1 (MQ=255)
                             gcggcggcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTg                                                                                                                       >  8:372660/1‑149 (MQ=255)
                             gcggcggcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTg                                                                                                                       <  8:356600/149‑1 (MQ=255)
                                   gcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                 <  8:63262/149‑1 (MQ=255)
                                    cgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAACAATGTCGGGCCGGACGGGCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttg                                                                                                                <  2:417382/149‑1 (MQ=255)
                                     gTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgt                                                                                                               >  7:84336/1‑149 (MQ=255)
                                      tGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgtt                                                                                                              <  3:169502/149‑1 (MQ=255)
                                       gAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTa                                                                                                             >  5:269161/1‑149 (MQ=255)
                                       gAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTa                                                                                                             >  7:173912/1‑149 (MQ=255)
                                        aaGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTAt                                                                                                            <  8:96445/149‑1 (MQ=255)
                                         aGGCGGCGATTAAATGCTATCGTATCTCCATTAAAACCGGTCCCGCGAAAAATGGCGGGCCGGAAGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATc                                                                                                           <  6:269161/149‑1 (MQ=255)
                                         aGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCcgcg                                                                                                                                 >  1:149981/1‑127 (MQ=255)
                                         aGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCcgcg                                                                                                                                 <  2:149981/127‑1 (MQ=255)
                                             ggcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAAc                                                                                                        >  7:68224/1‑148 (MQ=255)
                                                   tAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGATTATACAGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCg                                                                                                  >  7:72691/1‑148 (MQ=255)
                                                    aaaGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCg                                                                                                  >  1:378361/1‑147 (MQ=255)
                                                         cTATCGTATCGCCATTAAAACCGGTACCGCGAAATAGGTCGGGCCGGACGGGCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCgg                                                                                           <  1:261442/149‑1 (MQ=255)
                                                          tatcgtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTc                                                                                                                             <  5:146025/114‑1 (MQ=255)
                                                          tatcgtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTc                                                                                                                             >  6:146025/1‑114 (MQ=255)
                                                           atcgtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCggg                                                                                          >  2:463147/1‑148 (MQ=255)
                                                              gtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGGGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTaaa                                                                                      >  6:53044/1‑149 (MQ=255)
                                                              gtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGc                                                                                             <  7:261405/142‑1 (MQ=255)
                                                              gtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGc                                                                                             >  8:261405/1‑142 (MQ=255)
                                                               tatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGGCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAt                                                                                     <  6:245886/149‑1 (MQ=255)
                                                                atcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATcc                                                                                                   >  2:304020/1‑134 (MQ=255)
                                                                atcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATcc                                                                                                   <  1:304020/134‑1 (MQ=255)
                                                                 tcgCCATTAAAACCGGTACCTCGAAAAAGGTCGGGCCGGACTGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtac                                                                                   <  6:268424/149‑1 (MQ=255)
                                                                   gCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCa                                                                                                                                       >  4:211120/1‑95 (MQ=255)
                                                                   gCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCa                                                                                                                                       <  3:211120/95‑1 (MQ=255)
                                                                   gCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTAATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtacta                                                                                 <  2:125695/149‑1 (MQ=255)
                                                                    ccATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtactac                                                                                >  2:339576/1‑149 (MQ=255)
                                                                    ccATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtactac                                                                                <  8:84336/149‑1 (MQ=255)
                                                                     cATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcg                                                                               >  7:364368/1‑149 (MQ=255)
                                                                      aTTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcgg                                                                              <  7:372660/149‑1 (MQ=255)
                                                                      aTTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTAGCGCCTGCGAGTCAGCCTCGCTTCGCTCTGTTTTTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcgg                                                                              >  5:322627/1‑149 (MQ=255)
                                                                          aaaCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACggcggc                                                                          <  8:68224/149‑1 (MQ=255)
                                                                             ccGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCaa                                                                                                         >  7:76180/1‑115 (MQ=255)
                                                                             ccGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCaa                                                                                                         <  8:76180/115‑1 (MQ=255)
                                                                                   ccGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTcc                                                                 >  3:229530/1‑149 (MQ=255)
                                                                                    cgcgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                 >  4:422551/1‑100 (MQ=255)
                                                                                    cgcgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                 <  3:422551/100‑1 (MQ=255)
                                                                                     gcgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCcgc                                                               >  2:388435/1‑149 (MQ=255)
                                                                                      cgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgtt                                                                                                              >  4:364354/1‑101 (MQ=255)
                                                                                      cgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgtt                                                                                                              <  3:364354/101‑1 (MQ=255)
                                                                                       gAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGa                                                                                                                                           >  7:278677/1‑71 (MQ=255)
                                                                                       gAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGa                                                                                                                                           <  8:278677/71‑1 (MQ=255)
                                                                                       gAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGa                                                                                                      >  6:55538/1‑108 (MQ=255)
                                                                                       gAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGa                                                                                                      <  5:55538/108‑1 (MQ=255)
                                                                                          aaaGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCgg                                                          >  7:280225/1‑149 (MQ=255)
                                                                                            aGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGGGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTc                                                        <  5:53044/149‑1 (MQ=255)
                                                                                                 gggCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTg                                                    <  8:173912/148‑1 (MQ=255)
                                                                                                  ggCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGt                                                  <  1:339576/149‑1 (MQ=255)
                                                                                                      ggACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCca                                              <  5:48293/149‑1 (MQ=255)
                                                                                                       gACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCcat                                             >  5:93162/1‑149 (MQ=255)
                                                                                                           gTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATg                                         <  2:378361/149‑1 (MQ=255)
                                                                                                            tCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGGTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATgg                                        <  4:255389/149‑1 (MQ=255)
                                                                                                                taaaTCAATAAATATATTGCTTATACCGCAGGCGTTGCGCATGCGAGTCAGCCGCGCTTCGCTCTGGTTGTTGTTATCTACGATCCGCAGGCGGGTAAATACTACGGAGGCGCCGTTTCCGCGCCGGTGTTTGGTGCGATCATGggcgg                                    >  6:89403/4‑149 (MQ=255)
                                                                                                                   aTCAATAAATATATTTCTTATACCGCATGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGAGGCGCCGTTTCAGCGACGGTCTTTGGTGCCATCATGGgcggcg                                  >  5:396474/1‑148 (MQ=255)
                                                                                                                     cAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTAt                               >  1:385087/1‑149 (MQ=255)
                                                                                                                      aattaataTATTGCTTATAGCGCAGGAGTTGCGCCTGCGAGTCAGCCTGGCTTCGCGCTGGTTTTTGTTATCAACGCTCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTAtt                              <  6:322627/149‑1 (MQ=255)
                                                                                                                       ataaataTATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTg                             >  2:414067/1‑149 (MQ=255)
                                                                                                                         aaataTATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAATTACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCg                           >  7:444045/1‑149 (MQ=255)
                                                                                                                           atatatTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTa                         >  2:99648/1‑149 (MQ=255)
                                                                                                                            tatatTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTAc                        >  1:59939/1‑149 (MQ=255)
                                                                                                                            tatatTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTAc                        >  5:2234/1‑149 (MQ=255)
                                                                                                                               atTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCAt                     >  6:109890/1‑149 (MQ=255)
                                                                                                                                ttGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCAt                     >  4:367129/1‑148 (MQ=255)
                                                                                                                                ttGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCAt                     <  1:463147/148‑1 (MQ=255)
                                                                                                                                    ttATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGggcggc                                   <  6:270848/130‑1 (MQ=255)
                                                                                                                                    ttATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGggcggc                                   >  5:270848/1‑130 (MQ=255)
                                                                                                                                    ttATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACa                >  1:331170/1‑149 (MQ=255)
                                                                                                                                      ataCCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGGCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACAt               <  4:229530/148‑1 (MQ=255)
                                                                                                                                       taCCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATc              <  3:367129/148‑1 (MQ=255)
                                                                                                                                       taCCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATc              >  8:311164/1‑148 (MQ=255)
                                                                                                                                        cccGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGa            >  8:393181/2‑149 (MQ=255)
                                                                                                                                         tcGCAGGCGCAACGCCTGCGGGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtact                                                                                  >  3:224186/2‑78 (MQ=38)
                                                                                                                                         tcGCAGGCGCAACGCCTGCGGGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtact                                                                                  <  4:224186/77‑1 (MQ=38)
                                                                                                                                         tcGCAGGCGCAACGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAg                                                                                               >  2:270266/2‑65 (MQ=38)
                                                                                                                                         tcGCAGGCGCAACGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAg                                                                                               <  1:270266/64‑1 (MQ=38)
                                                                                                                                         tcGCAGGCGCAACGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAg           >  4:233295/2‑149 (MQ=255)
                                                                                                                                          cGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGc          >  6:238791/1‑149 (MQ=255)
                                                                                                                                           gCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGcc         <  1:414067/149‑1 (MQ=255)
                                                                                                                                           gCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGcc         >  3:324687/1‑149 (MQ=255)
                                                                                                                                           gCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCACGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGcc         >  6:66917/1‑149 (MQ=255)
                                                                                                                                              ggCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGa      >  6:60642/1‑149 (MQ=255)
                                                                                                                                               gCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCa                      >  2:379152/1‑132 (MQ=255)
                                                                                                                                               gCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCa                      <  1:379152/132‑1 (MQ=255)
                                                                                                                                                  ttGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATgcg  >  2:28667/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCG  >  NC_000913/92819‑93113

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