Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 92,967 G→A 5.1% A519T (GCG→ACG)  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,9670GA 5.1% 518.4 / 10.0 156A519T (GCG→ACG) ftsItranspeptidase involved in septal peptidoglycan synthesis; penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base G (74/74);  new base A (3/5);  total (77/79)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 161 total aligned reads displayed.

TATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGC  >  NC_000913/92819‑93112
                                                                                                                                                    |                                                                                                                                                 
tatGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTg                                                                                                                                                   >  2:1664286/1‑149 (MQ=255)
 atgatgGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAACCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTgc                                                                                                                                                  <  1:596241/149‑1 (MQ=255)
 atgatgGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTgc                                                                                                                                                  >  3:1508253/1‑149 (MQ=255)
     aggAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCt                                                                                                                                              >  1:528185/2‑149 (MQ=255)
       gAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGc                                                                                                                                            >  4:256217/1‑149 (MQ=255)
         aaGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGa                                                                                                                                          <  2:1577980/149‑1 (MQ=255)
              tGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGCAACGCCTGCggtataa                                                                                                                                     <  2:934009/149‑7 (MQ=255)
               ggCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGc                                                                                                                                    >  4:783064/1‑149 (MQ=255)
                      ccAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCtt                                                                                                                             >  3:771130/1‑149 (MQ=255)
                       cAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTc                                                                                                                            >  4:642608/1‑149 (MQ=255)
                        aGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTcg                                                                                                                           >  1:1530198/1‑149 (MQ=255)
                        aGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTcg                                                                                                                           >  4:1177041/1‑149 (MQ=255)
                            ggcggcggcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATGAATAAATCTATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCt                                                                                                                       >  1:425757/1‑149 (MQ=255)
                              cggcggcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTgg                                                                                                                     <  3:1177041/149‑1 (MQ=255)
                               ggcggcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgt                                                                                                                    <  4:586702/149‑1 (MQ=255)
                                   gcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                <  3:1223297/149‑1 (MQ=255)
                                   gcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                <  4:623402/149‑1 (MQ=255)
                                   gcgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgtt                                                                                                                <  3:642608/149‑1 (MQ=255)
                                    cgTGAAGGCGGCGATTAAAGGTTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttg                                                                                                               <  3:420158/149‑1 (MQ=255)
                                    cgTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttg                                                                                                               <  3:783064/149‑1 (MQ=255)
                                     gTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgt                                                                                                              >  3:352245/1‑149 (MQ=255)
                                      tGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgtt                                                                                                             <  4:277033/149‑1 (MQ=255)
                                      tGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgtt                                                                                                             <  3:468578/149‑1 (MQ=255)
                                         aGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATc                                                                                                          <  1:1078236/149‑1 (MQ=255)
                                          ggcggcGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCa                                                                                                         >  3:158516/1‑149 (MQ=255)
                                           gcggcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCaa                                                                                                        <  3:1183090/149‑1 (MQ=255)
                                           gcggcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCaa                                                                                                        <  1:717412/149‑1 (MQ=255)
                                            cggcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGTGCTGGTTGTTGTTATCAAc                                                                                                       <  2:1097727/149‑1 (MQ=255)
                                             ggcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCACAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACg                                                                                                      >  1:1326179/1‑149 (MQ=255)
                                              gcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGa                                                                                                     >  1:868896/1‑149 (MQ=255)
                                              gcgATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGa                                                                                                     >  1:573590/1‑149 (MQ=255)
                                                gATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATc                                                                                                   >  1:1683907/1‑149 (MQ=255)
                                                 aTTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATcc                                                                                                  <  1:1711889/149‑1 (MQ=255)
                                                 aTTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATcc                                                                                                  <  1:1395729/149‑1 (MQ=255)
                                                   tAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGc                                                                                                >  4:1504734/1‑149 (MQ=255)
                                                       ggCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGc                                                                                            >  4:118195/1‑149 (MQ=255)
                                                       ggCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGc                                                                                            <  4:1417050/149‑1 (MQ=255)
                                                       ggCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGc                                                                                            >  4:777562/1‑149 (MQ=255)
                                                       ggCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGc                                                                                            <  2:338993/149‑1 (MQ=255)
                                                              gtatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTaaa                                                                                     <  3:601235/149‑1 (MQ=255)
                                                               tatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAt                                                                                    <  1:1557380/149‑1 (MQ=255)
                                                               tatcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAt                                                                                    <  2:177437/149‑1 (MQ=255)
                                                                atcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAta                                                                                   <  1:1306737/149‑1 (MQ=255)
                                                                atcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAta                                                                                   <  4:352245/149‑1 (MQ=255)
                                                                 tcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtac                                                                                  >  2:1550890/1‑149 (MQ=255)
                                                                 tcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtac                                                                                  <  1:1443747/149‑1 (MQ=255)
                                                                 tcgCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtac                                                                                  >  4:301546/1‑149 (MQ=255)
                                                                   gCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtacta                                                                                <  3:212469/149‑1 (MQ=255)
                                                                    ccATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtactac                                                                               >  4:103252/1‑149 (MQ=255)
                                                                    ccATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtactac                                                                               >  4:1615978/1‑149 (MQ=255)
                                                                      aTTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcgg                                                                             >  3:708210/1‑149 (MQ=255)
                                                                       ttAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcggc                                                                            >  1:1612989/1‑149 (MQ=255)
                                                                       ttAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcggc                                                                            >  3:339201/1‑149 (MQ=255)
                                                                       ttAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcggc                                                                            >  2:457551/1‑149 (MQ=255)
                                                                        tAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcggcg                                                                           <  2:581365/149‑1 (MQ=255)
                                                                         aaaaCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTAcggcgg                                                                          >  2:716091/1‑149 (MQ=255)
                                                                          aaaCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACggcggc                                                                         <  3:1287116/149‑1 (MQ=255)
                                                                          aaaCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACggcggc                                                                         <  2:1277692/149‑1 (MQ=255)
                                                                           aaCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATc                                                                                                   >  1:494372/1‑122 (MQ=255)
                                                                           aaCCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATc                                                                                                   <  2:494372/122‑1 (MQ=255)
                                                                                   ccGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTcc                                                                >  3:409420/1‑149 (MQ=255)
                                                                                   ccGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTcc                                                                <  4:33782/149‑1 (MQ=255)
                                                                                   ccGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTcc                                                                <  3:645010/149‑1 (MQ=255)
                                                                                    cgcgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCcg                                                               >  3:786584/1‑149 (MQ=255)
                                                                                     gcgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGAt                                                                                                    >  3:1304353/1‑111 (MQ=255)
                                                                                     gcgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGAt                                                                                                    <  4:1304353/111‑1 (MQ=255)
                                                                                     gcgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCcgc                                                              >  1:1092368/1‑149 (MQ=255)
                                                                                     gcgAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCcgc                                                              >  2:582815/1‑149 (MQ=255)
                                                                                       gAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCgcgc                                                            <  3:700572/149‑1 (MQ=255)
                                                                                       gAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCgcgc                                                            >  4:1053039/1‑149 (MQ=255)
                                                                                            aGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTc                                                       <  3:1053039/149‑1 (MQ=255)
                                                                                              gTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCtt                                                     <  3:1039020/149‑1 (MQ=255)
                                                                                                cGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTg                                                   >  3:1017858/1‑149 (MQ=255)
                                                                                                   gCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTg                                                <  2:42579/149‑1 (MQ=255)
                                                                                                     cGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGcc                                              >  4:407825/1‑149 (MQ=255)
                                                                                                          ggTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCcatcat                                         <  4:708210/149‑1 (MQ=255)
                                                                                                           gTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATg                                        <  4:194619/149‑1 (MQ=255)
                                                                                                           gTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATg                                        >  3:1631491/1‑149 (MQ=255)
                                                                                                             cGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATggg                                      <  4:1329769/149‑1 (MQ=255)
                                                                                                                 aCATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGggcggc                                  <  4:1498819/149‑1 (MQ=255)
                                                                                                                   aTCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGt                                >  4:378469/1‑149 (MQ=255)
                                                                                                                    tCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTa                               >  4:488865/1‑149 (MQ=255)
                                                                                                                       ataaataTATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTg                            >  4:55296/1‑149 (MQ=255)
                                                                                                                       ataaataTATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTg                            >  1:1605773/1‑149 (MQ=255)
                                                                                                                          aataTATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATGCGTa                        <  2:1293246/149‑1 (MQ=255)
                                                                                                                                ttGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATg                   >  3:998982/1‑149 (MQ=255)
                                                                                                                                  gCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGaa                 <  1:582815/149‑1 (MQ=255)
                                                                                                                                    ttATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACa               >  1:1157374/1‑149 (MQ=255)
                                                                                                                                    ttATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACa               >  2:212045/1‑149 (MQ=255)
                                                                                                                                      gactcGCAGGCGCAACGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATc             >  1:1256508/5‑149 (MQ=255)
                                                                                                                                        ctcGCAGGCGCAACGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGcc                                              <  1:229944/112‑1 (MQ=255)
                                                                                                                                        aCCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGa           <  4:849603/149‑1 (MQ=255)
                                                                                                                                        aCCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGa           <  3:232637/149‑1 (MQ=255)
                                                                                                                                        aCCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGa           <  3:454662/149‑1 (MQ=255)
                                                                                                                                         tcGCAGGCGCAACGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGggcggc                                  >  3:498002/2‑125 (MQ=255)
                                                                                                                                            cAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCgcgc                                                            <  1:405978/96‑1 (MQ=255)
                                                                                                                                            cAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCgcgc                                                            >  2:405978/1‑96 (MQ=255)
                                                                                                                                            cAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCg       >  3:212708/1‑149 (MQ=255)
                                                                                                                                               gCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGAt    <  2:566121/149‑1 (MQ=255)
                                                                                                                                               gCGTTGCGCCTGCGAGTCAGCCGCACTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGAt    <  4:196688/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGC  >  NC_000913/92819‑93112

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