Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 93,173 A→C 100% R587S (AGA→AGC
D3A (GAT→GCT) 
ftsI →
murE →
transpeptidase involved in septal peptidoglycan synthesis; penicillin‑binding protein 3
UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091393,1730AC98.9% 272.1 / ‑5.1 92R587S (AGA→AGC
D3A (GAT→GCT) 
ftsI
murE
transpeptidase involved in septal peptidoglycan synthesis; penicillin‑binding protein 3
UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate:meso‑ diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (29/61);  minor base T (0/1);  total (30/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGATCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCGATAT  >  NC_000913/93032‑93318
                                                                                                                                             |                                                                                                                                                 
ctacGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAAt                                                                                                                                            >  2:1274426/1‑149 (MQ=255)
 tacGGCGGCGCTGTTTCCGCGCCGGCCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACATCGGGCGATAAAAATGAATTTGTGATTAATCAGGGCGAGGGGACAGGTGGCAGCTCGTAAtt                                                                                                                                           <  1:795245/149‑1 (MQ=255)
  acGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTCAttt                                                                                                                                          >  1:1244446/1‑149 (MQ=255)
  acGGCGGCGCCGCTTCCGCGCCGGTCTTTGAACGCGGCAGACGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAAttt                                                                                                                                          >  1:1283739/1‑149 (MQ=255)
    ggcggcGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCGGCCCGTAATTTgc                                                                                                                                        >  2:251728/1‑149 (MQ=255)
    ggcggcGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGACTAATCAAGGCGAGGGGACAGGTGGCAGCCCGTAATTTgc                                                                                                                                        >  2:666462/1‑149 (MQ=255)
     gcggcgCCGTTTCCGCGCCGGTCTTTGGTGCCATCCTGGGCGGCGTATTGCGGACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAATGCGAGGGGACAGGTGGCAGCTCGTAATTTgcg                                                                                                                                       <  2:544287/149‑1 (MQ=255)
        gcgcCGTTTCCGCGCCGGTCTTTGGTGCGATCATGGGCGGCGTATGGCGTACCATGAACATCGAGCCGGACGCGCTGACAGCGGGCGATAAAAATGAATTTGGGGGTAATCAACGCGAGGGGACAGGTGGCAGCTGGGAGTTTGCGGGa                                                                                                                                    >  1:694090/1‑149 (MQ=255)
              tttCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACcttct                                                                                                                              <  1:970793/149‑1 (MQ=255)
               ttCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACcttctg                                                                                                                             >  2:443010/1‑148 (MQ=255)
                tCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTg                                                                                                                            <  1:340909/149‑1 (MQ=255)
                   gcgcCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCGGGTGGTAATTTGGGCGACCTTCTTGCtt                                                                                                                         >  1:14568/1‑148 (MQ=255)
                   gcgcCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTc                                                                                                                         >  1:1122850/1‑149 (MQ=255)
                      ccGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGCGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGt                                                                                                                      <  2:1149808/149‑1 (MQ=255)
                       tggTTTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTg                                                                                                                     <  2:196994/148‑1 (MQ=255)
                        gggttgtGGGGCCAGCATGGGCGGCGTATTGCGTACCATGAACAGCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTgg                                                                                                                    <  2:185625/143‑1 (MQ=255)
                        ggTTTTTGGTGCCATCAGGGGCGGCGGATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTgg                                                                                                                    <  2:687536/149‑1 (MQ=255)
                        ggTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTgg                                                                                                                    <  2:75495/149‑1 (MQ=255)
                            tttGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGCGGGGACAGGTGGCAGCTCGTAGTTTGCGCGACCTTCCTGCTCCGTGGGTGc                                                                                                                >  2:302837/1‑149 (MQ=255)
                                 tGCCATCATGGGCGTCGTATTGCGGACCAGGAACATCGAGCCGGATGCGCTGACAACGGGCGAGAAAAATGTATTTGTGATTAATCAAGGCGTGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGAc                                                                                                           <  2:438960/149‑1 (MQ=255)
                                 tGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGAc                                                                                                           <  1:663654/149‑1 (MQ=255)
                                 ggCCATCATGTGCCGCGTATTGCGTACCACGAACACCCGGCCGCATGCGCTGACACCGGGCGATAACAATGAATTTGTGCTTAATCATGGCGTGGGGACAGGTGGCCGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGAc                                                                                                           <  2:112465/148‑1 (MQ=255)
                                  gCCATCATGGGCGGCGTAGTGCGTACCCTGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACg                                                                                                          <  2:1165250/149‑1 (MQ=255)
                                     atcatGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGGCCTTCTTGCTCCGTGGGTGCCAGACGCAc                                                                                                       >  1:1232201/1‑149 (MQ=255)
                                     atcatGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCAc                                                                                                       >  1:516591/1‑149 (MQ=255)
                                      tcatGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCAcc                                                                                                      >  1:702647/1‑149 (MQ=255)
                                      tcatAGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTGGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAGTTTGCGCGACCTTCTTGCCCCCTGGGTGCCAGCCGGAcc                                                                                                      >  2:666565/1‑149 (MQ=255)
                                         tGGGGGGCGTATTGCGTACCATGAACATCGAGCCGGAGGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTATTTTGCGCCACCTCCTTGCTCCGTGGGTGCCAGACGCACCTTc                                                                                                   <  1:955751/149‑1 (MQ=255)
                                         tGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTc                                                                                                   <  1:1025832/149‑1 (MQ=255)
                                           ggcggcGTACTGCGTCCCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTcgc                                                                                                 <  1:1077981/149‑1 (MQ=255)
                                              ggcgTATTGCGTACCATGAACATCGAGTCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAg                                                                                              <  1:904051/149‑1 (MQ=255)
                                              ggcgTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAg                                                                                              <  1:1269098/149‑1 (MQ=255)
                                              ggcgCATTGCGTGCCAGGAACATCGAGCCGGATGCGCTGACAACGGGCGCTAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCGCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAg                                                                                              <  1:909714/149‑1 (MQ=255)
                                              aacgTATTGCGTACCTTGAACAGCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGGATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAg                                                                                              <  2:322620/147‑1 (MQ=255)
                                                  tATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACt                                                                                          >  2:1076157/1‑149 (MQ=255)
                                                  tATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACt                                                                                          >  2:1270220/1‑149 (MQ=255)
                                                  tATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGGGCCAGACGCACCTTCGCGAGCACt                                                                                          >  2:337909/1‑149 (MQ=255)
                                                   aTTGCGTACCATGAGCATCGAGCCGAATGCGCTTACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACCGGTGGAAGCTCGTCATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTg                                                                                         <  1:792204/149‑1 (MQ=255)
                                                   aTTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGTATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCACCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACTCACCTTCGCGAGCACTg                                                                                         <  2:957226/149‑1 (MQ=255)
                                                    ttGCTTACCATGAACATCGCGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGc                                                                                        <  1:1154360/149‑1 (MQ=255)
                                                       cGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCgag                                                                                     <  2:1062664/149‑1 (MQ=255)
                                                       cGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCgag                                                                                     <  2:988436/149‑1 (MQ=255)
                                                       cGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCgag                                                                                     <  1:402597/149‑1 (MQ=255)
                                                       cGTACCAATAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGCAAGTTCGTTATTTGCGAAAACTTCGGGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCgag                                                                                     <  1:309570/149‑1 (MQ=255)
                                                           ccATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGGGTCAGACGCACCTTCGCGAGCACTTCGGAAGAt                                                                                 >  1:1162486/1‑149 (MQ=255)
                                                               gAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGaca                                                                             >  2:295881/1‑149 (MQ=255)
                                                                 aCATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACt                                                                           <  1:841340/149‑1 (MQ=255)
                                                                 aCAACGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAACTTGCGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACt                                                                           <  1:593798/149‑1 (MQ=255)
                                                                    tCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGGGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGa                                                                        >  2:721638/1‑149 (MQ=255)
                                                                       aGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAg                                                                     >  1:868977/1‑149 (MQ=255)
                                                                          cGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGGGAGATGACACTCGACAGCCg                                                                  >  2:1094942/1‑149 (MQ=255)
                                                                           ggATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCGGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGGCGCACCTTCGCGAGGACTGCGCGAGATGTTACTCGCCAGCCgt                                                                 >  1:369764/1‑149 (MQ=255)
                                                                                ggcTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACGCGACAGCCGTGTGGc                                                            >  2:484620/2‑149 (MQ=255)
                                                                                 gcTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCATCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCt                                                           >  1:862126/1‑149 (MQ=255)
                                                                                 gcTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACCCTCGACAGCCGTGTGGCt                                                           >  1:1058842/1‑149 (MQ=255)
                                                                                  cTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGCGTGGccg                                                          >  2:321502/1‑147 (MQ=255)
                                                                                          gggCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGGGAGATGACACTCGACAGCCGTGTGGCTGCggcgggc                                                  <  2:174479/149‑1 (MQ=255)
                                                                                           ggCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGTGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGCGGCTGCggcgggcg                                                 <  1:1158584/149‑1 (MQ=255)
                                                                                             cGATAACAATGAGTGTGTGATTAAGCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGAGGACACTCGACAGCCGTGTGGCTGCGGCGGGCGAt                                               <  1:322504/149‑1 (MQ=255)
                                                                                             cGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGAt                                               <  1:1274426/149‑1 (MQ=255)
                                                                                              gATAAAATGGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGAtc                                              <  1:594031/149‑1 (MQ=255)
                                                                                                    aaTGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTg                                        <  1:387434/149‑1 (MQ=255)
                                                                                                    aaTGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGGAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTg                                        <  1:1198065/149‑1 (MQ=255)
                                                                                                     attaaTTTGTGATTAAGCAAGGCGAGGGGACAGGTGGCAGCTCATAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGt                                       <  1:190521/146‑1 (MQ=255)
                                                                                                     aTGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGt                                       <  1:890235/149‑1 (MQ=255)
                                                                                                      tGAATTGGTGATTAAGCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCCCCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTa                                      <  2:148978/149‑1 (MQ=255)
                                                                                                      ggAAGTTGTGATCAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTGGCGCGACCTTCTTGCGCCGTGGGTGCCAGAGGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTa                                      <  2:174705/148‑1 (MQ=255)
                                                                                                       gAATTTGTGACTGATCAAGGCGAGGGGGCAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGGGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAg                                     <  2:1232112/149‑1 (MQ=255)
                                                                                                        aaTTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGc                                    >  1:333943/1‑149 (MQ=255)
                                                                                                         aTTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGCCAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCt                                   <  1:838830/149‑1 (MQ=255)
                                                                                                          tttGTGATTAATCAAGGTGAGGGAGCAGGAGGCAGCTCGTAACTTGCGCGACCTTCCTGCTCCGTGGGTGCCAGACGCCCCTTCGCGAGCACTGCGAGAGATGACACCCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTg                                  <  1:443010/149‑1 (MQ=255)
                                                                                                           ttGTGATTAATCAAGGCGGGGGGACAGGTGGCAGCTCGTAATTTGCGCGACTTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTgt                                 <  2:1122850/149‑1 (MQ=255)
                                                                                                           ttGCGATTAATCAAGGCGAGGGCACAGGTGGCAGCTCGTAATTTGCGCGACCTTCCTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCATTGTAGCTgt                                 <  2:14568/149‑1 (MQ=255)
                                                                                                            tgtgATTAATCAGGGCGAGGGGACAGGTGGCAGCTCGCAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTgta                                <  2:9360/149‑1 (MQ=255)
                                                                                                            tgtgATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTgta                                <  2:1286495/149‑1 (MQ=255)
                                                                                                              tgATTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTgtagt                              <  2:333943/149‑1 (MQ=255)
                                                                                                                actAATCAAGGCGAGGGGACAGGTGGCAGCTCGCAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGGCGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGtagtag                            <  1:377415/147‑1 (MQ=255)
                                                                                                                aTTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGtagtag                            <  1:376746/149‑1 (MQ=255)
                                                                                                                aTTAATCAAGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCAGTGTTTTAGCGGCGGGCGATCTCTTTGTAGCTGtagtag                            <  1:635469/149‑1 (MQ=255)
                                                                                                                   aaTCAAGGCGAGTGGACAGGTGGCAGCCCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGtc                         <  1:352090/149‑1 (MQ=255)
                                                                                                                       aaGGCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCAACCGCCGTGTGGCTGCGGCGGGCGATCTCCTTGTGGCGGTAGTAGGtcatca                     >  2:139384/1‑149 (MQ=255)
                                                                                                                          gCGAGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGc                  <  1:1015937/149‑1 (MQ=255)
                                                                                                                             gggggACGGGGGGCAGCTCGTAGTTTGCGCGACCTCCTTGCTCCGTGGGTGCCAGACGCCCCTTCGCGCGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGa               <  2:770195/148‑1 (MQ=255)
                                                                                                                             aGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGa               <  1:1009778/149‑1 (MQ=255)
                                                                                                                             aGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGa               <  1:23012/149‑1 (MQ=255)
                                                                                                                             aGGGGACAGGTGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGACAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGa               <  1:696439/149‑1 (MQ=255)
                                                                                                                             aGGGGACAGGCGGCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCGCTTTGTAGCTGTAGTAGGTCATCAGGCGGa               <  2:125489/149‑1 (MQ=255)
                                                                                                                                        ggCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCGat    <  1:781715/149‑1 (MQ=255)
                                                                                                                                        ggCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCGat    <  1:1215675/149‑1 (MQ=255)
                                                                                                                                        ggCAGCTCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCGat    <  1:1177501/149‑1 (MQ=255)
                                                                                                                                         gCAGCTCGTAATTTGCGCGACCTTCTTCCCCCGTCGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGCCTGCGGCGGGCGCTCTCTCCGCACCTGCAGCAGGCCACCAGCCCGCCGGGCGTCGcta   >  2:171181/1‑149 (MQ=255)
                                                                                                                                          cAGCTCGTAAGTTGCGCGACCTTGTTGCTCCGTGGGGGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACGCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCGatat  <  1:908001/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
CTACGGCGGCGCCGTTTCCGCGCCGGTCTTTGGTGCCATCATGGGCGGCGTATTGCGTACCATGAACATCGAGCCGGATGCGCTGACAACGGGCGATAAAAATGAATTTGTGATTAATCAAGGCGAGGGGACAGGTGGCAGATCGTAATTTGCGCGACCTTCTTGCTCCGTGGGTGCCAGACGCACCTTCGCGAGCACTGCGAGAGATGACACTCGACAGCCGTGTGGCTGCGGCGGGCGATCTCTTTGTAGCTGTAGTAGGTCATCAGGCGGACGGGCGTCGATAT  >  NC_000913/93032‑93318

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