Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,068,370 T→G intergenic (+116/+141) ymdF → / ← rutG KGG family protein/pyrimidine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,068,3700TG99.5% 643.3 / NA 182intergenic (+116/+141)ymdF/rutGKGG family protein/pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base G (96/84);  minor base A (0/1);  total (96/86)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

AATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTAC  >  NC_000913/1068273‑1068466
                                                                                                 |                                                                                                
aaTGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTa                                                                                                <  2:686384/100‑1 (MQ=255)
aaTGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTa                                                                                                <  2:2400331/100‑1 (MQ=255)
aaTGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTa                                                                                                <  2:4735944/100‑1 (MQ=255)
 aTGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTAt                                                                                               >  1:4543780/1‑100 (MQ=255)
  tGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATc                                                                                              <  2:653543/100‑1 (MQ=255)
   gCATTTTTGTCTGTAGCTCGTCAAAAAGCCTTCGCCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCa                                                                                             <  2:4207941/100‑1 (MQ=255)
   gCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCa                                                                                             <  2:2174333/100‑1 (MQ=255)
    cATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAg                                                                                            >  1:415234/1‑100 (MQ=255)
    cATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTGACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAg                                                                                            <  1:161049/100‑1 (MQ=255)
     aTTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAgg                                                                                           >  1:5528714/1‑100 (MQ=255)
      tttttGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGc                                                                                          <  2:969310/100‑1 (MQ=255)
      tttttGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGc                                                                                          <  1:293008/100‑1 (MQ=255)
      gttttGTGTGTAGTTCGTCAAAAACCCATCACCGCAGGTGACCCGGTGGTTGATACTGATGATAAATATAAGGATGCCTGCTGCGCTATGCGTATCGGGc                                                                                          <  2:3733363/99‑1 (MQ=255)
           gTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTAcc                                                                                     <  1:5452452/100‑1 (MQ=255)
           gTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTAcc                                                                                     >  2:3803922/1‑100 (MQ=255)
            tCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCa                                                                                    <  2:2002174/100‑1 (MQ=255)
                 aGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCagaaga                                                                               >  1:5395417/1‑100 (MQ=255)
                    tCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTg                                                                            >  1:264477/1‑100 (MQ=255)
                    tCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTg                                                                            <  1:4739756/100‑1 (MQ=255)
                     cGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGc                                                                           >  2:2115397/1‑100 (MQ=255)
                     cGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGc                                                                           <  1:5459034/100‑1 (MQ=255)
                      gTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCa                                                                          >  2:3518344/1‑100 (MQ=255)
                      gTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCa                                                                          <  2:1582278/100‑1 (MQ=255)
                      gTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCa                                                                          >  1:4975738/1‑100 (MQ=255)
                       tCAAAAAGCCATCACCGCCGGTTAGCCGGCGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCaa                                                                         >  1:5189904/1‑100 (MQ=255)
                       tCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCaa                                                                         <  2:2727278/100‑1 (MQ=255)
                        cAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAat                                                                        <  1:1274876/100‑1 (MQ=255)
                        cAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAat                                                                        >  2:3684107/1‑100 (MQ=255)
                         aaaaaGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAata                                                                       >  1:2277308/1‑100 (MQ=255)
                          aaaaGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAatat                                                                      >  2:4872766/1‑100 (MQ=255)
                          aaaaGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAatat                                                                      >  2:4872825/1‑100 (MQ=255)
                               ccATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGa                                                                 <  2:4475814/100‑1 (MQ=255)
                               ccATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGa                                                                 >  1:1788855/1‑100 (MQ=255)
                               ccATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGa                                                                 >  2:3971198/1‑100 (MQ=255)
                                cATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGaa                                                                >  1:1155855/1‑100 (MQ=255)
                                 aTCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAAt                                                               <  2:2529853/100‑1 (MQ=255)
                                  tCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAAtt                                                              <  2:1386031/100‑1 (MQ=255)
                                   cACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAAttt                                                             >  1:31898/1‑100 (MQ=255)
                                    acggtcgGTTACCGGCGGGTTGCTACTGATGACAAAGGTAGGCGTGCGTGAGGCGCGATGCGTAGCAGGCCTACCAGAAGATTGCAATATATTGAATTTg                                                            <  2:4474188/95‑1 (MQ=19)
                                     ccgccgGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGc                                                           >  1:635668/1‑100 (MQ=255)
                                       gccgGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCAc                                                         <  2:2237844/100‑1 (MQ=255)
                                       gccgGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCAc                                                         <  2:5322713/100‑1 (MQ=255)
                                       gccgGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCAc                                                         >  1:2817106/1‑100 (MQ=255)
                                       gccgGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCAc                                                         <  2:3281851/100‑1 (MQ=255)
                                         cgGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTg                                                       >  1:3731184/1‑100 (MQ=255)
                                          ggTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGCCATACCAGAAGATTGCAATATATTGAATTTGCACTGt                                                      <  2:793308/100‑1 (MQ=255)
                                               cccGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGt                                                 >  1:484191/1‑100 (MQ=255)
                                               cccGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTCATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGt                                                 <  2:2450104/100‑1 (MQ=255)
                                                 cGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAg                                               <  2:1834011/100‑1 (MQ=255)
                                                  ggtggtTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAgg                                              <  2:4280458/100‑1 (MQ=255)
                                                   gtggtTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGc                                             >  1:2058956/1‑100 (MQ=255)
                                                     ggtTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCg                                           >  1:3739151/1‑100 (MQ=255)
                                                      gtTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCgg                                          <  1:709540/100‑1 (MQ=255)
                                                      gtTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCgg                                          <  2:2084783/100‑1 (MQ=255)
                                                        tGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATCTTGAATTTGCACTGTTTTGTAGGCCggct                                        >  1:5089904/1‑98 (MQ=255)
                                                          aTACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATaa                                      <  1:3717280/100‑1 (MQ=255)
                                                           tACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAg                                     >  2:1958333/1‑100 (MQ=255)
                                                             cTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGc                                   >  2:1177083/1‑100 (MQ=255)
                                                             cTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGc                                   >  1:3805266/1‑100 (MQ=255)
                                                              tgatgaCAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCg                                  >  1:4070377/1‑100 (MQ=255)
                                                              tgatgaCAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCg                                  <  2:263662/100‑1 (MQ=255)
                                                               gatgaCAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGt                                 >  1:4381804/1‑100 (MQ=255)
                                                                 tgaCAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGttt                               >  2:4964284/1‑100 (MQ=255)
                                                                 tgaCAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGttt                               >  1:1056879/1‑100 (MQ=255)
                                                                 tgaCAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGttt                               >  1:3979624/1‑100 (MQ=255)
                                                                  gaCAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTa                              >  2:1855514/1‑100 (MQ=255)
                                                                  gaCAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTa                              >  1:3878663/1‑100 (MQ=255)
                                                                   aCAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGAGTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAc                             <  1:5292319/100‑1 (MQ=255)
                                                                    cAAATGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcg                            >  1:524074/1‑100 (MQ=255)
                                                                     aaaTGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCGGGATAAGGCGTTTAcgc                           <  2:2657493/100‑1 (MQ=255)
                                                                      aaTGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCGATACATTGAATTTGCACTGTTTTGTAGGCCGGATAACGCGTTTAcgcc                          <  1:4887148/100‑1 (MQ=25)
                                                                      aaTGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcgcc                          <  1:1816851/100‑1 (MQ=255)
                                                                      aaTGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcgcc                          <  1:3551183/100‑1 (MQ=255)
                                                                       aTGTAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcgccg                         <  2:5332935/100‑1 (MQ=255)
                                                                          tAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                      >  1:1033070/1‑100 (MQ=255)
                                                                          tAAGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCAt                      >  2:1571724/1‑100 (MQ=255)
                                                                            aGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                    <  1:2225135/100‑1 (MQ=255)
                                                                            aGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                    <  1:3044937/100‑1 (MQ=255)
                                                                            aGCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                    >  2:1738186/1‑100 (MQ=255)
                                                                             gCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                   >  2:318020/1‑100 (MQ=255)
                                                                             gCTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                   >  2:5088733/1‑100 (MQ=255)
                                                                             gCTTGCCTGAAGCGCGATGCGTATCAGGCCTACCAGAATATTACAATATAATGGATTTGCACTGTTGTGTAGCCGGGATAAGGCGCTTACTCCGCATCCg                   >  2:835327/1‑100 (MQ=21)
                                                                              cTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                  >  2:464443/1‑100 (MQ=255)
                                                                              cTTGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                  >  1:2529273/1‑100 (MQ=255)
                                                                              cTTGCCTGATGCGCGATGCGTATCAGGCATACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCTGGATATGGCGTTTACGCCGCAttcag                  >  2:145431/1‑96 (MQ=18)
                                                                               ttGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                 <  1:3015801/100‑1 (MQ=255)
                                                                               ttGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                 <  1:2137838/100‑1 (MQ=255)
                                                                                tGCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                >  1:3281317/1‑100 (MQ=255)
                                                                                 gCCTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAt               >  2:5191746/1‑100 (MQ=255)
                                                                                  ccTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGTAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATg              <  1:5330913/100‑1 (MQ=255)
                                                                                  ccTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATg              <  2:3770349/100‑1 (MQ=255)
                                                                                   cTGATGCGCGATGCGTATCAGGCCTACCGGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGa             <  1:5058774/100‑1 (MQ=255)
                                                                                   cTGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGa             >  2:3392293/1‑100 (MQ=255)
                                                                                    tGATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTTGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaa            >  2:3017442/1‑100 (MQ=255)
                                                                                     gATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAc           <  2:480535/100‑1 (MQ=255)
                                                                                     gATGCGCGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAc           <  1:813107/100‑1 (MQ=255)
                                                                                        gcgcGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAAt        <  2:3441257/100‑1 (MQ=255)
                                                                                        gcgcGATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAAt        <  1:3077478/100‑1 (MQ=255)
                                                                                         cgcgATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATg       <  2:4103361/100‑1 (MQ=255)
                                                                                          gcgATGCGTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGc      <  1:4299497/100‑1 (MQ=255)
                                                                                                 |                                                                                                
AATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTAC  >  NC_000913/1068273‑1068466

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