Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,750,696 (C)5→4 coding (318/669 nt) ydhX ← putative 4Fe‑4S ferredoxin‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,750,6920C.99.6% 1485.8 / ‑2.6 281coding (322/669 nt)ydhXputative 4Fe‑4S ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (130/150);  total (130/151)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 310 total aligned reads displayed.

AGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTATCGGTGAC  >  NC_000913/1750592‑1750788
                                                                                                    |                                                                                                
 ggTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGcc                                                                                                 >  2:3208465/1‑101 (MQ=255)
  gTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGccc                                                                                                >  1:3208466/1‑101 (MQ=255)
   tAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGcccc                                                                                               >  1:3208470/1‑101 (MQ=255)
    aGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGccccg                                                                                              <  1:3207794/101‑2 (MQ=255)
     gCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGccccgt                                                                                             >  2:3208474/1‑99 (MQ=255)
      cGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGccccgtc                                                                                            <  1:3207984/101‑4 (MQ=255)
       gAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGccccgtcg                                                                                           >  1:3208476/1‑97 (MQ=255)
            tGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCaa                                                                                     <  1:3207856/101‑1 (MQ=255)
             ggTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCaaa                                                                                    >  2:3208482/1‑101 (MQ=255)
             ggTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCaaa                                                                                    <  1:3208107/101‑1 (MQ=255)
               tACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACg                                                                                  <  2:3207949/101‑1 (MQ=255)
               tACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACg                                                                                  >  1:3208487/1‑101 (MQ=255)
                aCGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGt                                                                                 >  1:3208490/1‑101 (MQ=255)
                 cGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCa                                                                                      >  2:3208492/1‑95 (MQ=255)
                 cGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTc                                                                                <  2:3207837/101‑1 (MQ=255)
                  ggACAAGCGCCAATGCAATAACTGCAGCCGCTGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCg                                                                               <  1:3208013/101‑1 (MQ=255)
                   gACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGa                                                                              <  2:3208021/101‑1 (MQ=255)
                    aCAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGAt                                                                             <  1:3207998/101‑1 (MQ=255)
                    aCAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGAt                                                                             >  1:3208494/1‑101 (MQ=255)
                     cAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATg                                                                            <  1:3208051/101‑1 (MQ=255)
                      aaGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGc                                                                           >  1:3208499/1‑101 (MQ=255)
                      aaGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGc                                                                           >  2:3208502/1‑101 (MQ=255)
                         cgcCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACg                                                                        >  1:3208506/1‑101 (MQ=255)
                         cgcCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACg                                                                        <  2:3207712/101‑1 (MQ=255)
                           ccAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGa                                                                      >  2:3208508/1‑101 (MQ=255)
                           ccAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGa                                                                      <  2:3208081/101‑1 (MQ=255)
                            cAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAg                                                                     >  2:3208511/1‑101 (MQ=255)
                             aaTGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGc                                                                    >  2:3208514/1‑101 (MQ=255)
                              aTGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCa                                                                   <  2:3208318/101‑1 (MQ=255)
                                gCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATc                                                                 >  2:3208518/1‑101 (MQ=255)
                                 cAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCt                                                                >  1:3208520/1‑101 (MQ=255)
                                  aaTAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCtt                                                               <  2:3207747/101‑1 (MQ=255)
                                    tAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCg                                                             <  1:3207911/101‑1 (MQ=255)
                                    tAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCg                                                             >  2:3208524/1‑101 (MQ=255)
                                     aaCTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGc                                                            <  2:3208249/101‑1 (MQ=255)
                                     aaCTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGc                                                            >  1:3208526/1‑101 (MQ=255)
                                      aCTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCa                                                           <  1:3207926/101‑1 (MQ=255)
                                       cTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAg                                                          >  2:3208531/1‑101 (MQ=255)
                                       cTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAg                                                          <  1:3207971/101‑1 (MQ=255)
                                        tGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGt                                                         >  2:3208535/1‑101 (MQ=255)
                                        tGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGt                                                         <  1:3207557/101‑1 (MQ=255)
                                           aGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCt                                                      >  1:3208537/1‑101 (MQ=255)
                                            gCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTg                                                     <  1:3208213/101‑1 (MQ=255)
                                               gATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACa                                                  >  2:3208542/1‑101 (MQ=255)
                                               gATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACa                                                  <  1:3207765/101‑1 (MQ=255)
                                               gATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACa                                                  <  2:3207819/101‑1 (MQ=255)
                                                aTGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACAt                                                 <  2:3208131/101‑1 (MQ=255)
                                                 tGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATg                                                >  1:3208546/1‑101 (MQ=255)
                                                     cTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTg                                            <  2:3208135/101‑1 (MQ=255)
                                                     cTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTg                                            >  2:3208549/1‑101 (MQ=255)
                                                      tGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTgg                                           <  2:3207979/101‑1 (MQ=255)
                                                       gCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGc                                          <  1:3207943/101‑1 (MQ=255)
                                                       gCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGc                                          <  1:3208432/101‑1 (MQ=255)
                                                       gCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGc                                          >  2:3208552/1‑101 (MQ=255)
                                                       gCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGc                                          <  2:3208168/101‑1 (MQ=255)
                                                        cGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGc                                                                           >  1:3208555/1‑67 (MQ=255)
                                                        cGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCg                                         <  1:3208007/101‑1 (MQ=255)
                                                          aTTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGaa                                       >  2:3208558/1‑101 (MQ=255)
                                                          aTTTTACTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGaa                                       >  1:3208559/1‑101 (MQ=255)
                                                              ttCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGaaaaaa                                   <  1:3208065/101‑1 (MQ=255)
                                                              ttCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGaaaaaa                                   <  2:3207854/101‑1 (MQ=255)
                                                              ttCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGaaaaaa                                   >  2:3208563/1‑101 (MQ=255)
                                                               tCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAg                                  <  2:3208109/101‑1 (MQ=255)
                                                                cTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGt                                 >  1:3208566/1‑101 (MQ=255)
                                                                 tACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTg                                <  1:3207894/101‑1 (MQ=255)
                                                                  aCCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGa                               >  1:3208570/1‑101 (MQ=255)
                                                                   cccGCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGAt                              >  2:3208571/1‑101 (MQ=255)
                                                                      gCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACt                           <  1:3208161/101‑1 (MQ=255)
                                                                      gCACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACt                           >  2:3208575/1‑101 (MQ=255)
                                                                       cACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTg                          <  2:3208046/101‑1 (MQ=255)
                                                                       cACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTg                          >  1:3208577/1‑101 (MQ=255)
                                                                       cACGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTg                          <  1:3208502/101‑1 (MQ=255)
                                                                        aCGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGc                         <  2:3207895/101‑1 (MQ=255)
                                                                        aCGATCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGc                         <  1:3208226/101‑1 (MQ=255)
                                                                           aTCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGtt                      <  2:3208147/101‑1 (MQ=255)
                                                                            tCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGttt                     >  2:3208584/1‑101 (MQ=255)
                                                                            tCCCCTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGttt                     >  1:3208582/1‑101 (MQ=255)
                                                                             ccccTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGttt                     <  2:3207923/100‑1 (MQ=255)
                                                                              cccTGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCg                   >  1:3208587/1‑101 (MQ=255)
                                                                                 tGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTg                <  2:3207963/101‑1 (MQ=255)
                                                                                 tGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTg                >  1:3208593/1‑101 (MQ=255)
                                                                                 tGCTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTg                >  2:3208590/1‑101 (MQ=255)
                                                                                   cTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTc              <  2:3207994/101‑1 (MQ=255)
                                                                                   cTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTc              <  1:3208295/101‑1 (MQ=255)
                                                                                   cTCATCACGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTc              <  2:3208271/101‑1 (MQ=255)
                                                                                     catcaCGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGt            >  1:3208596/1‑101 (MQ=255)
                                                                                      atcaCGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTCGTCGtt           <  1:3208077/101‑1 (MQ=255)
                                                                                       tcaCGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCg                   <  1:3208023/92‑1 (MQ=255)
                                                                                       tcaCGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTa          <  1:3208139/101‑1 (MQ=255)
                                                                                       tcaCGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTa          <  1:3208100/101‑1 (MQ=255)
                                                                                       tcaCGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTa          <  1:3208154/101‑1 (MQ=255)
                                                                                        caCGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTAt         >  2:3208600/1‑101 (MQ=255)
                                                                                        caCGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTAt         <  2:3208039/101‑1 (MQ=255)
                                                                                         aCGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCg                                         <  2:3207914/68‑1 (MQ=255)
                                                                                          cGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTATCg       >  1:3208602/1‑101 (MQ=255)
                                                                                          cGCCACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTATCg       <  1:3208087/101‑1 (MQ=255)
                                                                                            ccACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGtt           <  2:3208177/95‑1 (MQ=255)
                                                                                             cACGACGCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTATCGGTg    <  2:3208064/101‑1 (MQ=255)
                                                                                              acgacgCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTATCGGTGa   >  1:3208607/1‑101 (MQ=255)
                                                                                               cgacgCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTATCGGTGAc  <  1:3208122/101‑1 (MQ=255)
                                                                                                    |                                                                                                
AGGTAGCGAACCTGGTACGGACAAGCGCCAATGCAATAACTGCAGCCGATGCACTGCGATTTTTCTACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACTGGCGAAAAAAGTGATACTGCGTTTCGTTGTCGTTATCGGTGAC  >  NC_000913/1750592‑1750788

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