Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 491,948 G→A 5.9% P179P (CCG→CCA)  apt → adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913491,9480GA 5.9% 217.8 / 2.3 68P179P (CCG→CCA) aptadenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCACCTTCCAGCGTTTCAGAGCCTGCCAAT  >  NC_000913/491808‑492093
                                                                                                                                            |                                                                                                                                                 
ggcACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCggg                                                                                                                                           <  7:34151/149‑1 (MQ=255)
     tATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCat                                                                                                                                       >  8:208267/1‑148 (MQ=255)
       tCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACACTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCattaa                                                                                                                                    <  3:60681/149‑1 (MQ=255)
         gAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCattaat                                                                                                                                   >  7:100426/1‑148 (MQ=255)
            gCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATc                                                                                                                               <  2:216327/149‑1 (MQ=255)
            gCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATc                                                                                                                               <  2:59970/149‑1 (MQ=255)
             cGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAACCAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCg                                                                                                                              <  3:73678/149‑1 (MQ=255)
               aCCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGGCATTAATTATCGc                                                                                                                             >  6:125095/1‑148 (MQ=255)
                ccGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGc                                                                                                                             >  7:152868/1‑147 (MQ=255)
                ccGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCa                                                                                                                           <  5:33114/149‑1 (MQ=255)
                      aaCTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTtg                                                                                                                     >  3:58481/1‑149 (MQ=255)
                      aaCTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTtg                                                                                                                     >  5:170668/1‑149 (MQ=255)
                      aaCTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTtg                                                                                                                     >  5:176066/1‑149 (MQ=255)
                       aCTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTtgt                                                                                                                    >  8:61828/1‑149 (MQ=255)
                       aCTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTtgt                                                                                                                    >  5:239098/1‑149 (MQ=255)
                       aCTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTtgt                                                                                                                    >  5:75214/1‑149 (MQ=255)
                       aCTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTtgt                                                                                                                    >  8:169410/1‑149 (MQ=255)
                           aTCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGtgct                                                                                                                 <  3:29493/148‑1 (MQ=255)
                              cgtcgtCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGcc                                                                                                              >  8:217992/1‑148 (MQ=255)
                              cgtcgtCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCGCGGCGGCGAACAGCGGCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGGCTTGTGCTGcc                                                                                                              <  7:204605/148‑1 (MQ=255)
                               gtcgtcTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCa                                                                                                            <  3:1829/149‑1 (MQ=255)
                               gtcgtcTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCa                                                                                                            >  6:14686/1‑149 (MQ=255)
                               gtcgtcTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCa                                                                                                            <  4:170797/149‑1 (MQ=255)
                                   tcTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGc                                                                                                         >  1:226268/1‑148 (MQ=255)
                                       ggtggtGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGTCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTAc                                                                                                      <  5:125095/147‑1 (MQ=255)
                                       ggtggtGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTAc                                                                                                      <  7:217992/147‑1 (MQ=255)
                                       ggtggtGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACgg                                                                                                    >  2:154871/1‑149 (MQ=255)
                                        gtggtgAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGa                                                                                                   <  7:208267/149‑1 (MQ=255)
                                         tggtgAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTtgtg                                                                                                                   <  4:45841/132‑1 (MQ=255)
                                         tggtgAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTtgtg                                                                                                                   >  3:45841/1‑132 (MQ=255)
                                             gAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGca                                                                                              <  8:152868/149‑1 (MQ=255)
                                                   gCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGAt                                                                                        >  7:90818/1‑149 (MQ=255)
                                                   gCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGAt                                                                                        <  6:75214/149‑1 (MQ=255)
                                                       aCGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCTCACGCTACGGACAGCACAAGAtgtg                                                                                     >  5:238887/1‑148 (MQ=255)
                                                       aCGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGc                                                                                    <  4:58481/149‑1 (MQ=255)
                                                        cGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCa                                                                                   <  5:82683/149‑1 (MQ=255)
                                                              gTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGc                                                                             <  3:186420/149‑1 (MQ=255)
                                                                     aTCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCggg                                                                                                                                           >  8:215804/1‑80 (MQ=255)
                                                                     aTCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCggg                                                                                                                                           <  7:215804/80‑1 (MQ=255)
                                                                        aaCCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGtt                                                                    <  1:39859/148‑1 (MQ=255)
                                                                         aCCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCt                                                                  >  2:96125/1‑149 (MQ=255)
                                                                           cTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGa                                                                <  2:226268/149‑1 (MQ=255)
                                                                           cTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGa                                                                <  6:176066/149‑1 (MQ=255)
                                                                             gTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACg                                                              >  7:7644/1‑149 (MQ=255)
                                                                                 gATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACgggg                                                           <  8:7644/148‑1 (MQ=255)
                                                                                  aTCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCt                                                         >  1:92451/1‑149 (MQ=255)
                                                                                        gcggcgAACAGCGTCTCGAAAAACAGGGCATTACCCGCTACAGCCTTGTCCCGTTCGCGGGCCATTATTTATCGCCAGTCTTGTGCTGCCCAGGCTACGGACTGCACTAGATGTGCATTCAGCCTCGCCGTTCTGTCGGGGCTGTGTTa                                                   <  3:177490/149‑1 (MQ=255)
                                                                                        gcggcgAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCTCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTa                                                   <  6:238887/149‑1 (MQ=255)
                                                                                        gcggcgAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTa                                                   <  7:61828/149‑1 (MQ=255)
                                                                                            cgAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCAt                                               <  2:92451/149‑1 (MQ=255)
                                                                                               aCAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTGGCATTAc                                            >  7:50911/1‑149 (MQ=255)
                                                                                                     tctcGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTc                                      <  1:96125/149‑1 (MQ=255)
                                                                                                       tcGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGt                                    >  3:100870/1‑149 (MQ=255)
                                                                                                         gAAAAACAGGGCATTACCAGCGACAGCCTTGTCCCGTTCCCGGGCCAGTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGGGTTAGCATTACCCCGTCGTg                                   >  6:156085/1‑148 (MQ=255)
                                                                                                             aaCAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATcc                              <  8:90818/149‑1 (MQ=255)
                                                                                                               cAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCAAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCAc                            >  5:9747/1‑149 (MQ=255)
                                                                                                                         ccAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCACCtt                         <  2:154417/142‑1 (MQ=255)
                                                                                                                         ccAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCACCtt                         >  1:154417/1‑142 (MQ=255)
                                                                                                                             cTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCACCTTCCAGCGTTTCa              <  7:169410/149‑1 (MQ=255)
                                                                                                                               aCAGCCTTGTCCCATTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACTGGGCTGTGTTAGCATTACCCCTTCGTGAATCCAc                            <  6:103527/133‑1 (MQ=255)
                                                                                                                               aCAGCCTTGTCCCATTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCAc                            <  4:115533/133‑1 (MQ=255)
                                                                                                                               aCAGCCTTGTCCCATTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCAc                            >  5:103527/1‑133 (MQ=255)
                                                                                                                               aCAGCCTTGTCCCATTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCAc                            >  3:115533/1‑133 (MQ=255)
                                                                                                                                cAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCACCTTCCAGCGTTTCagag           >  5:236730/1‑149 (MQ=255)
                                                                                                                                        tCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCACCTTCCAGCGTTTCAGAGCCTGCCaa   >  8:30073/1‑149 (MQ=255)
                                                                                                                                        tCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCACCTTCCAGCGTTTCAGAGCCTGCCaa   >  8:93457/1‑149 (MQ=255)
                                                                                                                                          ccGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCACCTTCCAGCGTTTCAGAGCCTGCCAAt  >  3:159803/1‑148 (MQ=255)
                                                                                                                                          ccGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCACCTTCCAGCGTTTCAGAGCCTGCCAAt  >  6:108554/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GGCACTATCGAAGCGACCGTTAAACTGATCCGTCGTCTGGGTGGTGAAGTGGCTGACGCTGCGTTCATTATCAACCTGTTCGATCTCGGCGGCGAACAGCGTCTCGAAAAACAGGGCATTACCAGCTACAGCCTTGTCCCGTTCCCGGGCCATTAATTATCGCCAGTCTTGTGCTGCCCACGCTACGGACAGCACAAGATGTGCATTCAGCCTCGCCGTTCTGACGGGGCTGTGTTAGCATTACCCCTTCGTGAATCCACCTTCCAGCGTTTCAGAGCCTGCCAAT  >  NC_000913/491808‑492093

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: ‑