Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,085,675 Δ1 bp intergenic (+38/+13) yiiG → / ← frvR hypothetical protein/putative frv operon regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,085,6750G.99.5% 938.6 / ‑3.7 202intergenic (+38/+13)yiiG/frvRhypothetical protein/putative frv operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (117/84);  minor base C (0/1);  total (117/85)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 219 total aligned reads displayed.

GCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085587‑4085762
                                                                                        |                                                                                       
gCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCg                                                                                    >  1:4483411/1‑93 (MQ=255)
  aTTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                  >  2:5416394/1‑93 (MQ=255)
  aTTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                  <  1:597319/93‑1 (MQ=255)
     aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGAcaca                                                                               <  1:81990/93‑1 (MQ=255)
      cgcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAg                                                                              >  2:3375294/1‑93 (MQ=255)
      cgcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAg                                                                              <  2:3694004/93‑1 (MQ=255)
       gcgAGTACAACGAAATGGTTGATGACTATAATAGCTTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGt                                                                             <  1:1144181/93‑1 (MQ=255)
       gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGt                                                                             >  1:1799254/1‑93 (MQ=255)
         gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTc                                                                           <  1:4015053/93‑1 (MQ=255)
           gTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAc                                                                         >  2:4432378/1‑93 (MQ=255)
             aCAACGAAATGGTTGATGACTATAATAGCCTGCTTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                       <  2:2290049/93‑1 (MQ=255)
             aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                       >  2:31908/1‑93 (MQ=255)
             aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                       >  2:1843532/1‑93 (MQ=255)
             aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                       <  2:4706767/93‑1 (MQ=255)
              cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcccc                                                                      >  1:3308277/1‑93 (MQ=255)
              cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcccc                                                                      <  2:2699048/93‑1 (MQ=255)
               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCt                                                                     <  1:663866/93‑1 (MQ=255)
               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCt                                                                     <  2:4646027/93‑1 (MQ=255)
               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCt                                                                     >  1:4957835/1‑93 (MQ=255)
               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCt                                                                     >  1:4327747/1‑93 (MQ=255)
               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGCATGCGGCCCATCCGACACAGTTCCCCCCt                                                                     <  1:5376139/93‑1 (MQ=255)
                aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCtt                                                                    <  2:2799284/93‑1 (MQ=255)
                   aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgtg                                                                 >  2:2523002/1‑93 (MQ=255)
                   aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgtg                                                                 >  2:3137202/1‑93 (MQ=255)
                    aaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                >  2:3938250/1‑93 (MQ=255)
                    aaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                >  1:5162227/1‑93 (MQ=255)
                    aaTGGTTAATAACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                <  1:3808759/93‑1 (MQ=255)
                     aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                               >  1:3128894/1‑93 (MQ=255)
                           gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACt                                                         >  1:4827713/1‑93 (MQ=255)
                              gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCCACACAGTTTACCCCTTGTGAAGTAACTCCa                                                      >  1:2243286/1‑93 (MQ=255)
                               aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCAt                                                     >  1:4894951/1‑93 (MQ=255)
                                cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATa                                                    >  1:2636725/1‑93 (MQ=255)
                                 tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGAAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGCAACGCCATaa                                                   >  2:4055845/1‑93 (MQ=255)
                                  ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                  <  2:2618835/93‑1 (MQ=255)
                                   taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAAt                                                 >  1:5024799/1‑93 (MQ=255)
                                   taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAAt                                                 >  2:5551432/1‑93 (MQ=255)
                                    aataGCCTGCGTTAAACGTTTTTGGTGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATg                                                <  1:3235822/93‑1 (MQ=255)
                                    aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATg                                                <  2:6115159/93‑1 (MQ=255)
                                     ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGc                                               >  1:1268926/1‑93 (MQ=255)
                                     ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGc                                               >  2:5330076/1‑93 (MQ=255)
                                      taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                              <  2:5834643/93‑1 (MQ=255)
                                      taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                              >  1:5583101/1‑93 (MQ=255)
                                      taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                              <  2:1749259/93‑1 (MQ=255)
                                        gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                                            <  2:4254299/93‑1 (MQ=255)
                                        gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                                            >  2:4816677/1‑93 (MQ=255)
                                         ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCa                                           >  2:187089/1‑93 (MQ=255)
                                         ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGCGAAGTAACTCCATAAATGCCTCa                                           >  2:1358120/1‑93 (MQ=255)
                                          cTGCGTTAACCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCat                                          <  2:216697/93‑1 (MQ=255)
                                          cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCat                                          <  2:3113519/93‑1 (MQ=255)
                                           tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCata                                         <  2:3592719/93‑1 (MQ=255)
                                           tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCata                                         >  2:5998648/1‑93 (MQ=255)
                                             cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                                       <  2:4735805/93‑1 (MQ=255)
                                             cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                                       >  2:3025158/1‑93 (MQ=255)
                                              gTTAAACTTTTTTGATAATAAATGCAGTGTGGGATGCGGCGCACCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                      <  2:4385869/93‑1 (MQ=255)
                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCTGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATCTcc                                      >  1:3701605/1‑93 (MQ=255)
                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                      <  2:75033/93‑1 (MQ=255)
                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                    >  1:1327812/1‑93 (MQ=255)
                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                    >  2:2931191/1‑93 (MQ=255)
                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                    >  1:4445637/1‑93 (MQ=255)
                                                 aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                   >  2:3963053/1‑93 (MQ=255)
                                                    cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                <  1:2041851/93‑1 (MQ=255)
                                                     gTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                               <  2:4105134/93‑1 (MQ=255)
                                                     gTTTTTGATAAAAAAGGCAGTGTCGGATGCGGCCCACCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                               <  1:115115/93‑1 (MQ=255)
                                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                             >  2:1193666/1‑93 (MQ=255)
                                                         ttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                           <  1:1521722/93‑1 (MQ=255)
                                                           gatgatAAATGCAGTGTCGTATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                         >  1:6291356/1‑93 (MQ=255)
                                                           gatgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                         >  2:807003/1‑93 (MQ=255)
                                                           gatgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                         >  1:73260/1‑93 (MQ=255)
                                                            atgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaa                        <  2:660292/93‑1 (MQ=255)
                                                               atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                     >  1:3945851/1‑93 (MQ=255)
                                                               atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                     <  1:6383677/93‑1 (MQ=255)
                                                                tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                    <  1:2832668/93‑1 (MQ=255)
                                                                tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                    <  2:4316396/93‑1 (MQ=255)
                                                                aaaaTGGAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                    <  2:1345717/92‑1 (MQ=255)
                                                                 aaaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAACCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc                   >  2:2049031/1‑93 (MQ=255)
                                                                  aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtct                  <  1:754419/93‑1 (MQ=255)
                                                                  aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtct                  >  2:581523/1‑93 (MQ=255)
                                                                  aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtct                  <  2:5162291/93‑1 (MQ=255)
                                                                    tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                <  2:5609536/93‑1 (MQ=255)
                                                                    tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                <  2:4620925/93‑1 (MQ=255)
                                                                    tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                >  2:39098/1‑93 (MQ=255)
                                                                     gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg               >  1:1488560/1‑93 (MQ=255)
                                                                     gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAACCGGCAATAATCTccg               >  1:4771447/1‑93 (MQ=255)
                                                                     gCAGCGTCGAAGGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg               <  1:3394109/93‑1 (MQ=255)
                                                                      cgtggtCGGAGGCGGCCCACCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgc              <  2:4483411/89‑1 (MQ=255)
                                                                       aGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc             >  2:6441628/1‑93 (MQ=255)
                                                                       aGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc             <  2:6360870/93‑1 (MQ=255)
                                                                           tCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt         <  1:4649956/93‑1 (MQ=255)
                                                                           tCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt         <  1:3775518/93‑1 (MQ=255)
                                                                             ggATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc       >  2:58909/1‑93 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      >  1:3494855/1‑93 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      <  2:1067105/93‑1 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      >  2:3792826/1‑93 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      >  2:140480/1‑93 (MQ=255)
                                                                               aTGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg     <  1:5759769/93‑1 (MQ=255)
                                                                                 gcggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac   >  2:1727523/1‑93 (MQ=255)
                                                                                 gcggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac   <  1:5840435/93‑1 (MQ=255)
                                                                                 gcggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac   <  2:1701087/93‑1 (MQ=255)
                                                                                 gcggcgCACCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac   <  1:5422587/93‑1 (MQ=255)
                                                                                  cggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca  >  2:5047643/1‑93 (MQ=255)
                                                                                        |                                                                                       
GCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085587‑4085762

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: