Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,518,507 G→A 78.7% R14H (CGC→CAC)  yfeC → DUF1323 family putative DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,518,5070GA78.7% 317.5 / 80.2 164R14H (CGC→CAC) yfeCDUF1323 family putative DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base A (68/60);  total (87/77)
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

Only 100 of 176 total aligned reads displayed.

TAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAATGAA  >  NC_000913/2518408‑2518604
                                                                                                   |                                                                                                 
taAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAgccg                                                                                                   >  2:2004664/1‑100 (MQ=255)
 aaaaTGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAgccgc                                                                                                  >  1:1937806/1‑100 (MQ=255)
  aaaTGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAgccacc                                                                                                 <  2:90835/100‑1 (MQ=255)
  aaaTAAACTACTAATAGACCCACATACATTAGGGGAATTGTTATGTTCAAGGCGAGTATGACACCAGATGAACTTGCCAGACTGACCGGCTATAgccacc                                                                                                 <  2:945115/100‑1 (MQ=255)
    aTGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAg                                                                                               <  1:1008962/100‑1 (MQ=255)
    aTGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAg                                                                                               <  1:2324797/100‑1 (MQ=255)
      gAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGAc                                                                                             <  2:2254768/100‑1 (MQ=255)
      gAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGAc                                                                                             <  1:2668597/100‑1 (MQ=255)
       aaCTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGAcc                                                                                            <  1:3897352/100‑1 (MQ=255)
         ctactaATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCAt                                                                                          >  1:1751039/1‑100 (MQ=255)
           actaATAGACCCACATGCATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATta                                                                                        >  1:3240178/1‑100 (MQ=255)
           actaATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATta                                                                                        >  1:1838320/1‑100 (MQ=255)
              aaTAGACGCACATACATTCAGGGAATTGTTATGTGCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATtaata                                                                                     <  2:2878806/100‑1 (MQ=255)
                 aGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAAt                                                                                  >  1:2218491/1‑100 (MQ=255)
                  gACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATg                                                                                 <  1:3142405/100‑1 (MQ=255)
                  gACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATg                                                                                 >  2:3746017/1‑100 (MQ=255)
                  gACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATg                                                                                 <  1:765623/100‑1 (MQ=255)
                    cccACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTTTAGCCGCCAGACCATTAATAAATggg                                                                               <  1:2539498/100‑1 (MQ=255)
                        cgtacatTCAGGGAATTGTTATGTGCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTTATTAATggtctgg                                                                           <  1:3419117/98‑6 (MQ=255)
                        catacatTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACg                                                                           <  1:3513254/100‑1 (MQ=255)
                         atacatTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACg                                                                           <  2:3393604/99‑1 (MQ=255)
                          tacatTCAGGGAATTGTTATGTTCCAGGAGAGGATGACAACAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACgta                                                                         >  2:2383914/1‑98 (MQ=255)
                          tacatTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCa                                                                         >  2:2069044/1‑100 (MQ=255)
                          tacatTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCa                                                                         >  1:3764935/1‑100 (MQ=255)
                            catTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCaag                                                                       >  1:2383444/1‑100 (MQ=255)
                            catTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCaag                                                                       >  1:2735644/1‑100 (MQ=255)
                            catTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCaag                                                                       <  1:393700/100‑1 (MQ=255)
                              ttCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCaagga                                                                     <  2:3635970/100‑1 (MQ=255)
                               tCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCaaggaa                                                                    <  2:3034621/100‑1 (MQ=255)
                                cAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCaaggaag                                                                   <  1:3841624/100‑1 (MQ=255)
                                cAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCaaggaag                                                                   <  1:1200481/100‑1 (MQ=255)
                                cAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCaaggaag                                                                   >  2:309833/1‑100 (MQ=255)
                                 aGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCaaggaagg                                                                  >  2:2467050/1‑100 (MQ=255)
                                   ggAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCt                                                                >  1:1560336/1‑100 (MQ=255)
                                     aaTTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTgg                                                              <  2:2403038/100‑1 (MQ=255)
                                      aTTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGa                                                             <  1:3446068/100‑1 (MQ=255)
                                      aTTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGa                                                             <  1:1879239/100‑1 (MQ=255)
                                       ttGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGAc                                                            <  1:51818/100‑1 (MQ=255)
                                       ttGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGAc                                                            <  1:938171/100‑1 (MQ=255)
                                        tGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGa                                                             <  2:3361792/98‑1 (MQ=255)
                                         gTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACaa                                                          >  1:1409783/1‑100 (MQ=255)
                                         gTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACaa                                                          >  2:3138743/1‑100 (MQ=255)
                                         gTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACaa                                                          <  1:1261264/100‑1 (MQ=255)
                                         gTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACaa                                                          >  1:2121498/1‑100 (MQ=255)
                                         gTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTAAGCAAGGAAGGCTGGACaa                                                          >  2:1842174/1‑100 (MQ=255)
                                         gTTATGTTCAAGGAGAGAATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACaa                                                          >  2:1406656/1‑100 (MQ=255)
                                          ttATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAAc                                                         <  2:3806760/100‑1 (MQ=255)
                                          ttATGTTCAACGAGAGGATGACACCAGATGAACGTGCCAGACAGACCGGCTATAGCCACCAGACCACTATTAAATGGGTACGCAAGGAAGGCTGGACAAc                                                         <  2:2746083/100‑1 (MQ=255)
                                           tATGTTCAAGGAGAGGTTGACACCAGATGAACTTGCCAGACTGACAGGATATAGCCGCCAGACCATAAATAAATGGGTAGGCAAGGAAGGGTGGCCAAc                                                         >  2:553770/1‑99 (MQ=255)
                                           tATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACg                                                        <  2:3270862/100‑1 (MQ=255)
                                             tGTTCAAGGAGAAGACGACCCCCAATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAACTGGGTACGCAAGGAAGGCTGGACAACGTc                                                      <  2:4159768/100‑1 (MQ=255)
                                              gTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGGCCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCa                                                     >  2:1781743/1‑100 (MQ=255)
                                                tCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCAcc                                                   >  2:585424/1‑100 (MQ=255)
                                                tCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCAcc                                                   >  2:2160043/1‑100 (MQ=255)
                                                tCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCAcc                                                   >  1:1965407/1‑100 (MQ=255)
                                                tCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCAcc                                                   <  2:4033235/100‑1 (MQ=255)
                                                  aaGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCaa                                                 >  2:1109951/1‑100 (MQ=255)
                                                   nggAGAGGATGACACCAGATGAACTTGCCAGACTGAGCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCaca                                                >  2:3196871/2‑98 (MQ=255)
                                                        agGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTg                                           <  1:101704/100‑1 (MQ=255)
                                                         ggATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTgg                                          >  1:137508/1‑100 (MQ=255)
                                                           aTGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCg                                        <  1:1617613/100‑1 (MQ=255)
                                                            tGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGt                                       >  1:1439394/1‑100 (MQ=255)
                                                             gatacCAGATGAACTTGCAAGAAAGACCGGCTCTAGCCCACAGACCAATAAGAAATGAGTACGAAAGGAAGGCTGGCCAACAGCACCAacaca                                             >  2:3512870/4‑92 (MQ=255)
                                                             gACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACAAAAACCTGGCGt                                       >  2:1027503/1‑99 (MQ=255)
                                                              acacCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTcc                                     >  1:1077857/1‑100 (MQ=255)
                                                               cacCAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCa                                    >  2:3549919/1‑100 (MQ=255)
                                                                 ccAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAACGAAGGCTGGACAACTTCACCAAAACCTGGCGTCCAg                                   <  2:2970909/99‑1 (MQ=255)
                                                                  cAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAggg                                 <  1:930069/100‑1 (MQ=255)
                                                                  cAGATGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAggg                                 >  2:3999690/1‑100 (MQ=255)
                                                                     aTGAACTTGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTgg                              <  2:1531260/100‑1 (MQ=255)
                                                                       gAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCa                            <  2:1168894/100‑1 (MQ=255)
                                                                            tGCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAGAGcc                       >  1:3372286/1‑100 (MQ=255)
                                                                             gCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGccc                      >  1:2583979/1‑100 (MQ=255)
                                                                             gCCAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGccc                      >  1:2057761/1‑100 (MQ=255)
                                                                              ccAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCg                     >  2:4007836/1‑100 (MQ=255)
                                                                              ccAGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCg                     >  1:543120/1‑100 (MQ=255)
                                                                                aGACTGAGCGGCTATAGCCACCAGACCATTAATAAATGGGTCCGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGAAAAGCCCGAc                   >  2:3597670/1‑100 (MQ=255)
                                                                                aGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGAc                   <  1:1612986/100‑1 (MQ=255)
                                                                                aGACTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGAc                   <  1:2155322/100‑1 (MQ=255)
                                                                                  aCTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTg                 <  2:40627/100‑1 (MQ=255)
                                                                                   cTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACt                  >  2:3619634/1‑98 (MQ=255)
                                                                                   cTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTgg                <  2:3240157/100‑1 (MQ=255)
                                                                                   cTGACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTgg                <  1:4045705/100‑1 (MQ=255)
                                                                                     gACCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGtt              >  1:2876081/1‑100 (MQ=255)
                                                                                      aCCGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTc             <  1:4125043/100‑1 (MQ=255)
                                                                                       ccGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCa            <  1:3602332/100‑1 (MQ=255)
                                                                                        cGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCAc           <  1:1223787/100‑1 (MQ=255)
                                                                                        cGGCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCAGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCAc           >  2:2303068/1‑100 (MQ=255)
                                                                                         ggCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACg          >  1:705865/1‑100 (MQ=255)
                                                                                         ggCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACg          <  2:1427916/100‑1 (MQ=255)
                                                                                          gCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGATAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACg          <  1:1499554/99‑1 (MQ=255)
                                                                                          gCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGt         <  2:507261/100‑1 (MQ=255)
                                                                                          gCTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTAACg          >  1:637634/1‑99 (MQ=255)
                                                                                           cTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTc        <  1:589620/100‑1 (MQ=255)
                                                                                           cTATAGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTc        >  1:800510/1‑100 (MQ=255)
                                                                                             ataGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCaa      <  1:1537084/100‑1 (MQ=255)
                                                                                              taGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAAt     >  2:1046498/1‑100 (MQ=255)
                                                                                               aGCCACCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAATg    <  2:1464628/100‑1 (MQ=255)
                                                                                                gccaccAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAATGa   >  1:361215/1‑100 (MQ=255)
                                                                                                 ccgccAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAATGaa  <  1:3106963/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
TAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACCGGCTATAGCCGCCAGACCATTAATAAATGGGTACGCAAGGAAGGCTGGACAACGTCACCAAAACCTGGCGTCCAGGGTGGCAAAGCCCGACTGGTTCACGTCAATGAA  >  NC_000913/2518408‑2518604

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: