Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,469,180 A→T 5.8% L273H (CTC→CAC)  mrdB ← murein transpeptidase factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,469,1800AT 5.8% 247.2 / 2.6 80L273H (CTC→CAC) mrdBmurein transpeptidase factor
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base T (2/3);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CAAAAAAAGTCATGGTCAGGCTGCCCGCAAGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGTGCCCAACGG  >  NC_002947/5469037‑5469322
                                                                                                                                                |                                                                                                                                              
cAAAAAAAGTCATGGTCAGGCTGCCCGCAAGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACg                                                                                                                                           >  5:230002/1‑149 (MQ=255)
 aaaaaaaGTCATGGTCAGGCTGCCCGCAAGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACgg                                                                                                                                          >  2:181384/1‑149 (MQ=255)
        gTCATGGTCAGGCTGCCCGCAAGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGTGCTCGGCgaggaa                                                                                                                                   <  5:217871/149‑2 (MQ=255)
              gTCAGGCTGCCCGCAAGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTc                                                                                                                             >  7:92936/1‑149 (MQ=255)
                    cTGCCCGCAAGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCCGCCGAGCACGGCGATGATGAAGTCGgtgt                                                                                                                        <  1:187932/149‑1 (MQ=255)
                        ccGCAAGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGTGCTCGGCgagg                                                                                                                                     >  1:1304/1‑129 (MQ=255)
                        ccGCAAGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGTGCTCGGCgagg                                                                                                                                     <  2:1304/131‑3 (MQ=255)
                           cAAGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTcc                                                                                                                <  5:203276/149‑1 (MQ=255)
                            aaGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCg                                                                                                               <  1:181384/149‑1 (MQ=255)
                                     gCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCaaaaa                                                                                                        <  8:92936/147‑1 (MQ=255)
                                     gCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCaaaaa                                                                                                        <  2:167493/147‑1 (MQ=255)
                                       cGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCaaaaa                                                                                                        <  7:108134/145‑1 (MQ=255)
                                       cGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCaaaaa                                                                                                        >  8:108134/1‑145 (MQ=255)
                                         aaCAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAg                                                                                                  <  1:131385/149‑1 (MQ=255)
                                                     gcCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGt                                                                                      <  1:50943/149‑1 (MQ=255)
                                                           gCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAag                                                                                <  2:277505/149‑1 (MQ=255)
                                                               tGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCaaaaa                                                                                                        <  7:219746/121‑1 (MQ=255)
                                                               tGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCaaaaa                                                                                                        >  8:219746/1‑121 (MQ=255)
                                                                    aCCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCagccagcc                                                                       >  6:115387/1‑149 (MQ=255)
                                                                    aCCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCAAGGTGCGACTGCGTGCCAAGcagccagc                                                                        >  8:157891/1‑148 (MQ=255)
                                                                      cAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAagcagc                                                                            >  7:300418/1‑142 (MQ=255)
                                                                      cAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAagcagc                                                                            <  8:300418/142‑1 (MQ=255)
                                                                       aGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCttt                                                                    <  6:230002/149‑1 (MQ=255)
                                                                       aGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCttt                                                                    >  6:153798/1‑149 (MQ=255)
                                                                         gCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCaaaaa                                                                                                        >  3:201754/1‑111 (MQ=255)
                                                                         gCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCaaaaa                                                                                                        <  4:201754/111‑1 (MQ=255)
                                                                         gCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGc                                                                  >  8:18002/1‑149 (MQ=255)
                                                                             cGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAATCAGCCAGCCTTTGCCaaa                                                              <  6:31947/149‑1 (MQ=255)
                                                                              ggCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCaaa                                                              >  5:181818/1‑148 (MQ=255)
                                                                              ggCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAac                                                             <  5:115387/149‑1 (MQ=255)
                                                                              ggCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAac                                                             <  7:157891/149‑1 (MQ=255)
                                                                                 cGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCaaa                                                              <  1:310103/145‑1 (MQ=255)
                                                                                 cGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCaaa                                                              >  2:310103/1‑145 (MQ=255)
                                                                                  gATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACAccc                                                         >  4:226226/1‑149 (MQ=255)
                                                                                       gcagcagGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCTAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCAAGCCTTTGCCAAACACcccgc                                                      >  6:161216/1‑147 (MQ=255)
                                                                                         agcagGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGACGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGa                                                  >  4:285568/1‑149 (MQ=255)
                                                                                          gcagGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAg                                                 >  8:166228/1‑149 (MQ=255)
                                                                                          gcagGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAg                                                 <  7:18002/149‑1 (MQ=255)
                                                                                          gcagGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGACCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCGTTGCCAAACACCCCGCCCGAg                                                 >  4:100524/1‑149 (MQ=255)
                                                                                              gTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGa                                                                                            >  7:279613/1‑102 (MQ=255)
                                                                                              gTAGACAATCAGCAGCAGGCAGATGCCTACCAGGACGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGa                                                                                            <  8:279613/102‑1 (MQ=255)
                                                                                               tAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCATTTGCCAACACCCCCGCCGGACCCgtt                                            >  3:191056/1‑147 (MQ=255)
                                                                                                 gACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGcagcca                                                                          <  4:259188/117‑1 (MQ=255)
                                                                                                 gACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGcagcca                                                                          >  3:259188/1‑117 (MQ=255)
                                                                                                 gACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATcg                                          <  5:153798/149‑1 (MQ=255)
                                                                                                   cAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATcgcc                                        <  3:226226/149‑1 (MQ=255)
                                                                                                      tCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGGAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGc                                                <  8:179566/138‑1 (MQ=255)
                                                                                                      tCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGc                                                >  7:179566/1‑138 (MQ=255)
                                                                                                        agcagcagGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCt                                    <  7:166228/148‑1 (MQ=255)
                                                                                                        agcagcagGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGAAGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCtt                                   >  6:154123/1‑149 (MQ=255)
                                                                                                         gcagcagGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTg                                  >  3:244810/1‑149 (MQ=255)
                                                                                                            gcaggcagATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGcc                                               <  4:145085/133‑1 (MQ=255)
                                                                                                            gcaggcagATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGcc                                               >  3:145085/1‑133 (MQ=255)
                                                                                                             caggcagATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACt                              >  3:190895/1‑149 (MQ=255)
                                                                                                               ggcagATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTg                             >  1:200733/1‑148 (MQ=255)
                                                                                                                gcagATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGga                           >  3:184942/1‑149 (MQ=255)
                                                                                                                  agATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGgat                          >  7:164918/1‑148 (MQ=255)
                                                                                                                   gATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGga                           <  6:112344/146‑1 (MQ=255)
                                                                                                                   gATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGga                           >  5:112344/1‑146 (MQ=255)
                                                                                                                     tGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGatgatg                      <  2:200733/149‑1 (MQ=255)
                                                                                                                        cTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGa                                >  3:130486/1‑136 (MQ=255)
                                                                                                                        cTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGa                                <  4:130486/136‑1 (MQ=255)
                                                                                                                        cTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTc                   <  8:164918/149‑1 (MQ=255)
                                                                                                                         tACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGcagcc                                                                           >  2:85640/1‑92 (MQ=255)
                                                                                                                         tACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGcagcc                                                                           <  1:85640/92‑1 (MQ=255)
                                                                                                                         tACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTc                   >  4:152306/1‑148 (MQ=255)
                                                                                                                           ccAGGCCGAATTCCTCGCCGAGCACGGCGAGGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCACGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCaa                >  1:107289/1‑149 (MQ=255)
                                                                                                                             aGGCCGAATTCCTCGCCGAGCACGGCGAAGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAAcc              <  5:154123/149‑1 (MQ=255)
                                                                                                                                ccGAATTCCTCGCCGAGCACGGCGATGATGAGGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGt           >  1:35142/1‑149 (MQ=255)
                                                                                                                                 cGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGTg          <  4:244810/149‑1 (MQ=255)
                                                                                                                                      tcatCGCCGTGCTCGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGTg          <  6:183861/141‑1 (MQ=255)
                                                                                                                                      tcatCGCCGTGCTCGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGTg          >  5:183861/4‑144 (MQ=255)
                                                                                                                                        cTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGTGCCCAACg   >  1:264843/1‑149 (MQ=255)
                                                                                                                                        cTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGTGCCCAACg   <  4:190895/149‑1 (MQ=255)
                                                                                                                                        cTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGTGCCCAACg   >  2:77327/1‑149 (MQ=255)
                                                                                                                                        cTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGTGCCCAACg   >  2:45183/1‑149 (MQ=255)
                                                                                                                                         tCGCAGAGCACGTCGATGATGAAGTCGGTGTGGCTTTCCGGAAAAAAGTCGAGGTGCGTCTGGGTGCCAAGCAGGCAGCCTTTGACAAACACCCCGCCGGAGCCGATCGCGGCCTTGGACTGgat                          <  7:255563/125‑1 (MQ=255)
                                                                                                                                            gCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGtt                    >  7:39637/1‑129 (MQ=255)
                                                                                                                                            gCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGtt                    <  8:39637/129‑1 (MQ=255)
                                                                                                                                            gCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGTGCCCAACgg  <  4:184942/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CAAAAAAAGTCATGGTCAGGCTGCCCGCAAGCAGCTTGCCGAACAGGGTCTGCGCCTGGGCGGTGATCACCAGGCCACGGCCGATCAGCAGCAGGTAGACAATCAGCAGCAGGCAGATGCCTACCAGGCCGAATTCCTCGCCGAGCACGGCGATGATGAAGTCGGTGTGGCTTTCCGGCAAAAAGTCCAGGTGCGACTGCGTGCCAAGCAGCCAGCCTTTGCCAAACACCCCGCCCGAGCCGATCGCCGCCTTGGACTGGATGATGTTCCAACCGGTGCCCAACGG  >  NC_002947/5469037‑5469322

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: