Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 467,350 A→C intergenic (‑38/‑62) ybaE ← / → cof putative ABC transporter periplasmic binding protein/thiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913467,3500AC99.4% 640.9 / ‑5.4 178intergenic (‑38/‑62)ybaE/cofputative ABC transporter periplasmic binding protein/thiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (99/78);  total (99/79)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 179 total aligned reads displayed.

TTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCAC  >  NC_000913/467254‑467446
                                                                                                |                                                                                                
tttCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCtt                                                                                               >  2:2042589/1‑100 (MQ=255)
 ttCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTc                                                                                              <  1:1701764/100‑1 (MQ=255)
  tCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTAttccttcc                                                                                             <  2:3735961/100‑1 (MQ=255)
   ccGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCg                                                                                            <  1:2902676/100‑1 (MQ=255)
   ccGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCg                                                                                            <  2:480026/100‑1 (MQ=255)
    cggcggAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGt                                                                                           <  2:2200618/100‑1 (MQ=255)
     ggcggAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGTATATAAACCCGGAGCAATATTATTTAACTAGTCACTATTCCTTCCGta                                                                                          <  2:3340222/100‑1 (MQ=255)
       cggAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGtata                                                                                        <  1:74637/100‑1 (MQ=255)
         gAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGtatata                                                                                      >  1:4322966/1‑100 (MQ=255)
         gAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGtatata                                                                                      >  1:3244810/1‑100 (MQ=255)
           aGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATc                                                                                    >  1:5555657/1‑100 (MQ=255)
            ggTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCa                                                                                   <  1:3294847/100‑1 (MQ=255)
            ggTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCa                                                                                   >  2:32125/1‑100 (MQ=255)
             gTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAg                                                                                  <  2:1798521/100‑1 (MQ=255)
               tGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGt                                                                                <  1:4001378/100‑1 (MQ=255)
                gCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTg                                                                               >  2:3112714/1‑100 (MQ=255)
                  cacaGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGAt                                                                             <  1:3679733/100‑1 (MQ=255)
                    caGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATAc                                                                           >  2:3461796/1‑100 (MQ=255)
                     aGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACt                                                                          <  1:1104824/100‑1 (MQ=255)
                     aGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACt                                                                          >  2:2512851/1‑100 (MQ=255)
                      gACGTTGATACTGGTTAAGACGGTTGTGAAATCGCATATAATCCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCACGTGATACTc                                                                         <  2:4734760/100‑1 (MQ=255)
                      gACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTc                                                                         <  2:4806201/100‑1 (MQ=255)
                        cGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCaa                                                                       >  1:2183221/1‑100 (MQ=255)
                         gTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAAt                                                                      >  1:3355425/1‑100 (MQ=255)
                          ttGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATc                                                                     >  1:3876733/1‑100 (MQ=255)
                          ttGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATc                                                                     <  1:2268368/100‑1 (MQ=255)
                          ttGATACTGGTTAAGACGGTTGAGCAAGCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATc                                                                     <  1:3041262/100‑1 (MQ=255)
                            gATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCAc                                                                   >  2:3635469/1‑100 (MQ=255)
                              tACTGGTTAAGACGGTTGAGCGATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCa                                                                 >  1:5393595/1‑100 (MQ=255)
                              tACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCa                                                                 >  2:1452945/1‑100 (MQ=255)
                               aCTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCAt                                                                >  2:2777347/1‑100 (MQ=255)
                                cTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCAtt                                                               >  1:479214/1‑100 (MQ=255)
                                 tGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGACACTCAATCACCATTa                                                              <  2:4174415/100‑1 (MQ=255)
                                  ggTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTaa                                                             >  1:3722557/1‑100 (MQ=255)
                                  ggTTAAGACGGTTGAGCAATCACATATAAACCCGGAACAATATTATTTAACTATCCACTATTCCTTGCGTATATATCAGGTGATACTCAATCACTATTaa                                                             >  2:1325511/1‑100 (MQ=255)
                                   gTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAAGCACCATTAAc                                                            >  2:770124/1‑100 (MQ=255)
                                    ttAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAAcc                                                           <  2:3368803/100‑1 (MQ=255)
                                    ttAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAAcc                                                           <  2:2105556/100‑1 (MQ=255)
                                    ttAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAAcc                                                           <  2:2896672/100‑1 (MQ=255)
                                    ttAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAAcc                                                           <  2:2763846/100‑1 (MQ=255)
                                     tAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCg                                                          >  2:2053327/1‑100 (MQ=255)
                                      aaGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCgt                                                         >  2:4861529/1‑100 (MQ=255)
                                       ngACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCgtg                                                        >  1:4421469/2‑100 (MQ=255)
                                       aGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCgtg                                                        >  1:5003591/1‑100 (MQ=255)
                                       aGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCgtg                                                        >  1:1135533/1‑100 (MQ=255)
                                        gACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCgtgt                                                       >  1:4139159/1‑100 (MQ=255)
                                         aCGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTc                                                      <  1:2663023/100‑1 (MQ=255)
                                          nggTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTca                                                     >  1:1507554/2‑100 (MQ=255)
                                           ggTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTcac                                                    >  1:4852436/1‑100 (MQ=255)
                                             ttGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACag                                                  >  2:44844/1‑100 (MQ=255)
                                                aGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGt                                               <  1:4574803/100‑1 (MQ=255)
                                                 gCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTg                                              <  2:3942791/100‑1 (MQ=255)
                                                 gCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTg                                              >  2:545223/1‑100 (MQ=255)
                                                  cAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTgg                                             <  1:887829/100‑1 (MQ=255)
                                                  cAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTgg                                             <  1:3553408/100‑1 (MQ=255)
                                                   aaTCGCATATAAACCCGTAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGga                                            <  1:3140787/100‑1 (MQ=255)
                                                    aTCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGgag                                           >  1:1759310/1‑100 (MQ=255)
                                                    aTCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGgag                                           >  1:1558270/1‑100 (MQ=255)
                                                    aTCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGgag                                           >  2:3148472/1‑100 (MQ=255)
                                                     tCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGgaga                                          >  2:4340952/1‑100 (MQ=255)
                                                     tCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGgaga                                          <  2:4283180/100‑1 (MQ=255)
                                                      cGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGCGTGGacaa                                         >  2:1228173/1‑100 (MQ=255)
                                                      cGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaa                                         <  2:4347804/100‑1 (MQ=255)
                                                      cGCATATAAACCCGGAACAATATTATTTAACTATCCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaa                                         >  2:5372448/1‑100 (MQ=255)
                                                       gCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaaa                                        >  1:49075/1‑100 (MQ=255)
                                                       gCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaaa                                        >  2:2295251/1‑100 (MQ=255)
                                                       gCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaaa                                        <  1:338416/100‑1 (MQ=255)
                                                       gCATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaaa                                        >  1:1339522/1‑100 (MQ=255)
                                                        cATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaaag                                       <  2:5192092/100‑1 (MQ=255)
                                                        cATATAAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaaag                                       >  2:4194276/1‑100 (MQ=255)
                                                          tataAGCCCGGAACAATATTATTTAACTACTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaaagaa                                     <  1:5082915/100‑1 (MQ=255)
                                                          tataAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaaagaa                                     >  2:2413214/1‑100 (MQ=255)
                                                          tataAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGagaaagaa                                     >  1:2034615/1‑100 (MQ=255)
                                                           ataAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTATCACAGAGTGGAgaaagaaa                                    <  2:1120175/100‑1 (MQ=255)
                                                            taAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAAt                                   <  2:1514322/100‑1 (MQ=255)
                                                            taAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAAt                                   >  1:1783218/1‑100 (MQ=255)
                                                            taAACCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAAt                                   >  1:2030962/1‑100 (MQ=255)
                                                             aaaCCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATg                                  >  2:3618204/1‑100 (MQ=255)
                                                             aaaCCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATg                                  <  2:3795354/100‑1 (MQ=255)
                                                               aCCCGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGc                                >  1:5235391/1‑100 (MQ=255)
                                                                cccGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCt                               >  1:2981928/1‑100 (MQ=255)
                                                                  cGGAACAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCg                             <  2:878461/100‑1 (MQ=255)
                                                                       cAATATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTgg                        >  2:1315359/1‑100 (MQ=255)
                                                                        aaTATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGgc                       <  2:1538954/100‑1 (MQ=255)
                                                                        aaTATTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGgc                       <  1:3494191/100‑1 (MQ=255)
                                                                         atatTATTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGgca                      <  1:3914115/100‑1 (MQ=255)
                                                                          tattatTTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGgcag                     >  2:1917733/1‑100 (MQ=255)
                                                                              attTAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCAtt                 <  1:3218315/100‑1 (MQ=255)
                                                                               tttAACTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCAttt                <  2:1120441/100‑1 (MQ=255)
                                                                                  aaCTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGat             <  2:1902050/100‑1 (MQ=255)
                                                                                  aaCTATTCACTATTCCGTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGat             <  1:1669232/100‑1 (MQ=255)
                                                                                    cTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGatat           >  1:2365701/1‑100 (MQ=255)
                                                                                    cTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGatat           >  2:3580766/1‑100 (MQ=255)
                                                                                    cTATTCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGatat           >  2:5342572/1‑100 (MQ=255)
                                                                                       ttCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATatgga        <  1:973138/100‑1 (MQ=255)
                                                                                        tCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATatggat       <  1:3213123/100‑1 (MQ=255)
                                                                                        tCACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATatggat       <  2:2701638/100‑1 (MQ=255)
                                                                                         cACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGTAAGAAATGGCGCGTCTGGCAGCATTTGACatggatg      >  2:5350909/1‑100 (MQ=255)
                                                                                         cACTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATatggatg      >  2:1653561/1‑100 (MQ=255)
                                                                                           cTATTCCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGc    <  2:4574357/100‑1 (MQ=255)
                                                                                                |                                                                                                
TTTCCGGCGGAAGGTTGCCACAGACGTTGATACTGGTTAAGACGGTTGAGCAATCGCATATAAACCCGGAACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCAC  >  NC_000913/467254‑467446

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