Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 145,691 C→T A204V (GCC→GTC)  yadE → putative polysaccharide deacetylase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913145,6910CT98.9% 307.4 / NA 93A204V (GCC→GTC) yadEputative polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (51/41);  minor base G (1/0);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGAT  >  NC_000913/145595‑145787
                                                                                                |                                                                                                
tGCTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTg                                                                                               <  2:631998/100‑1 (MQ=255)
 gCTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTgg                                                                                              >  2:1088994/1‑100 (MQ=255)
 gCTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTgg                                                                                              >  1:1600728/1‑100 (MQ=255)
      ccTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCg                                                                                         <  1:1407795/100‑1 (MQ=255)
      ccTATCACCATATTCTGCCCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCg                                                                                         <  2:1681181/100‑1 (MQ=255)
       cTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGt                                                                                        <  2:967329/100‑1 (MQ=255)
        tATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTg                                                                                       <  1:315653/100‑1 (MQ=255)
        tATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTg                                                                                       <  2:1506097/100‑1 (MQ=255)
        tATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGACGGTCTGGCTGCGTg                                                                                       >  2:695524/1‑100 (MQ=255)
          tCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGAc                                                                                     <  1:1204370/100‑1 (MQ=255)
           cACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACa                                                                                    >  1:2220504/1‑100 (MQ=255)
           cACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACa                                                                                    <  1:113667/100‑1 (MQ=255)
            aCCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAg                                                                                   >  1:1953730/1‑100 (MQ=255)
             ccATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGCCAgg                                                                                  >  2:1612107/1‑100 (MQ=255)
             ccATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAgg                                                                                  >  1:636755/1‑100 (MQ=255)
              cagatTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGGCTGGCTGCGTGACAggg                                                                                 >  2:933677/4‑100 (MQ=255)
                 atTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGAt                                                                              >  2:2077176/1‑100 (MQ=255)
                     tGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATAcgc                                                                          <  2:327697/100‑1 (MQ=255)
                     tGCGCGACGAAGAAAACACCCGTTTCCGCCATACTTCGACGACCACATCGGTACGCGCTCTCACTACCCAGATGGTCTGGCTGCGTGCCAGGGGATAcgc                                                                          >  2:866418/1‑100 (MQ=255)
                      gcgcGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATAcgcg                                                                         <  1:1683898/100‑1 (MQ=255)
                      gcgcGACGAAGAAAACACCCGTTTTCGCCATACTTCCACCACCACATCCGTACGCGCTTTCAATAACCAGATGTTCTGGCTGCGTGACAGGTGATAcgcg                                                                         <  2:1167556/100‑1 (MQ=255)
                       cgcgACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGa                                                                        >  1:1305112/1‑100 (MQ=255)
                       cgcgACGAAGAAAACACACGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGa                                                                        >  2:2065848/1‑100 (MQ=255)
                          gacgaAGACAGCACCCGGTTTCGCCATACATCGACGACCACATGGGTACGCGCTTTCAATAACCAGATGGTATGGCTGCGTGACAGGGGATAAGCGacac                                                                     <  2:1202279/100‑1 (MQ=255)
                          gacgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGacac                                                                     <  2:1196814/100‑1 (MQ=255)
                          gacgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGacac                                                                     >  2:801935/1‑100 (MQ=255)
                          gacgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGacac                                                                     >  2:2134212/1‑100 (MQ=255)
                           acgaAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACt                                                                    >  2:2063706/1‑100 (MQ=255)
                              aagaaAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAg                                                                 <  1:1611849/100‑1 (MQ=255)
                              aagaaAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAg                                                                 <  2:2140139/100‑1 (MQ=255)
                                gaaAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCa                                                               <  1:1804498/100‑1 (MQ=255)
                                   aaCACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATgg                                                            <  1:2089624/100‑1 (MQ=255)
                                    acacCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGt                                                           <  2:1529672/100‑1 (MQ=255)
                                     cacCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTg                                                          <  1:2194238/100‑1 (MQ=255)
                                     cacCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTg                                                          >  1:1450544/1‑100 (MQ=255)
                                         cGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGc                                                      >  2:276689/1‑100 (MQ=255)
                                           ttttCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTg                                                    <  1:1005858/100‑1 (MQ=255)
                                             ttCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGa                                                  >  1:1935771/1‑100 (MQ=255)
                                               cGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAg                                                <  1:880559/100‑1 (MQ=255)
                                               cGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAg                                                >  2:599758/1‑100 (MQ=255)
                                               cGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAg                                                <  2:2375693/100‑1 (MQ=255)
                                               cGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAg                                                <  2:1817508/100‑1 (MQ=255)
                                                gCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGCTACGCGACACTGAGCATGGTGCAGCTGGAAgg                                               >  2:67726/1‑100 (MQ=255)
                                                gCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAgg                                               <  2:2390467/100‑1 (MQ=255)
                                                gCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGATGGAAgg                                               >  2:1363341/1‑100 (MQ=255)
                                                 ccATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGc                                              >  2:1819221/1‑100 (MQ=255)
                                                  cATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGCTGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCt                                             >  2:1025735/1‑100 (MQ=255)
                                                  cATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCt                                             <  1:738045/100‑1 (MQ=255)
                                                    tACTTCTACGAGCACATCGGTAAGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACAGATCATAGTGGAGCTAGAAgacaac                                           >  2:2135657/1‑95 (MQ=255)
                                                    tACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTAc                                           <  2:2348348/100‑1 (MQ=255)
                                                    tACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTAc                                           >  2:933618/1‑100 (MQ=255)
                                                      cTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGt                                         >  1:1717039/1‑100 (MQ=255)
                                                      cTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGt                                         <  2:1601935/100‑1 (MQ=255)
                                                      cTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGt                                         >  1:331247/1‑100 (MQ=255)
                                                      cTTCGACGAACACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGt                                         >  1:2580101/1‑100 (MQ=255)
                                                       ttCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTg                                        <  2:26013/100‑1 (MQ=255)
                                                          gacgacCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTgaag                                     <  1:750844/100‑1 (MQ=255)
                                                           acgacCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTgaaga                                    <  2:2102485/100‑1 (MQ=255)
                                                            cgacCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTgaagaa                                   >  2:1253956/1‑100 (MQ=255)
                                                             gacCACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAAt                                  >  1:1448658/1‑100 (MQ=255)
                                                               ccACATCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATaa                                >  1:1679605/1‑100 (MQ=255)
                                                                  caccGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGAt                             >  2:1428328/4‑100 (MQ=255)
                                                                  caTCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGAt                             >  2:2490539/1‑100 (MQ=255)
                                                                   aTCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTCCGTGAAGAATAAGATc                            >  1:2172547/1‑100 (MQ=255)
                                                                   aTCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATc                            <  1:2333225/100‑1 (MQ=255)
                                                                   aTCGGTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATc                            <  2:623894/100‑1 (MQ=255)
                                                                      ggTACGCGCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAAt                         >  2:733979/1‑100 (MQ=255)
                                                                       gTACGCGCGTTCAATAACCAGATGGTCTGGCTGCGCGACAGGAGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGGGAAGAATAAGATCaacc                        <  2:1440281/100‑3 (MQ=255)
                                                                          cgcgCTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTcc                     >  1:1465998/1‑100 (MQ=255)
                                                                           gcgcTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATGTccc                    >  2:1645476/1‑100 (MQ=255)
                                                                           gcgcTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTccc                    <  2:2203268/100‑1 (MQ=255)
                                                                           gcgcTTTCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTccc                    >  2:1645559/1‑100 (MQ=255)
                                                                             gcTTTCAATAGCCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTg                  <  2:1610265/100‑1 (MQ=255)
                                                                               tttCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGACGGCTACGTGAAGAATAAGATCAATCTCCCTgcg                >  2:1163199/1‑100 (MQ=255)
                                                                               tttCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTgcg                <  1:2454254/100‑1 (MQ=255)
                                                                               tttCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTgcg                >  2:1617171/1‑100 (MQ=255)
                                                                                ttCAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTgcgc               >  2:835434/1‑100 (MQ=255)
                                                                                 tCAATAACCAGATGGTCTGGCTTCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTTcgcg              >  2:435020/1‑100 (MQ=255)
                                                                                 tCAATAACCAGAGGGTCTGGCTGCGTGACAGGGGATACGCTACACTGAGCATGGTGCACTTGGAAGGCTACGTGTAGAATAAGATCAATCACCCTGcgcg              <  2:702495/100‑1 (MQ=255)
                                                                                  cAATAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGa             >  2:719081/1‑100 (MQ=255)
                                                                                    aTAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGc           >  2:1738378/1‑100 (MQ=255)
                                                                                     tAACCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCg          >  1:1669079/1‑100 (MQ=255)
                                                                                       aCCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCggt        >  1:1079584/1‑100 (MQ=255)
                                                                                       aCCAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCggt        >  1:1731110/1‑100 (MQ=255)
                                                                                        ccAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCggtg       >  2:788611/1‑100 (MQ=255)
                                                                                         cAGATGGTCTGGCTGCGTGACGGGGGATACGCGACACTGAGCATGGTGCAGCTGCGAGTCTACGTGAAGAATAAGATCAATCTCACTGCGCGAGtggtgg      >  2:2472704/1‑100 (MQ=255)
                                                                                         cAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCggtgg      <  1:1395365/100‑1 (MQ=255)
                                                                                         cAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCggtgg      <  2:1929044/100‑1 (MQ=255)
                                                                                         cAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCggtgg      >  1:1672564/1‑100 (MQ=255)
                                                                                         cAGATGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCggtgg      >  1:441375/1‑100 (MQ=255)
                                                                                            aTGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGGTCAATCTCCCTGCGCGAGCGGTGGTGa   <  2:1323526/100‑1 (MQ=255)
                                                                                             tGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGAt  <  1:889200/100‑1 (MQ=255)
                                                                                             tGGTCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGAt  >  2:462249/1‑100 (MQ=255)
                                                                                                |                                                                                                
TGCTAACCTATCACCATATTCTGCGCGACGAAGAAAACACCCGTTTTCGCCATACTTCGACGACCACATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGAT  >  NC_000913/145595‑145787

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