Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,380,060 C→T G153D (GGT→GAT)  frdB ← fumarate reductase iron‑sulfur protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,380,0600CT99.7% 1381.0 / ‑4.4 326G153D (GGT→GAT) frdBfumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (158/167);  total (159/167)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 326 total aligned reads displayed.

CGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGAT  >  NC_000913/4379962‑4380157
                                                                                                  |                                                                                                 
 gcgGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATcac                                                                                                <  2:8135012/101‑1 (MQ=255)
 gcgGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATcac                                                                                                <  1:8135208/101‑1 (MQ=255)
   ggCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAg                                                                                              >  1:8135665/1‑101 (MQ=255)
   ggCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAg                                                                                              >  2:8135667/1‑101 (MQ=255)
    gCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGt                                                                                             <  2:8135195/101‑1 (MQ=255)
    gCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGt                                                                                             <  1:8135081/101‑1 (MQ=255)
     cTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGtt                                                                                            <  2:8135052/101‑1 (MQ=255)
     cTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGtt                                                                                            <  1:8135126/101‑1 (MQ=255)
       aTCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGa                                                                                          >  1:8135674/1‑101 (MQ=255)
        tCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATcaca                                                                                               >  2:8135675/1‑95 (MQ=255)
          ttCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGc                                                                                       >  1:8135680/1‑101 (MQ=255)
          ttCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGc                                                                                       <  2:8135020/101‑1 (MQ=255)
          ttCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGc                                                                                       <  2:8135371/101‑1 (MQ=255)
          ttCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGc                                                                                       >  2:8135681/1‑101 (MQ=255)
            cGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCaa                                                                                     >  2:8135687/1‑101 (MQ=255)
            cGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCaa                                                                                     <  2:8135405/101‑1 (MQ=255)
             gTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGc                                                                                       <  2:8134965/98‑1 (MQ=255)
              ttATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAAcc                                                                                   >  1:8135690/1‑101 (MQ=255)
               tataACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCg                                                                                  <  2:8135303/101‑1 (MQ=255)
                ataACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCgg                                                                                 <  1:8135255/101‑1 (MQ=255)
                 taACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGga                                                                                >  1:8135696/1‑101 (MQ=255)
                 taACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGga                                                                                >  2:8135695/1‑101 (MQ=255)
                   aCGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGgaga                                                                              <  2:8135120/101‑1 (MQ=255)
                   aCGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGgaga                                                                              >  1:8135700/1‑101 (MQ=255)
                    cGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGtt                                                                                            >  1:8135706/1‑86 (MQ=255)
                      aTGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACt                                                                           >  1:8135710/1‑101 (MQ=255)
                      aTGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACt                                                                           <  1:8135317/101‑1 (MQ=255)
                      aTGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACt                                                                           >  2:8135709/1‑101 (MQ=255)
                       tGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTg                                                                          <  2:8135186/101‑1 (MQ=255)
                        gcgcCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTgg                                                                         <  2:8135435/101‑1 (MQ=255)
                        gcgcCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTgg                                                                         >  2:8135715/1‑101 (MQ=255)
                         cgcCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGt                                                                        <  1:8135103/101‑1 (MQ=255)
                          gcCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTg                                                                       >  2:8135718/1‑101 (MQ=255)
                           ccAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGa                                                                      <  1:8134810/101‑1 (MQ=255)
                              gCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATAc                                                                   >  2:8135722/1‑101 (MQ=255)
                              gCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATAc                                                                   >  1:8135724/1‑101 (MQ=255)
                                gTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGAt                                                                                         >  1:8135728/1‑77 (MQ=255)
                                gTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACtt                                                                 <  2:8135078/101‑1 (MQ=255)
                                 tAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTc                                                                >  1:8135731/1‑101 (MQ=255)
                                    tGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGcc                                                             >  2:8135733/1‑101 (MQ=255)
                                     ggCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCa                                                            <  1:8134989/101‑1 (MQ=255)
                                      gCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCAt                                                           <  2:8134858/101‑1 (MQ=255)
                                      gCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCAt                                                           >  2:8135738/1‑101 (MQ=255)
                                       cAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATc                                                          >  1:8135739/1‑101 (MQ=255)
                                       cAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATc                                                          <  1:8135037/101‑1 (MQ=255)
                                        aGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCt                                                         <  1:8135114/101‑1 (MQ=255)
                                         gCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTg                                                        >  1:8135744/1‑101 (MQ=255)
                                          ccGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTgc                                                       <  1:8135589/101‑1 (MQ=255)
                                           cGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTgcg                                                      <  2:8135329/101‑1 (MQ=255)
                                           cGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTgcg                                                      >  1:8135749/1‑101 (MQ=255)
                                             gACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGcc                                                    >  2:8135751/1‑101 (MQ=255)
                                              aCCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCg                                                   <  1:8135226/101‑1 (MQ=255)
                                              aCCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCg                                                   >  1:8135755/1‑101 (MQ=255)
                                               ccGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCgg                                                  <  2:8135102/101‑1 (MQ=255)
                                                cGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCggg                                                 <  1:8134934/101‑1 (MQ=255)
                                                cGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCggg                                                 >  1:8135758/1‑101 (MQ=255)
                                                 gATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCgggg                                                >  2:8135760/1‑101 (MQ=255)
                                                  aTGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGt                                               <  2:8135285/101‑1 (MQ=255)
                                                   tGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTc                                              <  1:8135171/101‑1 (MQ=255)
                                                    gAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCt                                             >  1:8135766/1‑101 (MQ=255)
                                                      aCTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAGCCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTgg                                           >  2:8135768/1‑101 (MQ=255)
                                                      aCTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTgg                                           >  1:8135771/1‑101 (MQ=255)
                                                       ctctGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGa                                          >  2:8135773/1‑101 (MQ=255)
                                                       ctctGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGa                                          <  1:8134958/101‑1 (MQ=255)
                                                        tctGGGTTCAGGTCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGAt                                         <  1:8135119/101‑1 (MQ=255)
                                                         ctGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATg                                        <  2:8134739/101‑1 (MQ=255)
                                                         ctGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATg                                        <  1:8135065/101‑1 (MQ=255)
                                                            ggTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTa                                     <  2:8135267/101‑1 (MQ=255)
                                                            ggTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTa                                     <  1:8135291/101‑1 (MQ=255)
                                                                   gCCAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCt                              >  1:8135778/1‑101 (MQ=255)
                                                                     cAAACTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGa                            <  1:8135139/101‑1 (MQ=255)
                                                                      aaaCTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGAt                           <  2:8135227/101‑1 (MQ=255)
                                                                      aaaCTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGAt                           >  2:8135781/1‑101 (MQ=255)
                                                                      aaaCTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGAt                           >  1:8135782/1‑101 (MQ=255)
                                                                        aCTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATcc                         <  2:8135466/101‑1 (MQ=255)
                                                                        aCTGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATcc                         <  2:8135253/101‑1 (MQ=255)
                                                                          tGCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCcgc                       <  1:8135472/101‑1 (MQ=255)
                                                                           gCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATcc                         <  1:8135155/98‑1 (MQ=255)
                                                                           gCGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCcgcg                      >  1:8135788/1‑101 (MQ=255)
                                                                            cGGGCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCgg                     >  2:8135790/1‑101 (MQ=255)
                                                                             gggCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCt                              <  1:8135179/91‑1 (MQ=255)
                                                                             gggCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGt                    >  1:8135793/1‑101 (MQ=255)
                                                                               gCACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGc                  <  2:8135148/101‑1 (MQ=255)
                                                                                cACGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCg                 <  2:8135210/101‑1 (MQ=255)
                                                                                 aCGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCgg                <  2:8134693/101‑1 (MQ=255)
                                                                                 aCGCGGCGTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCgg                >  1:8135797/1‑101 (MQ=255)
                                                                                    cggcgTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAg             >  2:8135801/1‑101 (MQ=255)
                                                                                     ggcgTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGt            >  1:8135803/1‑101 (MQ=255)
                                                                                      gcgTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGtt           <  1:8134894/101‑1 (MQ=255)
                                                                                      gcgTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGtt           <  1:8135233/101‑1 (MQ=255)
                                                                                       cgTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTg          >  1:8135807/1‑101 (MQ=255)
                                                                                        gTAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGc         <  1:8135330/101‑1 (MQ=255)
                                                                                         tAGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGcc        >  2:8135811/1‑101 (MQ=255)
                                                                                          aGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCg       >  1:8135816/1‑101 (MQ=255)
                                                                                          aGCACAAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCg       >  1:8135813/1‑101 (MQ=255)
                                                                                            cacaAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCgat     <  2:8135616/101‑1 (MQ=255)
                                                                                            cacaAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCgat     <  1:8135411/101‑1 (MQ=255)
                                                                                            cacaAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCgat     >  2:8135820/1‑101 (MQ=255)
                                                                                             acaAATCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCgatg    <  1:8135341/101‑1 (MQ=255)
                                                                                               aaaTCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCgatgat  <  2:8135167/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
CGCGGCTATCTTCGTTATAACGATGCGCCAGCGTAATGGCAGCCGGACCGATGAACTCTGGGTTCAGGCCAAACTGCGGGCACGCGGCGTAGCACAAACCACAGTTGATGCAACCGGAGAACTGGTGATACTTCGCCATCTGCGCCGGGGTCTGGATGTTAGTACCCTGATCCGCGGTGCGGGAGTTGCCGATGAT  >  NC_000913/4379962‑4380157

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: ‑