Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,719,448 C→T 5.1% C52Y (TGC→TAC)  anmK ← anhydro‑N‑acetylmuramic acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,719,4480CT 5.1% 254.7 / 2.1 79C52Y (TGC→TAC) anmKanhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base C (41/33);  new base T (2/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGCGGCC  >  NC_000913/1719305‑1719593
                                                                                                                                               |                                                                                                                                                 
tttGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAata                                                                                                                                              >  5:333285/1‑149 (MQ=255)
      cGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGt                                                                                                                                        <  2:178915/149‑1 (MQ=255)
       gTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCa                                                                                                                                          <  6:21817/146‑1 (MQ=255)
       gTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCa                                                                                                                                          >  5:21817/1‑146 (MQ=255)
        tGACAACCGATCGCAACTAAATCTCTTGCCAGCAGATTTTGTTCCTTAAGAAAGACATTGACAGCATCAGCAAAAAGATGACCGAGTTGAGTATCAAGCTTACCAAACTCCGAAAGTGTAAGCAGCTGGCCCAGGCAAATATCCAGTa                                                                                                                                       >  4:94606/1‑148 (MQ=255)
         gACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCCGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTa                                                                                                                                       >  1:311445/1‑147 (MQ=255)
            aaCCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGcc                                                                                                                                  <  5:217422/149‑1 (MQ=255)
            aaCCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGAAAAAAGTTGCCAGAGTTGAGTATCAAGCTGTCAAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGAACAGcc                                                                                                                                  >  5:221784/1‑149 (MQ=255)
              ccGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATcc                                                                                                                                           <  8:109092/138‑1 (MQ=255)
              ccGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTg                                                                                                                                >  1:197133/1‑149 (MQ=255)
                     cAACTATATCTCTTGCCTTCAGATTTTGATCCTTAAGCAAGGCATTGACCGCATCAGAAAAAAGTTGCCCGAGTTCAGTATCAAGCTGTCCAAACCGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTCTTTCaga                                                                                                                         <  6:34443/149‑1 (MQ=255)
                     cAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTTAGCTGCCGGCCCTGGCAAATATCCAGTACAGCCTGTTTCaga                                                                                                                         <  7:70862/149‑1 (MQ=255)
                      aaCTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCagag                                                                                                                        <  8:126302/149‑1 (MQ=255)
                       aCTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAgaga                                                                                                                       <  6:206944/149‑1 (MQ=255)
                           tatCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCGTCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATAc                                                                                                                    <  5:85665/148‑1 (MQ=255)
                              ctctTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCgg                                                                                                                 >  2:79107/1‑148 (MQ=255)
                               tctTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGc                                                                                                                                   <  1:133547/129‑1 (MQ=255)
                               tctTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGc                                                                                                                                   >  2:133547/1‑129 (MQ=255)
                                       gCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGTAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCaa                                                                                                       >  6:294134/1‑149 (MQ=255)
                                       gCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGTAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCa                                                                                                        >  6:299453/1‑148 (MQ=255)
                                           aTTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTc                                                                                                    >  5:54093/1‑148 (MQ=255)
                                            ttttGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCTGCCAACTCa                                                                                                   >  5:155820/1‑148 (MQ=255)
                                             tttGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCaaa                                                                                                 >  7:245208/1‑149 (MQ=255)
                                              ttGTTCCTTAATCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTTAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTTGCAAATATCCAGTACAGCCTGTTTCAGTGATACCGTGATCGGCCAACTCaaa                                                                                                 <  2:319176/148‑1 (MQ=255)
                                              ttGTTCCTTAAGCAAGGCATTGACCGCATCAGCTAAAAGTTGCCCGATTTTAGTATCAAGCTGTCCAAACTGCGAATGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACATCCTTTTTCAGAGATACCGGGATCGGCCAACTCaaa                                                                                                 <  4:113631/148‑1 (MQ=255)
                                                 ttCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTg                                                                                             >  7:53824/1‑149 (MQ=255)
                                                     ttAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCaga                                                                                                                         >  5:297681/1‑117 (MQ=255)
                                                     ttAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCaga                                                                                                                         <  6:297681/117‑1 (MQ=255)
                                                     ttAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAg                                                                                         >  2:256120/1‑149 (MQ=255)
                                                        aGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCt                                                                                       <  6:54093/148‑1 (MQ=255)
                                                         gCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGcc                                                                                                         <  4:212950/129‑1 (MQ=255)
                                                         gCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGcc                                                                                                         >  3:212950/1‑129 (MQ=255)
                                                         gCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCtgt                                                                                     <  1:79107/149‑1 (MQ=255)
                                                                    aCCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCtgt                                                                                     >  6:293731/1‑138 (MQ=255)
                                                                    aCCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCtgt                                                                                     <  5:293731/138‑1 (MQ=255)
                                                                    aCCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgtg                                                                          >  3:12635/1‑149 (MQ=255)
                                                                    aCCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCCGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgtg                                                                          <  2:311445/149‑1 (MQ=255)
                                                                       gCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCACACTTGCCCGCTGTGCGACCCGGTGTtc                                                                       >  6:96250/1‑149 (MQ=255)
                                                                            cgCAAAAAGTTGCCGGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCGGTTTCAGTGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCTTCaa                                                                  <  7:359834/148‑1 (MQ=255)
                                                                            aGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCaa                                                                  <  8:245208/149‑1 (MQ=255)
                                                                              cAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAAtt                                                                >  1:152191/1‑149 (MQ=255)
                                                                              cAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAAtt                                                                >  3:154300/1‑149 (MQ=255)
                                                                               aaaaaGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTg                                                               >  4:295436/1‑149 (MQ=255)
                                                                                       gCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGcaaca                                                       >  6:232044/1‑145 (MQ=255)
                                                                                        cccGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaaca                                                      <  2:197133/149‑1 (MQ=255)
                                                                                         ccGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAGCTTTCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaaca                                                      >  1:382195/1‑148 (MQ=255)
                                                                                         ccGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGAccc                                                                              >  8:282947/1‑124 (MQ=255)
                                                                                         ccGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGAccc                                                                              <  7:282947/124‑1 (MQ=255)
                                                                                         ccGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGATGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACAGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaaca                                                      >  2:302770/1‑148 (MQ=255)
                                                                                           gAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacacaa                                                   >  1:164120/1‑149 (MQ=255)
                                                                                            aGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacacaac                                                  <  6:333285/149‑1 (MQ=255)
                                                                                            aGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacacaa                                                   >  5:8585/1‑148 (MQ=255)
                                                                                                gAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGTAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATc                                               <  5:299453/148‑1 (MQ=255)
                                                                                                gAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGTAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATc                                               <  5:294134/148‑1 (MQ=255)
                                                                                                gAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATc                                               >  2:45535/1‑148 (MQ=255)
                                                                                                gAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGATGCTGGAGATTGCAAATATCAAGTACAGCCTGTTTCTTAGTTTCCGGGTTCGGCCAACTCTAACTTGCCATCTGTGCGAGCGGGTGTTCATCATTTGTCGCCAACACAAGATc                                               <  6:221784/148‑1 (MQ=255)
                                                                                                    aTCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACGGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGCCCCGGTGTTCATCAATTGTCGCCAACACAACATCAacac                                          >  2:198228/1‑149 (MQ=255)
                                                                                                          cTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCAAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCaa                                    <  3:94606/149‑1 (MQ=255)
                                                                                                                aaaCTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTgg                              >  1:266220/1‑149 (MQ=255)
                                                                                                                aaaCTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTgg                              >  7:385839/1‑149 (MQ=255)
                                                                                                                   cTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGc                           <  4:154300/149‑1 (MQ=255)
                                                                                                                      cGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTg                        <  3:295436/149‑1 (MQ=255)
                                                                                                                       gAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGCCCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTgc                       >  5:33735/1‑148 (MQ=255)
                                                                                                                       gAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGc                                                           <  7:257668/113‑1 (MQ=255)
                                                                                                                       gAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGc                                                           >  8:257668/1‑113 (MQ=255)
                                                                                                                           gtgtAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACAt                    >  7:316458/1‑148 (MQ=255)
                                                                                                                              tAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACg                >  6:157736/1‑149 (MQ=255)
                                                                                                                                aggtgctgGCCCTGGAAAATATCCAGTACAGCCTGTTTCAGAGATACCTGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGc               <  8:385839/145‑1 (MQ=255)
                                                                                                                                aGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGc               >  2:211474/1‑148 (MQ=255)
                                                                                                                                aGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGc               <  1:211474/148‑1 (MQ=255)
                                                                                                                                       ggCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGc                                                                                            >  4:9040/1‑64 (MQ=255)
                                                                                                                                       ggCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGc                                                                                            <  3:9040/64‑1 (MQ=255)
                                                                                                                                        gCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGc      >  5:165569/1‑149 (MQ=255)
                                                                                                                                          ccTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGa                       >  5:202170/1‑130 (MQ=255)
                                                                                                                                          ccTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGa                       <  6:202170/130‑1 (MQ=255)
                                                                                                                                          ccTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGc               >  5:51207/1‑138 (MQ=255)
                                                                                                                                          ccTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGc               <  6:51207/138‑1 (MQ=255)
                                                                                                                                          ccTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGCgg    >  5:41103/1‑149 (MQ=255)
                                                                                                                                              gCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGCGGcc  <  1:256120/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGCGGCC  >  NC_000913/1719305‑1719593

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: