Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,080,195 A→G C49R (TGT→CGT)  fdoH ← formate dehydrogenase‑O, iron‑sulfur subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,080,1950AG98.9% 545.0 / ‑3.5 174C49R (TGT→CGT) fdoHformate dehydrogenase‑O, iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (85/87);  total (86/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 186 total aligned reads displayed.

CGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGG  >  NC_000913/4080106‑4080284
                                                                                         |                                                                                         
cGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACgcgc                                                                                        >  1:18117/1‑93 (MQ=255)
 gTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGCACGCGcc                                                                                       >  2:2371025/1‑93 (MQ=255)
 gTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGcc                                                                                       >  2:4333301/1‑93 (MQ=255)
  tCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCa                                                                                      >  2:2491476/1‑93 (MQ=255)
   ccACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCAc                                                                                     <  2:688267/93‑1 (MQ=255)
   ccACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCAc                                                                                     >  1:4146852/1‑93 (MQ=255)
     aCGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCACTGTTATTGCCGACGGTATAGCGACTGTCGTTCCACTCTGAACGCGACACCt                                                                                   >  1:1069280/1‑93 (MQ=255)
     aCGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCt                                                                                   >  2:3938709/1‑93 (MQ=255)
        aTTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGAc                                                                                <  1:1068745/93‑1 (MQ=255)
        aTTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGAc                                                                                >  2:4651607/1‑93 (MQ=255)
         tttGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACa                                                                               <  1:168469/93‑1 (MQ=255)
         tttGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACa                                                                               <  2:1382229/93‑1 (MQ=255)
         tttGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACa                                                                               <  1:4112788/93‑1 (MQ=255)
          ttGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAg                                                                              <  1:2960813/93‑1 (MQ=255)
          ttGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAg                                                                              <  2:365740/93‑1 (MQ=255)
            ggCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGc                                                                            >  2:4304377/1‑93 (MQ=255)
             gcgcTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCt                                                                           <  2:142720/93‑1 (MQ=255)
              cgcTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCtt                                                                          <  2:2803298/93‑1 (MQ=255)
               gcTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCttt                                                                         <  1:2589647/93‑1 (MQ=255)
                  tAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACa                                                                      <  1:4402576/93‑1 (MQ=255)
                  tAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACa                                                                      <  1:2060167/93‑1 (MQ=255)
                   aaaTCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAg                                                                     <  1:3000067/93‑1 (MQ=255)
                       cATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGa                                                                 >  1:4002118/1‑93 (MQ=255)
                        aTTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGAt                                                                <  2:2813382/93‑1 (MQ=255)
                        aTTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGAt                                                                <  1:2397175/93‑1 (MQ=255)
                         ttGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATa                                                               >  1:4644301/1‑93 (MQ=255)
                          tGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATAc                                                              <  1:232445/93‑1 (MQ=255)
                          tGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATAc                                                              >  2:1071391/1‑93 (MQ=255)
                           ggggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACa                                                             <  1:2029021/93‑1 (MQ=255)
                            gggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAg                                                            <  1:1847674/93‑1 (MQ=255)
                             ggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAgg                                                           <  1:2078057/93‑1 (MQ=255)
                             ggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAgg                                                           <  2:2681315/93‑1 (MQ=255)
                              gTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAggt                                                          >  2:4731512/1‑93 (MQ=255)
                               ttGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAggtg                                                         <  2:343118/93‑1 (MQ=255)
                                tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAggtgg                                                        <  1:4290428/93‑1 (MQ=255)
                                 gTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAggtggt                                                       >  2:4145212/1‑93 (MQ=255)
                                  tCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTa                                                      >  1:4507485/1‑93 (MQ=255)
                                   cGTACACCCCAATGTCATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTaa                                                     >  2:858435/1‑93 (MQ=255)
                                    gTACACCCGAATTTTATTGCCGACGGTATCGCGAATGTCGGTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAAc                                                    <  2:1781594/93‑1 (MQ=21)
                                     tACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACg                                                   >  2:1707721/1‑93 (MQ=255)
                                     tACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACg                                                   >  2:2649247/1‑93 (MQ=255)
                                       cacCCCAATGTGATTGCCGACGGTATCGCGAATGTCGTACCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTc                                                 <  2:2180741/93‑1 (MQ=255)
                                        acCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGGAAGGCTTTAAAGCCGATACAGGTGGTCACGTCg                                                >  2:4090286/1‑93 (MQ=25)
                                        acCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCg                                                >  1:3486016/1‑93 (MQ=255)
                                         cgccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGACCGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGa                                               <  2:3298340/91‑1 (MQ=255)
                                         ccccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGa                                               <  2:30373/93‑1 (MQ=255)
                                          cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGAt                                              >  1:306274/1‑93 (MQ=255)
                                          cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGAt                                              <  1:2208874/93‑1 (MQ=255)
                                          cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGAt                                              <  1:1789634/93‑1 (MQ=255)
                                          cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGAt                                              <  2:1680547/93‑1 (MQ=255)
                                          cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGAt                                              <  2:3867869/93‑1 (MQ=255)
                                           ccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGACGATACAGGTGGTAACGTCGATg                                             >  2:2983270/1‑93 (MQ=255)
                                            cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGa                                            <  1:2291317/93‑1 (MQ=255)
                                            cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGa                                            <  1:3646569/93‑1 (MQ=255)
                                             aaTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAg                                           >  2:1176510/1‑93 (MQ=255)
                                              aTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGt                                          >  2:2147258/1‑93 (MQ=255)
                                              aTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGt                                          >  2:312574/1‑93 (MQ=255)
                                               tGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGtt                                         >  2:3874793/1‑93 (MQ=255)
                                               tGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGtt                                         >  2:4570266/1‑93 (MQ=255)
                                               tGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGtt                                         >  1:2240573/1‑93 (MQ=255)
                                                gTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGttt                                        >  2:1829827/1‑93 (MQ=255)
                                                 ttATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTc                                       <  2:4716461/93‑1 (MQ=255)
                                                 ttATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTc                                       >  2:2283863/1‑93 (MQ=255)
                                                 ttATTGCCGAAGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTc                                       >  2:3555093/1‑93 (MQ=255)
                                                    ttGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGcc                                    <  1:2918020/93‑1 (MQ=255)
                                                      tccGAGGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCAc                                  <  2:3357973/92‑1 (MQ=255)
                                                      gCCGACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCAc                                  <  1:2926156/93‑1 (MQ=255)
                                                         gACGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCActtc                               <  1:1398009/93‑1 (MQ=255)
                                                           cGGTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCActtctt                             >  2:1038078/1‑93 (MQ=255)
                                                            ggTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTTGATGAGTTTCGCCACttcttc                            <  1:378744/93‑1 (MQ=255)
                                                            ggTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACttcttc                            >  2:3250443/1‑93 (MQ=255)
                                                            ggTATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACttcttc                            <  1:4686143/93‑1 (MQ=255)
                                                              tATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCt                          >  2:2250371/1‑93 (MQ=255)
                                                              tATCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGAGGTAACGTCGATGAGTTTCGCCACTTCTTCCt                          >  2:1932931/1‑93 (MQ=18)
                                                              tATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCt                          <  1:3827400/93‑1 (MQ=255)
                                                               aTCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTg                         <  2:1016874/93‑1 (MQ=255)
                                                               aTCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTg                         <  1:3568447/93‑1 (MQ=255)
                                                                tCGCGAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTgg                        >  2:3624097/1‑93 (MQ=255)
                                                                 cgcgAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGa                       >  1:1559835/1‑93 (MQ=255)
                                                                  gcgAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGaa                      <  2:4049398/93‑1 (MQ=255)
                                                                   cgAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAg                     <  2:1329965/93‑1 (MQ=255)
                                                                    gAATGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGt                    >  1:3984593/1‑93 (MQ=255)
                                                                     aaTGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTc                   >  2:929986/1‑93 (MQ=255)
                                                                      aTGTCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTcc                  >  1:1405718/1‑93 (MQ=255)
                                                                         tCGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGCAGGTccgcg               >  2:2410416/1‑90 (MQ=18)
                                                                          cGTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCcgcc              >  1:898429/1‑92 (MQ=255)
                                                                           gTTCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCgcgc             <  1:4540983/93‑1 (MQ=255)
                                                                             tCCACTCTGAACGCGCCACCTGACAGGCTTTCCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCt           <  1:515924/93‑1 (MQ=255)
                                                                             tCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCt           <  1:3251791/93‑1 (MQ=255)
                                                                             tCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCt           <  1:2169517/93‑1 (MQ=255)
                                                                             tCCACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCt           >  2:252677/1‑93 (MQ=255)
                                                                              ccACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGcccc          >  1:947586/1‑91 (MQ=255)
                                                                               cACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGa         <  1:4257272/93‑1 (MQ=255)
                                                                               cACTCTGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGa         <  1:1911986/93‑1 (MQ=255)
                                                                                  tctGAACGCGCTACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc      <  2:4512944/93‑1 (MQ=255)
                                                                                  tctGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc      >  1:1583142/1‑93 (MQ=255)
                                                                                  tctGAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACCTCTTCCTGGAAGTACCGCGCCTGAggg      >  1:4428414/1‑92 (MQ=25)
                                                                                     gAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGcgcg   <  1:4192181/93‑1 (MQ=255)
                                                                                     gAACGCGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGcgcg   >  1:3226345/1‑93 (MQ=255)
                                                                                     gAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCAGCGCCTGAGGcgcg   >  2:3924238/1‑93 (MQ=255)
                                                                                         |                                                                                         
CGTCCACGATTTGGCGCTTAAATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGG  >  NC_000913/4080106‑4080284

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