Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 702,544 C→A 9.5% D70Y (GAC→TAC)  nagA ← N‑acetylglucosamine‑6‑phosphate deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913702,5440CA 9.5% 235.5 / 9.9 72D70Y (GAC→TAC) nagAN‑acetylglucosamine‑6‑phosphate deacetylase
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base A (4/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CAGGTACTCGCGCATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAA  >  NC_000913/702398‑702690
                                                                                                                                                  |                                                                                                                                                  
cAGGTACTCGCGCATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCg                                                                                                                                                   <  5:98458/148‑1 (MQ=255)
 aGGTACTCGCGCATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTAGt                                                                                                                                                  >  8:116147/1‑148 (MQ=255)
    tACTGGCGCATAACGCGCAAGCCCGGTTTAATCTGCTCATCGCTGGTGGTGATAAGCGTCTTCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTAGTTaa                                                                                                                                               <  1:268993/148‑1 (MQ=255)
       tCGCGCATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTAAAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTAAGAGGTGTCGTTAAACt                                                                                                                                            >  2:160050/1‑148 (MQ=255)
         gcgcATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGt                                                                                                                                          <  5:277390/148‑1 (MQ=255)
         gcgcATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGt                                                                                                                                          <  6:118189/148‑1 (MQ=255)
          cgcATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGtt                                                                                                                                                 <  8:18990/140‑1 (MQ=255)
          cgcATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGtt                                                                                                                                                 >  7:18990/1‑140 (MQ=255)
          cgcATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTa                                                                                                                                         <  4:199187/148‑1 (MQ=255)
            cATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTAcg                                                                                                                                       >  7:321348/1‑148 (MQ=255)
              tAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTAcgcc                                                                                                                                     >  2:94762/1‑148 (MQ=255)
                         cgcTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGt                                                                                                                          <  2:207392/146‑1 (MQ=255)
                               ttCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTAGTTAAACTGTACGCCGCCGCAGCCGTTTAACt                                                                                                                    >  7:196326/1‑148 (MQ=255)
                                 catcaGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTg                                                                                                                   <  2:292050/147‑1 (MQ=255)
                                   tcaGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGcc                                                                                                                            >  4:4733/1‑136 (MQ=255)
                                   tcaGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGcc                                                                                                                            <  3:4733/136‑1 (MQ=255)
                                    caGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTATACTGTACGCCGCCGCAGCCGTTTAACTGcaca                                                                                                               >  6:40701/1‑148 (MQ=255)
                                            cGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATaa                                                                                                        >  2:96260/1‑147 (MQ=255)
                                              cTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTaaa                                                                                                                                              >  2:308825/1‑107 (MQ=255)
                                              cTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTaaa                                                                                                                                              <  1:308825/107‑1 (MQ=255)
                                                  tggtgATAAGCGTCGTCAGATAGTTATTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCATCCGTTTAACTGCACATCGATAAAACCggg                                                                                                 <  5:252259/148‑1 (MQ=255)
                                                   ggtgATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCgggg                                                                                                >  3:272094/1‑148 (MQ=255)
                                                     tgATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCgggg                                                                                                >  1:220174/1‑146 (MQ=255)
                                                     tgATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCgggg                                                                                                >  8:94534/1‑146 (MQ=255)
                                                       aTAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTAGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGgaga                                                                                            <  8:196326/148‑1 (MQ=255)
                                                        tAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGga                                                                                              >  2:63669/1‑145 (MQ=255)
                                                         aaGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGcaca                                                                                                               >  2:311583/1‑127 (MQ=255)
                                                         aaGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGcaca                                                                                                               <  1:311583/127‑1 (MQ=255)
                                                          aGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCGTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGgagaga                                                                                          <  2:164864/147‑1 (MQ=255)
                                                              tCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGc                                                                                     <  3:71591/148‑1 (MQ=255)
                                                                ggCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGccc                                                                                   <  8:321348/148‑1 (MQ=255)
                                                                  cAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCg                                                                                 >  8:248851/1‑148 (MQ=255)
                                                                     aTAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTc                                                                              <  1:94762/148‑1 (MQ=255)
                                                                       agttagtACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTAGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAg                                                                            >  2:237634/1‑148 (MQ=255)
                                                                        gttagtACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGt                                                                           >  5:196656/1‑148 (MQ=255)
                                                                              aCAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACgt                                                                     >  5:263272/1‑148 (MQ=255)
                                                                               cAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACgtt                                                                    >  2:314242/1‑148 (MQ=255)
                                                                                aGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACgttg                                                                   <  1:160050/148‑1 (MQ=255)
                                                                                 gCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACgttgt                                                                  <  2:220174/148‑1 (MQ=255)
                                                                                  ccTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGgagaga                                                                                          <  2:220428/123‑1 (MQ=255)
                                                                                  ccTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGgagaga                                                                                          >  1:220428/1‑123 (MQ=255)
                                                                                  ccTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACgttgtt                                                                 <  7:248851/148‑1 (MQ=255)
                                                                                        ttCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGAtct                                                           >  2:179686/1‑148 (MQ=255)
                                                                                           tcATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCggggg                                                                                               <  7:293367/109‑1 (MQ=255)
                                                                                           tcATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCggggg                                                                                               >  8:293367/1‑109 (MQ=255)
                                                                                              ttGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTggcg                                                     >  6:136714/1‑148 (MQ=255)
                                                                                               tGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGaa                                                                        >  7:323261/1‑128 (MQ=255)
                                                                                               tGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGaa                                                                        <  8:323261/128‑1 (MQ=255)
                                                                                                 gCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCTGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCa                                                  <  1:96260/148‑1 (MQ=255)
                                                                                                 gCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTggcggc                                                   <  5:40701/147‑1 (MQ=255)
                                                                                                   tttCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTAGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGt                                                >  8:25788/1‑148 (MQ=255)
                                                                                                    ttCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGgagaga                                                                                          >  7:112689/1‑105 (MQ=255)
                                                                                                    ttCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGgagaga                                                                                          <  8:112689/105‑1 (MQ=255)
                                                                                                         cATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTAGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCg                                            <  1:237634/146‑1 (MQ=255)
                                                                                                           tGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTAAAACTGTACGCTGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTAc                                        <  8:304484/148‑1 (MQ=255)
                                                                                                             aTTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTAcc                                       <  6:263272/147‑1 (MQ=255)
                                                                                                              tttCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCgg                                     <  3:146007/148‑1 (MQ=255)
                                                                                                                tCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGa                                    <  1:59768/147‑1 (MQ=255)
                                                                                                                     cGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACg                              <  1:63669/148‑1 (MQ=255)
                                                                                                                       tttCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCt                            >  5:77068/1‑148 (MQ=255)
                                                                                                                         tCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCttt                          >  3:205581/1‑148 (MQ=255)
                                                                                                                           cACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCtttt                         >  1:78164/1‑147 (MQ=255)
                                                                                                                           cACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCtttt                         >  4:51241/1‑147 (MQ=255)
                                                                                                                            aCGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTaa                       <  6:196656/148‑1 (MQ=255)
                                                                                                                              gCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATc                     <  1:179686/148‑1 (MQ=255)
                                                                                                                              gCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATc                     <  5:136714/148‑1 (MQ=255)
                                                                                                                               cTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTg                                                        >  6:104056/1‑112 (MQ=255)
                                                                                                                               cTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTg                                                        <  5:104056/112‑1 (MQ=255)
                                                                                                                                  aCCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAg                   <  7:33875/146‑1 (MQ=255)
                                                                                                                                  aCCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACGGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGc                 >  8:244651/1‑148 (MQ=255)
                                                                                                                                                gtCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTa                        <  4:23641/127‑1 (MQ=255)
                                                                                                                                                gtCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTa                        >  3:23641/1‑127 (MQ=255)
                                                                                                                                                  cGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATaacaa  >  6:120657/1‑147 (MQ=255)
                                                                                                                                                  cGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATaacaa  <  5:120657/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGGTACTCGCGCATAACGCGCACGCCCTGTTTCATCAGCTCATCGCTGGTGGTGATAAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAA  >  NC_000913/702398‑702690

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