Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,080,616 C→A 5.8% G205C (GGC→TGC)  fdhE ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,6160CA 5.8% 619.4 / 21.1 171G205C (GGC→TGC) fdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base C (85/76);  new base A (6/4);  total (91/80)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 174 total aligned reads displayed.

GCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTA  >  NC_000913/4080474‑4080764
                                                                                                                                              |                                                                                                                                                    
gCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTg                                                                                                                                                <  2:787623/149‑1 (MQ=255)
 cTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGc                                                                                                                                               <  2:992293/149‑1 (MQ=255)
  tttAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa                                                                                                                                              >  4:781464/1‑149 (MQ=255)
  tttAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCa                                                                                                                                              >  3:721567/1‑149 (MQ=255)
   ttAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCAc                                                                                                                                             >  4:1223949/1‑149 (MQ=255)
     aaTCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCa                                                                                                                                           <  3:829958/149‑1 (MQ=255)
       tCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATg                                                                                                                                         >  1:767884/1‑149 (MQ=255)
       tCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCAAATTGGCACCac                                                                                                                                          >  1:1287218/1‑147 (MQ=255)
         gcgGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCt                                                                                                                                       <  1:1223945/149‑1 (MQ=255)
            gCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGa                                                                                                                                    <  2:297493/149‑1 (MQ=255)
              cTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGaca                                                                                                                                  >  3:813871/1‑149 (MQ=255)
               tGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGacac                                                                                                                                 >  4:300866/1‑149 (MQ=255)
                 ttCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGTCAATTTGCAc                                                                                                                                             <  4:1354351/135‑1 (MQ=255)
                 ttCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGTCAATTTGCAc                                                                                                                                             >  3:1354351/1‑135 (MQ=255)
                    gTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCa                                                                                                                            <  4:182469/149‑1 (MQ=255)
                     tCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCAt                                                                                                                           <  4:1105201/149‑1 (MQ=255)
                     tCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCAt                                                                                                                           >  2:921963/1‑149 (MQ=255)
                     tCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCAt                                                                                                                           <  3:1223949/149‑1 (MQ=255)
                       aTCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAg                                                                                                                         <  4:335182/149‑1 (MQ=255)
                        tCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGa                                                                                                                        >  2:263034/1‑149 (MQ=255)
                         ccAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGaa                                                                                                                       <  3:781464/149‑1 (MQ=255)
                          cAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAAc                                                                                                                      <  2:750022/149‑1 (MQ=255)
                           aGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGaa                                                                                                                       >  3:168528/1‑147 (MQ=255)
                                 cAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAAc                                                                                                               <  2:39170/149‑1 (MQ=255)
                                   gTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACca                                                                                                                    <  2:48680/142‑1 (MQ=255)
                                   gTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCACGCTGGACACCGGCATAGAACCACAAACCg                                                                                                             <  1:178880/149‑1 (MQ=255)
                                      aTGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGc                                                                                                          <  4:1157066/149‑1 (MQ=255)
                                       tGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCa                                                                                                         >  4:715356/1‑149 (MQ=255)
                                           gTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAatat                                                                                                     >  3:1252435/1‑149 (MQ=255)
                                            tttGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATAtt                                                                                                    >  3:768378/1‑149 (MQ=255)
                                            tttGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATAtt                                                                                                    >  2:719702/1‑149 (MQ=255)
                                                        cTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGtt                                                                                              <  1:823340/143‑1 (MQ=255)
                                                         tCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGcc                                                                                       >  3:392604/1‑149 (MQ=255)
                                                         tCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTACTCGcc                                                                                       >  1:1111188/1‑149 (MQ=255)
                                                           gCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGt                                                                                     >  4:745585/1‑149 (MQ=255)
                                                           gCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACATTGCTCGCCGt                                                                                     >  4:519193/1‑149 (MQ=255)
                                                                tgctgcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCa                                                                                <  1:51792/149‑1 (MQ=255)
                                                                  ctgcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAg                                                                                                                         <  2:328440/106‑1 (MQ=255)
                                                                   tgcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCgg                                                                                                                              >  4:419607/1‑100 (MQ=255)
                                                                   tgcATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCgg                                                                                                            <  3:273605/118‑1 (MQ=255)
                                                                       ttttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCggg                                                                         >  2:578292/1‑149 (MQ=255)
                                                                       ttttACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCACAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCggg                                                                         >  1:845062/1‑149 (MQ=255)
                                                                               gcACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGGGTGGTGCCAAttt                                                                                                                                                 <  4:875323/69‑1 (MQ=255)
                                                                               gcACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAAttt                                                                                                                                                 >  3:875323/1‑69 (MQ=255)
                                                                               gcACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAAttt                                                                                                                                                 <  2:1199604/69‑1 (MQ=255)
                                                                               gcACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACATTGCTCGCCGTATTCAGCGCGGGCTTTGCCg                                                                 <  3:519193/149‑1 (MQ=255)
                                                                                cACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgg                                                                >  3:350941/1‑149 (MQ=255)
                                                                                 aCTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCggg                                                               <  1:1418550/149‑1 (MQ=255)
                                                                                  cTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                              <  4:813871/149‑1 (MQ=255)
                                                                                     cGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCACAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATc                                                           <  2:845062/149‑1 (MQ=255)
                                                                                        gCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGGCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGa                                                        >  4:853951/1‑149 (MQ=255)
                                                                                        gCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGa                                                        >  3:11385/1‑149 (MQ=255)
                                                                                         ccATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAt                                                       >  1:586553/1‑149 (MQ=255)
                                                                                           aTTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTg                                                     >  4:16937/1‑149 (MQ=255)
                                                                                           aTTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTg                                                     <  1:921963/149‑1 (MQ=255)
                                                                                            ttCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTgg                                                    <  3:300866/149‑1 (MQ=255)
                                                                                             tCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGc                                                   <  1:719702/149‑1 (MQ=255)
                                                                                                 tttCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgg                                                                <  2:465878/132‑1 (MQ=255)
                                                                                                 tttCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgg                                                                >  1:465878/1‑132 (MQ=255)
                                                                                                 tttCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATc                                               >  4:75571/1‑149 (MQ=255)
                                                                                                   tCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATc                                                           >  1:664516/1‑135 (MQ=255)
                                                                                                   tCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATc                                                           <  2:664516/135‑1 (MQ=255)
                                                                                                     acacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGa                                                             <  3:684682/131‑1 (MQ=255)
                                                                                                     acacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGa                                                             >  4:684682/1‑131 (MQ=255)
                                                                                                     acacaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTgcg                                           >  1:393522/1‑149 (MQ=255)
                                                                                                       acaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTg                                                                                                   >  1:28228/1‑91 (MQ=255)
                                                                                                       acaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                              >  4:381843/1‑128 (MQ=255)
                                                                                                       acaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCgggg                                                              <  3:562089/128‑1 (MQ=255)
                                                                                                       acaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTg                                                     <  4:522004/137‑1 (MQ=255)
                                                                                                       acaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTg                                                     >  3:522004/1‑137 (MQ=255)
                                                                                                       acaGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTGGCACCATGCTGGAGACCGGCAGAGAACCACAAACCGGGCAATATTGACGTTGCGCGCCGTATTGAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGcc                                         >  3:911631/1‑149 (MQ=255)
                                                                                                          ggTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAg                                      <  3:745585/149‑1 (MQ=255)
                                                                                                            ttGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTa                                    >  4:1035715/1‑149 (MQ=255)
                                                                                                                       gTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAg                                      >  3:1749/1‑136 (MQ=255)
                                                                                                                       gTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAg                                      <  4:1749/136‑1 (MQ=255)
                                                                                                                       gTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTg                         >  1:147687/1‑149 (MQ=255)
                                                                                                                       gTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTg                         <  3:50355/149‑1 (MQ=255)
                                                                                                                          aCGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAg                      >  4:1329408/1‑149 (MQ=255)
                                                                                                                            gTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGcc                    >  1:906561/1‑149 (MQ=255)
                                                                                                                             tAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGccc                   >  4:986384/1‑149 (MQ=255)
                                                                                                                               gCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGtt                                                                                              <  4:364192/72‑1 (MQ=255)
                                                                                                                               gCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGtt                                                                                              >  3:364192/1‑72 (MQ=255)
                                                                                                                                  ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATa              >  4:1051845/1‑149 (MQ=255)
                                                                                                                                  ccTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATa              >  2:1299556/1‑149 (MQ=255)
                                                                                                                                   cagcatggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATaa             >  1:1141685/6‑149 (MQ=255)
                                                                                                                                   cTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATaa             >  4:1470188/1‑149 (MQ=255)
                                                                                                                                    tGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACACACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGCGATCAGATTGGCCATCTGCGCCCAGCAGTGCGACAGTGCAGCCCAGATaaa            <  4:911631/149‑1 (MQ=255)
                                                                                                                                     gcatggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTgg                                                    >  2:133817/4‑108 (MQ=255)
                                                                                                                                     gcatggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTgg                                                    <  1:133817/105‑1 (MQ=255)
                                                                                                                                     gcatggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAAc           >  3:991505/4‑149 (MQ=255)
                                                                                                                                     gAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATaaaa           <  4:1266474/149‑2 (MQ=255)
                                                                                                                                         ggtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGgcg       <  1:1352409/149‑1 (MQ=255)
                                                                                                                                          gtgCAAATTGGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGgcgc      >  2:1448865/1‑149 (MQ=255)
                                                                                                                                           tgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGt                                     <  2:450000/117‑1 (MQ=255)
                                                                                                                                           tgCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGt                                     >  1:450000/1‑117 (MQ=255)
                                                                                                                                            gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGcc                                                                                       >  3:1050697/1‑66 (MQ=255)
                                                                                                                                            gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGAtt                                                      <  4:538555/99‑1 (MQ=255)
                                                                                                                                            gCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTACTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCtt    <  2:1111188/149‑1 (MQ=255)
                                                                                                                                              cAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTgg                                                    <  1:332256/99‑1 (MQ=255)
                                                                                                                                              cAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTa  >  2:1248400/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GCTTTAATCGCGGCCTGTTCGTCATCCAGCGACCAGTAATGCAGTTTGCCGCTCTGCTCGCAGTTGCTGCATTTTACGCGCACTACGTGCCATTCAGTTTCACACAGGTTGCAGTGCAGGTAACGTAGCCCCTGAGTGGTGCCAATTTGCACCATGCTGGACACCGGCATAGAACCACAAACCGGGCAATATTGACGTTGCTCGCCGTATTCAGCGCGGGCTTTGCCGGGGATCAGATTGGCCATCTGCGCCCAGTAGAGCGACAGTGCAGCCCAGATAAACGGCGCTTTA  >  NC_000913/4080474‑4080764

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