Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,130,060 C→T 15.1% Q76* (CAA→TAA)  metL → Bifunctional aspartokinase/homoserine dehydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,130,0600CT15.1% 2217.1 / 302.6 838Q76* (CAA→TAA) metLBifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (396/314);  new base T (72/55);  total (468/370)
Fisher's exact test for biased strand distribution p-value = 9.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 887 total aligned reads displayed.

GTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAGT  >  NC_000913/4129912‑4130208
                                                                                                                                                    |                                                                                                                                                    
 tgtTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTcaaca                                                                                                                                                     >  2:2883764/1‑149 (MQ=255)
   ttATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATaaa                                                                                                                                                   >  2:42366/1‑149 (MQ=255)
         gCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGc                                                                                                                                             <  4:3349354/149‑1 (MQ=255)
         gCGTGTCGAGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGc                                                                                                                                             >  4:616974/1‑149 (MQ=255)
                  gggCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATc                                                                                                                                    >  1:335180/1‑149 (MQ=255)
                  gggCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATc                                                                                                                                    >  4:5174249/1‑149 (MQ=255)
                  gggCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATc                                                                                                                                    >  2:169105/1‑149 (MQ=255)
                     cATTATGGCGGAGTACTCTTAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGt                                                                                                                                 <  4:2206784/149‑1 (MQ=255)
                       ttATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGc                                                                                                                               >  4:551402/1‑149 (MQ=255)
                          tGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGAt                                                                                                                            <  4:5836641/149‑1 (MQ=255)
                              ggAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGa                                                                                                                        <  2:4866314/149‑1 (MQ=255)
                              ggAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGa                                                                                                                        <  1:2597029/149‑1 (MQ=255)
                                 gTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTa                                                                                                                     <  4:2481075/149‑1 (MQ=255)
                                  tACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAg                                                                                                                    >  3:5819884/1‑149 (MQ=255)
                                   aCTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGc                                                                                                                   <  1:5426974/149‑1 (MQ=255)
                                     tctcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGTTAGCACCACTAACCAGTTGATTACCTGGTTGAAACTAAGCCATACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGGCGCTATCATTGCGATCTGATTAGCgg                                                                                                                 <  2:1819772/149‑1 (MQ=255)
                                     tctcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTa                                                                                                                     >  2:5049004/1‑145 (MQ=255)
                                     tctcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCgg                                                                                                                 >  1:4924140/1‑149 (MQ=255)
                                     tctcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCgg                                                                                                                 >  3:4347605/1‑149 (MQ=255)
                                       tcAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTc                                                                                                               >  3:3141594/1‑149 (MQ=255)
                                         aGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTg                                                                                                             <  2:1109266/149‑1 (MQ=255)
                                         aGCCTGACGATATGATGGCGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTg                                                                                                             >  1:2360773/1‑149 (MQ=255)
                                          gCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGc                                                                                                            <  2:3372085/149‑1 (MQ=255)
                                          gCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGc                                                                                                            >  4:6665926/1‑149 (MQ=255)
                                          gCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGc                                                                                                            >  1:6291014/1‑149 (MQ=255)
                                           ccTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCt                                                                                                           >  1:1290311/1‑149 (MQ=255)
                                               aCGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTAccc                                                                                                       >  1:5581209/1‑149 (MQ=255)
                                                   tatGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTg                                                                                                   <  3:2164784/149‑1 (MQ=255)
                                                   tatGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTg                                                                                                   <  3:160825/149‑1 (MQ=255)
                                                     tgatgGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTgaa                                                                                                 <  1:1083338/149‑1 (MQ=255)
                                                        tggtggTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTgaagaa                                                                                              <  3:6194865/149‑1 (MQ=255)
                                                              tttCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGc                                                                                                     >  1:7121038/1‑136 (MQ=255)
                                                              tttCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGAt                                                                                        >  1:1914940/1‑149 (MQ=255)
                                                               ttCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATa                                                                                       <  3:396499/149‑1 (MQ=255)
                                                               ttCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATa                                                                                       <  2:5668036/149‑1 (MQ=255)
                                                               ttCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATa                                                                                       <  4:4799560/149‑1 (MQ=255)
                                                                 ccgccgccgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGc                                                                                     <  1:3520732/149‑1 (MQ=255)
                                                                   gccgccgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCt                                                                                   >  2:1327753/1‑149 (MQ=255)
                                                                        cgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTa                                                                              <  3:5489205/149‑1 (MQ=255)
                                                                        cgGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTa                                                                              <  3:5265534/149‑1 (MQ=255)
                                                                          gTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAgc                                                                            <  2:6977905/149‑1 (MQ=255)
                                                                          gTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAgc                                                                            >  2:4506087/1‑149 (MQ=255)
                                                                               accacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCttt                                                                       >  2:5375362/1‑149 (MQ=255)
                                                                                 cacTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCa                                                                                 >  1:2783044/1‑137 (MQ=255)
                                                                                   cTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGAt                                                                                                                            <  2:6150/92‑1 (MQ=255)
                                                                                    tAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCa                                                                  <  1:4566597/149‑1 (MQ=255)
                                                                                    tAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCa                                                                  <  4:3972135/149‑1 (MQ=255)
                                                                                    tAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCa                                                                  >  2:2117325/1‑149 (MQ=255)
                                                                                      aCCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGc                                                                >  4:1577487/1‑149 (MQ=255)
                                                                                      aCCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGc                                                                <  1:6835581/149‑1 (MQ=255)
                                                                                       ccAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCg                                                               <  4:6809057/149‑1 (MQ=255)
                                                                                        cAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGa                                                              >  3:475654/1‑149 (MQ=255)
                                                                                            tGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCtt                                                          <  2:756267/149‑1 (MQ=255)
                                                                                              aTTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTgaagaa                                                                                              <  3:6803261/111‑1 (MQ=255)
                                                                                               ttAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAg                                                       >  3:6519083/1‑149 (MQ=255)
                                                                                                  aCTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAgcgc                                                    >  4:1834093/1‑149 (MQ=255)
                                                                                                  aCTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAgcgc                                                    >  3:2818833/1‑149 (MQ=255)
                                                                                                   cTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCAt                                                                                <  1:1701637/120‑1 (MQ=255)
                                                                                                      gTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTgg                                                >  2:2992655/1‑149 (MQ=255)
                                                                                                      gTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTgg                                                >  2:194753/1‑149 (MQ=255)
                                                                                                       ttGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTggc                                               >  4:1707422/1‑149 (MQ=255)
                                                                                                         gAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAATGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCAGTAGCGCTTTTGTCAGCGACCTTGAGCGCCTggcgg                                             >  3:4874205/1‑149 (MQ=255)
                                                                                                            atcaaGCCAGACCGATCGTATCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGgcgc                                          >  2:2278097/4‑149 (MQ=255)
                                                                                                            aCTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGgcgc                                          >  2:3579997/1‑149 (MQ=255)
                                                                                                             cTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCtt                                                                        <  1:421741/118‑1 (MQ=255)
                                                                                                              tAAGCCATACCGATCGTCTCTCTGAGCATCAGGTTCAACAAACGCGGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTCCCCGGTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCGGCGACCTTCACCGCCTGGCGGCgctg                                        >  4:6459635/1‑149 (MQ=255)
                                                                                                              tAAGCCAGACCGATCGTATCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCgctg                                        >  1:4376873/1‑149 (MQ=255)
                                                                                                               aaGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCgctgc                                       >  4:742878/1‑149 (MQ=255)
                                                                                                               aaGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCgctgc                                       <  1:1536118/149‑1 (MQ=255)
                                                                                                                aGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTga                                                                                                  <  3:4171776/89‑1 (MQ=255)
                                                                                                                 gCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTc                                     <  1:3834991/149‑1 (MQ=255)
                                                                                                                  ccAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATc                                                                                                                           <  4:31663/62‑1 (MQ=255)
                                                                                                                  ccAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCg                                    >  1:6577681/1‑149 (MQ=255)
                                                                                                                    aGACCGATCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGAc                                  >  3:4527599/1‑149 (MQ=255)
                                                                                                                    aGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGAc                                  >  2:2420689/1‑149 (MQ=255)
                                                                                                                        cGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCg                              >  2:6408935/1‑149 (MQ=255)
                                                                                                                        cGATCGTCTCTCTGCGCAGCAGGTTCAACAAACGCTGCTTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAACCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCGCGACAGCg                              <  2:5763964/149‑1 (MQ=255)
                                                                                                                          aTCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGt                            >  2:3941492/1‑149 (MQ=255)
                                                                                                                          aTCGTCTCTCTGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGt                            <  4:4527599/149‑1 (MQ=255)
                                                                                                                          aTCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGt                            <  4:4209346/149‑1 (MQ=255)
                                                                                                                             gTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTAtt                         >  3:1417103/1‑149 (MQ=255)
                                                                                                                                tctctGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAac                      <  4:1131844/149‑1 (MQ=255)
                                                                                                                                tctctGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAac                      <  1:2167091/149‑1 (MQ=255)
                                                                                                                                   ctGCGCATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAacgac                   >  4:1970501/1‑149 (MQ=255)
                                                                                                                                      cgcATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCAATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTGTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAGCGACGCa                <  2:659697/149‑1 (MQ=255)
                                                                                                                                       gcATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGcc                                                   >  3:3540882/1‑113 (MQ=255)
                                                                                                                                        cATCAGGTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGc                                                                                                                   >  2:2685427/1‑48 (MQ=255)
                                                                                                                                          tCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAgtgt            <  3:713460/149‑1 (MQ=255)
                                                                                                                                          tCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAgtgt            <  3:1918113/149‑1 (MQ=255)
                                                                                                                                              gTTCAATAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGc        >  1:3209831/1‑149 (MQ=255)
                                                                                                                                               ttCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAacgac                   <  3:244056/137‑1 (MQ=255)
                                                                                                                                               ttCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCg       >  4:4946986/1‑149 (MQ=255)
                                                                                                                                                  aacaaACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGcc                                                   <  4:2692386/102‑1 (MQ=255)
                                                                                                                                                  aacaaACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGaa    <  1:773092/149‑1 (MQ=255)
                                                                                                                                                  aacaaACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGaa    >  1:5138325/1‑149 (MQ=255)
                                                                                                                                                   acaaACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAg   >  4:3677833/1‑149 (MQ=255)
                                                                                                                                                   acaaACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAg   <  3:1096274/149‑1 (MQ=255)
                                                                                                                                                   acaaACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAg   >  1:2901342/1‑149 (MQ=255)
                                                                                                                                                    caaACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAgt  >  4:2994699/1‑149 (MQ=255)
                                                                                                                                                    caaACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAgt  >  1:5787282/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGTTATTTGCGTGTCGCGGGCATTATGGCGGAGTACTCTCAGCCTGACGATATGATGGTGGTTTCCGCCGCCGGTAGCACCACTAACCAGTTGATTAACTGGTTGAAACTAAGCCAGACCGATCGTCTCTCTGCGCATCAGGTTCAACAAACGCTGCGTCGCTATCAGTGCGATCTGATTAGCGGTCTGCTACCCGCTGAAGAAGCCGATAGCCTCATTAGCGCTTTTGTCAGCGACCTTGAGCGCCTGGCGGCGCTGCTCGACAGCGGTATTAACGACGCAGTGTATGCGGAAGT  >  NC_000913/4129912‑4130208

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