Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,720,210 G→A 11.3% L86L (CTG→CTA) ‡ pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,720,2100GA11.3% 208.2 / 15.7 71L86L (CTG→CTA) ‡pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base A (4/4);  total (35/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.68e-01

TTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGG  >  NC_000913/2720071‑2720351
                                                                                                                                           |                                                                                                                                             
ttAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGca                                                                                                                                       >  2:676530/1‑147 (MQ=255)
 tAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                      <  1:676530/148‑3 (MQ=255)
     ggACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCg                                                                                                                                  <  4:958531/148‑1 (MQ=255)
          ggTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGaaa                                                                                                                             >  3:965698/1‑148 (MQ=255)
                    gTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTaaaa                                                                                                                   <  4:965698/148‑1 (MQ=255)
                                  gAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCg                                                                                                                                  <  3:75311/119‑1 (MQ=255)
                                  gAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCg                                                                                                                                  >  4:75311/1‑119 (MQ=255)
                                      gCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                      >  3:5200/1‑109 (MQ=255)
                                      gCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                      <  4:5200/111‑3 (MQ=255)
                                           gtTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGca                                                                                                                                       >  1:1075693/1‑104 (MQ=255)
                                           gtTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGca                                                                                                                                       >  3:496848/1‑104 (MQ=255)
                                           gtTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGca                                                                                                                                       <  4:496848/105‑2 (MQ=255)
                                           gtTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGca                                                                                                                                       <  2:1075693/105‑2 (MQ=255)
                                                gtgtgtTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCt                                                                                                                >  3:653371/1‑123 (MQ=255)
                                                gtgtgtTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCt                                                                                                                <  4:653371/123‑1 (MQ=255)
                                                     ttGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAg                                                                                                                            <  1:144840/106‑1 (MQ=255)
                                                     ttGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAg                                                                                                                            >  2:144840/1‑106 (MQ=255)
                                                           tggccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGa                                                                                                                                             >  4:62448/1‑83 (MQ=255)
                                                           tggccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGa                                                                                                                                             <  3:62448/83‑1 (MQ=255)
                                                              ccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTaaaa                                                                                                                   >  4:1076039/1‑106 (MQ=255)
                                                              ccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTaaaa                                                                                                                   <  3:1076039/106‑1 (MQ=255)
                                                                ggATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCat                                                                                                            <  2:710179/111‑1 (MQ=255)
                                                                ggATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCat                                                                                                            >  1:710179/1‑111 (MQ=255)
                                                                ggATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGAtctc                                                                       >  1:158377/1‑148 (MQ=255)
                                                                ggATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGAtctc                                                                       >  2:194354/1‑148 (MQ=255)
                                                                  atatTGCCAGATTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCcg                                                                     >  3:651671/1‑148 (MQ=255)
                                                                   tatTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTaaaa                                                                                                                   <  1:422795/101‑1 (MQ=255)
                                                                   tatTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTaaaa                                                                                                                   >  2:422795/1‑101 (MQ=255)
                                                                             ttGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCg                                                                                                 >  1:309295/1‑109 (MQ=255)
                                                                             ttGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCg                                                                                                 <  2:309295/109‑1 (MQ=255)
                                                                                 ccTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTg                                                                                                               >  3:317797/1‑91 (MQ=255)
                                                                                 ccTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTg                                                                                                               <  4:317797/91‑1 (MQ=255)
                                                                                  cTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgcgccctgcgc                                                     >  4:233478/1‑148 (MQ=255)
                                                                                   tttACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                    >  2:195008/1‑100 (MQ=255)
                                                                                   tttACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                    <  1:195008/100‑1 (MQ=255)
                                                                                     tACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTTTTCTGGAAGAGCTCGGCGAGAAAGGCTg                                                                                                                        <  2:870305/78‑1 (MQ=255)
                                                                                     tACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTg                                                                                                                        >  1:870305/1‑78 (MQ=255)
                                                                                       cacCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGt                                                                                                                       >  2:572051/1‑77 (MQ=255)
                                                                                       cacCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGt                                                                                                                       <  1:572051/77‑1 (MQ=255)
                                                                                          ccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACAt                                             >  2:22307/1‑148 (MQ=255)
                                                                                           cGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgcgccctgcg                                                      >  2:71607/1‑138 (MQ=255)
                                                                                           cGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgcgccctgcg                                                      <  1:71607/138‑1 (MQ=255)
                                                                                           cGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATg                                            <  4:651671/148‑1 (MQ=255)
                                                                                             aCCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATgcg                                          <  2:158377/148‑1 (MQ=255)
                                                                                                 tGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgcgccctg                                                        >  2:555715/1‑130 (MQ=255)
                                                                                                 tGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgcgccctg                                                        <  1:555715/130‑1 (MQ=255)
                                                                                                 tGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTg                                      >  4:677662/1‑148 (MQ=255)
                                                                                                       tttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                    <  4:332842/80‑1 (MQ=255)
                                                                                                       tttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                    >  3:332842/1‑80 (MQ=255)
                                                                                                           tGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCaa                            <  3:233478/148‑1 (MQ=255)
                                                                                                           tGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCaa                            <  3:677662/148‑1 (MQ=255)
                                                                                                              aCCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACa                         >  1:794205/1‑148 (MQ=255)
                                                                                                               ccAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGa                                                                           <  3:791783/97‑1 (MQ=255)
                                                                                                               ccAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGa                                                                           >  4:791783/1‑97 (MQ=255)
                                                                                                                cAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGt                       <  2:794205/148‑1 (MQ=255)
                                                                                                                 aaTGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGa                                                                           <  2:300839/95‑1 (MQ=255)
                                                                                                                 aaTGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGa                                                                           >  1:300839/1‑95 (MQ=255)
                                                                                                                   tGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGa                                                                              >  3:902308/1‑90 (MQ=255)
                                                                                                                   tGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGa                                                                              <  4:902308/90‑1 (MQ=255)
                                                                                                                   tGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTg                                  <  2:496516/134‑1 (MQ=255)
                                                                                                                   tGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTg                                  >  1:496516/1‑134 (MQ=255)
                                                                                                                   tGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTg                    <  1:22307/148‑1 (MQ=255)
                                                                                                                        gCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                    <  2:538163/63‑1 (MQ=255)
                                                                                                                        gCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                    >  1:538163/1‑63 (MQ=255)
                                                                                                                              aTCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTg                    <  2:361907/137‑1 (MQ=255)
                                                                                                                              aTCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTg                    >  1:361907/1‑137 (MQ=255)
                                                                                                                                  tGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTa              <  3:949604/139‑1 (MQ=255)
                                                                                                                                  tGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTa              >  4:949604/1‑139 (MQ=255)
                                                                                                                                     tcttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCGCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTgg  <  1:194354/148‑1 (MQ=255)
                                                                                                                                       ttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGa                                                                              >  2:1109985/1‑70 (MQ=255)
                                                                                                                                       ttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGa                                                                              <  1:1109985/70‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
TTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGG  >  NC_000913/2720071‑2720351

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