Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 961,355 T→G 5.7% V52G (GTT→GGT)  cmk → cytidylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913961,3550TG 5.7% 212.9 / 6.3 70V52G (GTT→GGT) cmkcytidylate kinase
Reads supporting (aligned to +/- strand):  ref base T (35/31);  new base G (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

TTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGC  >  NC_000913/961211‑961499
                                                                                                                                                |                                                                                                                                                
ttGCCCCGGTTATTACCATTGATGGCCCACGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCatgttga                                                                                                                                               <  5:226978/148‑1 (MQ=255)
ttGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCatgttga                                                                                                                                               <  5:424705/148‑1 (MQ=255)
     ccGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCatgttgatgttg                                                                                                                                          <  5:125116/148‑1 (MQ=255)
       ggTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCg                                                                                                                                        >  7:228719/1‑148 (MQ=255)
             taCCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCatgttgatgt                                                                                                                                            >  5:345255/1‑138 (MQ=255)
             taCCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCatgttgatgt                                                                                                                                            <  6:345255/138‑1 (MQ=255)
                  ttGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATgc                                                                                                                             >  1:380618/1‑148 (MQ=255)
                     aTGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGATGCGTTGCAATGGCATCTGCTGGACTCGTGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCt                                                                                                                          <  6:372346/148‑1 (MQ=255)
                          ccAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGACGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTAc                                                                                                                     >  5:321285/1‑148 (MQ=255)
                          ccAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTAc                                                                                                                     >  8:259425/1‑148 (MQ=255)
                           cAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTAcc                                                                                                                    >  7:25360/1‑148 (MQ=255)
                           cAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGAGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTAcc                                                                                                                    >  1:79940/1‑148 (MQ=255)
                            aaGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCatgttgatg                                                                                                                                             >  4:358280/1‑122 (MQ=255)
                            aaGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCatgttgatg                                                                                                                                             <  3:358280/122‑1 (MQ=255)
                              gCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAg                                                                                                                                 >  5:210972/1‑132 (MQ=255)
                              gCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAg                                                                                                                                 <  6:210972/132‑1 (MQ=255)
                                       tgAAATGCACATTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCGTCTGTACTCGGGTGCAATTTATCGCGTAATGGCATTGGCTGCATTACATCACCATGTTGATGTTGCGTCGTAATATGCGCTGGTACAGCTGTCATccc                                                                                                         >  2:338896/2‑147 (MQ=255)
                                              cACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCatgttg                                                                                                                                                >  8:42502/1‑101 (MQ=255)
                                              cACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCatgttg                                                                                                                                                <  7:42502/101‑1 (MQ=255)
                                                ccTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATg                                                                                                >  3:197065/1‑147 (MQ=255)
                                                 cTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTa                                                                                              >  8:449532/1‑148 (MQ=255)
                                                 cTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTa                                                                                              <  5:94094/148‑1 (MQ=255)
                                                  ttGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGCACTCGGGTGCAATTTATAGCGTACTGGCATTCGCGGCATTACATAACAATGTTGATGTTGCGTCGGAATATGCGCTGGTACCGCTGGCATCCAAACTGGATGTAc                                                                                             >  8:360169/1‑148 (MQ=255)
                                                     tgtAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGtt                                                                                          >  4:189615/1‑148 (MQ=255)
                                                       tAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCATTTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCTCTGGTACCGCTTGCATCCCATCTGGATGTACGttt                                                                                         <  7:360169/147‑1 (MQ=255)
                                                       tAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGtttt                                                                                        <  7:259425/148‑1 (MQ=255)
                                                         aGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATTTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTtg                                                                                       <  1:338896/147‑1 (MQ=255)
                                                          ggCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTtgtg                                                                                     <  2:380618/148‑1 (MQ=255)
                                                              aTGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGa                                                                                 >  3:6576/1‑148 (MQ=255)
                                                                  cGGAAGCGTTGCAATGGCATGTGATGGACTCGTGTGCAATTTATCGCGTACTGGCATTTGCGTCATTAAATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCaa                                                                             <  8:25360/148‑1 (MQ=255)
                                                                  cGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCaa                                                                             <  4:197065/148‑1 (MQ=255)
                                                                           tGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCt                                                                    >  6:437269/1‑148 (MQ=255)
                                                                              aaTGGCTTCTGCTGGACTCGGGTGCAATGTATCGCGTAGTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCTGAAGATTCGCTGTTACCGCTGGCATCGCATCTGGATGTACTTTTTTTGTCGACCAATGGCAATCTGGa                                                                 <  6:321285/148‑1 (MQ=255)
                                                                              aaTGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGa                                                                 <  8:95241/148‑1 (MQ=255)
                                                                              aaTGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGa                                                                 <  3:152483/148‑1 (MQ=255)
                                                                               aTGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTgg                                                                  >  7:288523/1‑146 (MQ=255)
                                                                                tGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAg                                                               >  6:390239/1‑148 (MQ=255)
                                                                                    aTCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt                                                           >  8:137689/1‑148 (MQ=255)
                                                                                     tCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCa                                                                                                             <  7:197379/97‑1 (MQ=255)
                                                                                     tCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCa                                                                                                             >  8:197379/1‑97 (MQ=255)
                                                                                     tCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATc                                                          >  8:15138/1‑148 (MQ=255)
                                                                                       tgctgGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCt                                                        >  7:364351/1‑148 (MQ=255)
                                                                                         ctgGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTc                                                       >  8:3007/1‑147 (MQ=255)
                                                                                               tCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAgggg                                                <  8:364351/148‑1 (MQ=255)
                                                                                                  ggTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAg                                             <  5:437269/148‑1 (MQ=255)
                                                                                                       aaTTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTc                                        <  4:6576/148‑1 (MQ=255)
                                                                                                              cgcgTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGAt                                           >  7:293475/1‑138 (MQ=255)
                                                                                                              cgcgTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGAt                                           <  8:293475/138‑1 (MQ=255)
                                                                                                                 gTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTg                                                                   <  1:3533/111‑1 (MQ=255)
                                                                                                                 gTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTg                                                                   >  2:3533/1‑111 (MQ=255)
                                                                                                                    cTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAAtt                            <  8:288523/147‑1 (MQ=255)
                                                                                                                       gCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTa                        <  4:290713/148‑1 (MQ=255)
                                                                                                                         aTTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACt                      <  7:3007/148‑1 (MQ=255)
                                                                                                                           tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACt                      >  8:422121/1‑146 (MQ=255)
                                                                                                                           tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCGGGATGTACGTTTTGTGTCGACCAATGGAAATCTCGAAGGGATCCGAGAAGGGGAAGATGTCAGCGGACAAATACGTACt                      >  2:238077/1‑146 (MQ=255)
                                                                                                                            ggcggcATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCaa                                                                       <  4:141319/96‑1 (MQ=255)
                                                                                                                            ggcggcATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCaa                                                                       >  3:141319/1‑96 (MQ=255)
                                                                                                                            ggcggcATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCa                    <  8:228719/147‑1 (MQ=255)
                                                                                                                              cggcATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAgg                  >  3:18581/1‑147 (MQ=255)
                                                                                                                               ggcATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGaa                <  7:422121/148‑1 (MQ=255)
                                                                                                                                    tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTgtgt                                                                                    >  5:206917/1‑75 (MQ=255)
                                                                                                                                    tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTgtgt                                                                                    <  6:206917/75‑1 (MQ=255)
                                                                                                                                    tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAg                                             >  8:71964/1‑114 (MQ=255)
                                                                                                                                    tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAg                                             <  7:71964/114‑1 (MQ=255)
                                                                                                                                    tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGc           >  2:279350/1‑148 (MQ=255)
                                                                                                                                    tACATCACCATGTGGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGc           >  5:40042/1‑148 (MQ=255)
                                                                                                                                     aCATCAACATGGTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGc           >  3:10544/1‑147 (MQ=255)
                                                                                                                                     aCATCAACATGGTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGc           >  7:181678/1‑147 (MQ=255)
                                                                                                                                      cATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCaa                                                                       <  8:293765/86‑1 (MQ=255)
                                                                                                                                      cATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCaa                                                                       >  7:293765/1‑86 (MQ=255)
                                                                                                                                      cATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCg          >  1:273641/1‑147 (MQ=255)
                                                                                                                                      cATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGa         <  1:279350/148‑1 (MQ=255)
                                                                                                                                         cACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCa                                                                                                             <  4:107107/45‑1 (MQ=255)
                                                                                                                                         cACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCa                                                                                                             >  3:107107/1‑45 (MQ=255)
                                                                                                                                          aacATGGTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATg      <  8:181678/145‑1 (MQ=255)
                                                                                                                                          aacATGGTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATg      <  4:10544/145‑1 (MQ=255)
                                                                                                                                          aCCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATg      <  7:15138/147‑1 (MQ=255)
                                                                                                                                           ccATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAgg                  >  5:169226/1‑134 (MQ=255)
                                                                                                                                           ccATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAgg                  <  6:169226/134‑1 (MQ=255)
                                                                                                                                             atgttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGc  >  4:177391/1‑148 (MQ=255)
                                                                                                                                               gttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAg                                             <  3:257973/103‑1 (MQ=255)
                                                                                                                                               gttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAg                                             >  4:257973/1‑103 (MQ=255)
                                                                                                                                                ttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTg                                                                   >  6:260440/1‑80 (MQ=255)
                                                                                                                                                ttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTg                                                                   <  5:260440/80‑1 (MQ=255)
                                                                                                                                                ttgatgttgCGTCGGAAGATGCGCCGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGt                                         <  2:195127/106‑1 (MQ=255)
                                                                                                                                                ttgatgttgCGTCGGAAGATGCGCCGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGt                                         >  1:195127/1‑106 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTGCCCCGGTTATTACCATTGATGGCCCAAGCGGTGCAGGGAAAGGCACCTTGTGTAAGGCTATGGCGGAAGCGTTGCAATGGCATCTGCTGGACTCGGGTGCAATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGC  >  NC_000913/961211‑961499

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: