Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,906 A→T 5.2% M201K (ATG→AAG)  glyA ← serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,684,9060AT 5.2% 254.2 / 3.3 77M201K (ATG→AAG) glyAserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base T (2/2);  total (38/39)
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.15e-01

TACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGAT  >  NC_000913/2684762‑2685051
                                                                                                                                                |                                                                                                                                                 
tACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCata                                                                                                                                                 <  8:368774/147‑1 (MQ=255)
  cAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCa                                                                                                                                              <  4:6961/148‑1 (MQ=255)
   aGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCaa                                                                                                                                             <  6:340455/148‑1 (MQ=255)
           gCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAgg                                                                                                                                     >  7:71492/1‑148 (MQ=255)
                     ttCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGa                                                                                                                           >  4:28456/1‑148 (MQ=255)
                     ccgttctggATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCg                                                                                                                            <  5:202524/140‑1 (MQ=255)
                      tCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCg                                                                                                                            <  6:294728/146‑1 (MQ=255)
                      tCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGGTTTGTGTGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCg                                                                                                                            <  2:268506/146‑1 (MQ=255)
                        gCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCTGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATg                                                                                                                         <  1:74662/147‑1 (MQ=255)
                          taGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAACCACCGATGCtc                                                                                                                      >  8:336858/2‑147 (MQ=255)
                          cAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTg                                                                                                                      >  2:163021/1‑148 (MQ=255)
                           aGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGt                                                                                                                     >  8:139436/1‑148 (MQ=255)
                             gATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGt                                                                                                                     >  1:194142/1‑146 (MQ=255)
                               tCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGTGTAGACGCAAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGtgcgg                                                                                                                 >  6:188295/1‑144 (MQ=255)
                                 aGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCAACGATGCTGTCAGcct                                                                                                               >  2:215862/1‑146 (MQ=255)
                                  ggCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTCGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGCTGCTGTCAGCGAt                                                                                                              >  5:91166/1‑148 (MQ=255)
                                   gccgccGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCata                                                                                                                                                 >  1:203598/1‑112 (MQ=255)
                                   gccgccGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCata                                                                                                                                                 <  2:203598/112‑1 (MQ=255)
                                   gccgccGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGAtt                                                                                                             <  3:28456/148‑1 (MQ=255)
                                    ccgccgCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGAttt                                                                                                            >  7:191117/1‑148 (MQ=255)
                                       ccgcGCGGACCCGCCAGGGTTTTGTGAGTGTTGGTAGTAACAACGTGAGCATGAGGAACCGGGGTCGGGTAGACGCCAGCAGCAACCATGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCAc                                                                                                         <  2:64917/148‑1 (MQ=255)
                                        cgcgcgGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACg                                                                                                        <  1:163021/148‑1 (MQ=255)
                                              gACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCAtttt                                                                                                  >  1:326285/1‑148 (MQ=255)
                                                   gCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGccc                                                                                             <  3:46184/148‑1 (MQ=255)
                                                     cAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAg                                                                                           >  4:224432/1‑148 (MQ=255)
                                                          ttttGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccac                                                                                      <  7:139436/148‑1 (MQ=255)
                                                           tttGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGCACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccacc                                                                                     >  1:41396/1‑148 (MQ=255)
                                                            ttGTGAGTGGTGGTATTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGACAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAATCACCGATTCTGTAAGCTATTTCACGCATTTTAGCCCAGTccacc                                                                                     <  7:336858/147‑1 (MQ=255)
                                                            ttGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccacc                                                                                     <  3:224432/147‑1 (MQ=255)
                                                                gAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGccc                                                                                             >  8:38290/1‑135 (MQ=255)
                                                                gAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGccc                                                                                             <  7:38290/135‑1 (MQ=255)
                                                                 aGTGGTGTAGTAACAACGTGAGCATGAGGAACCGGGTGCGGGTAGACGCCAGCGGCAACCAGGCCCGCAACGTGCGCCAGTTCAACGATCAGGTAAGCAACGAGGCTGTCAGCGATTTCACGCTTTTTCGCCCAGTCCACCccgccg                                                                               >  2:258299/1‑147 (MQ=255)
                                                                   tggtggTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGa                                                                             <  2:261806/148‑1 (MQ=255)
                                                                     gtggtAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACATGTAAGCACAGATGATGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAat                                                                           <  5:124771/148‑1 (MQ=255)
                                                                        gtagtaACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAatat                                                                         <  2:194142/147‑1 (MQ=255)
                                                                        gtagtaACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAatat                                                                         <  1:215862/147‑1 (MQ=255)
                                                                           gtaACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCa                                                                      <  4:294221/147‑1 (MQ=255)
                                                                            taACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCaga                                                                    <  7:88062/148‑1 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCTTATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAat                                                                           >  4:175192/1‑140 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCTTATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAat                                                                           >  4:169434/1‑140 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCTTATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAat                                                                           <  3:169434/140‑1 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCTTATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAat                                                                           <  3:175192/140‑1 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCAcc                                                                                                                             >  7:128753/1‑90 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCAcc                                                                                                                             <  8:128753/90‑1 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCCGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCa                                                                      >  7:173753/1‑145 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccacc                                                                                     >  7:51478/1‑130 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTccacc                                                                                     <  8:51478/130‑1 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCa                                                                      >  5:50423/1‑145 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCa                                                                      >  7:39917/1‑145 (MQ=255)
                                                                             aacaacGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTAGCCCAGTCCACCACGCCGGAATATGCa                                                                      >  5:134870/1‑145 (MQ=255)
                                                                                           gAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATaa                                                      >  4:122093/1‑147 (MQ=255)
                                                                                           gAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAAt                                                     >  3:93051/1‑148 (MQ=255)
                                                                                           gAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAAt                                                     >  8:277426/1‑148 (MQ=255)
                                                                                              gAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCAt                                                  >  4:274188/1‑148 (MQ=255)
                                                                                                aCCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCa                                                                                                          >  2:38783/1‑90 (MQ=255)
                                                                                                aCCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCa                                                                                                          <  1:38783/90‑1 (MQ=255)
                                                                                                 ccGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCAtttt                                               <  8:71492/148‑1 (MQ=255)
                                                                                                  cGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAATAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTcc                                                                                        >  1:171679/1‑106 (MQ=255)
                                                                                                  cGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAATAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTcc                                                                                        <  2:171679/106‑1 (MQ=255)
                                                                                                       tCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCtt                                         >  6:59867/1‑148 (MQ=255)
                                                                                                             aGACGCCAGCAGCAACCAGGCACTCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCTCCCAGTCCACCACGCCGGAATATGCAGAGAAACAACCGATATTCAGTTTAGGCTTGTGTTc                                   >  2:339174/1‑148 (MQ=255)
                                                                                                               aCGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCtt                                 <  6:91166/148‑1 (MQ=255)
                                                                                                                cGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCttt                                >  8:419674/1‑148 (MQ=255)
                                                                                                                cGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCttt                                <  8:191117/148‑1 (MQ=255)
                                                                                                                  ccAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACATGTAAGCACAGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATATTCGGCTTGTGTTCTTTgg                              <  5:159521/148‑1 (MQ=255)
                                                                                                                  ccAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCtt                                 <  8:155922/145‑1 (MQ=255)
                                                                                                                  ccAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCtt                                 >  7:155922/1‑145 (MQ=255)
                                                                                                                       agcaACCAGTCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACAGAAGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCTACCACGCCGGAATATGCCGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTgg                         <  6:270704/148‑2 (MQ=255)
                                                                                                                        gcaACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGtt                        <  5:59867/148‑1 (MQ=255)
                                                                                                                           aCCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGttttt                     <  2:326285/148‑1 (MQ=255)
                                                                                                                            ccAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTc                    >  2:140364/1‑148 (MQ=255)
                                                                                                                              aGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCt                                          >  6:131333/1‑124 (MQ=255)
                                                                                                                              aGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCt                                          <  5:131333/124‑1 (MQ=255)
                                                                                                                                gCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTc                    >  4:241915/1‑144 (MQ=255)
                                                                                                                                 cccGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTc                    >  6:163856/1‑143 (MQ=255)
                                                                                                                                    gCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCgg            <  3:241915/148‑1 (MQ=255)
                                                                                                                                         gTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGAACGGCGTAg       >  2:288061/1‑148 (MQ=255)
                                                                                                                                            cgcCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCg    <  3:122093/148‑1 (MQ=255)
                                                                                                                                            cgcCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCg    <  4:93051/148‑1 (MQ=255)
                                                                                                                                              ccATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGAATTCAAGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCAACGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGat  <  1:339174/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TACAGCTCTTCGCTACCACCTTTCGCCAGGATCAGGCCGCCGCGCGGACCCGCCAGGGTTTTGTGAGTGGTGGTAGTAACAACGTGAGCATGAGGAACCGGGTTCGGGTAGACGCCAGCAGCAACCAGGCCCGCAACGTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAATATGCAGAGAAACCACCGATAATCATTTTCGGCTTGTGTTCTTTGGCTTGTTTTTCCAGATCGGCGTAGTCGAT  >  NC_000913/2684762‑2685051

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