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.4% 716.2 / ‑2.3 161intergenic (+38/+13)yiiG/frvRhypothetical protein/putative frv operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (105/54);  total (106/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

GCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085587‑4085762
                                                                                        |                                                                                       
 cATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGa                                                                                   >  1:2679044/1‑93 (MQ=255)
  aTTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                  >  2:2762766/1‑93 (MQ=255)
  aTTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                  <  2:270964/93‑1 (MQ=255)
   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                 >  1:119164/1‑93 (MQ=255)
   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                 >  2:820176/1‑93 (MQ=255)
   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                 >  2:4469037/1‑93 (MQ=255)
   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                 <  2:3980339/93‑1 (MQ=255)
   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                 <  2:3331179/93‑1 (MQ=255)
   ttACGCGACTACAACGAAATGGTTGATGACTAGAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCAGCCGaca                                                                                 <  2:2975884/93‑1 (MQ=255)
    tACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGacac                                                                                <  2:1241954/93‑1 (MQ=255)
     aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGAcaca                                                                               >  1:852002/1‑93 (MQ=255)
      cgcgAGAACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAg                                                                              >  1:1547822/1‑93 (MQ=255)
       gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGt                                                                             <  2:2925885/93‑1 (MQ=255)
       gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGt                                                                             >  2:2360972/1‑93 (MQ=255)
        cgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGtt                                                                            <  2:2488110/93‑1 (MQ=255)
          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                          >  2:1076866/1‑93 (MQ=255)
          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCa                                                                          <  1:4327110/93‑1 (MQ=255)
            tACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                        <  1:2856062/93‑1 (MQ=255)
            tACAACGAAATGGTTGATGACTAAAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                        <  2:1924107/93‑1 (MQ=255)
              cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcccc                                                                      >  2:1799643/1‑93 (MQ=255)
               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccccc                                                                     >  2:4021375/1‑92 (MQ=255)
               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccccc                                                                     >  2:793937/1‑92 (MQ=255)
                aCGAAAGGGGTGATGACGAGAATACCCTGGATTACACGTTTTTAACGACAAGTGCCGTGTCGGATGCGGCGCATCCGACACAGTTCAccccct                                                                    <  2:3748980/93‑3 (MQ=255)
                 cGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtg                                                                   <  1:2116469/93‑1 (MQ=255)
                   aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCtgtg                                                                 >  1:3712563/1‑93 (MQ=255)
                     aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGaa                                                               <  1:669731/93‑1 (MQ=255)
                      tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAg                                                              <  1:2429035/93‑1 (MQ=255)
                      tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAg                                                              >  2:1475152/1‑93 (MQ=255)
                       ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGt                                                             >  1:2771917/1‑93 (MQ=255)
                          tgatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAAc                                                          <  1:3591515/93‑1 (MQ=255)
                           gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACt                                                         >  1:904822/1‑93 (MQ=255)
                            atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTc                                                        >  1:2793367/1‑93 (MQ=255)
                            atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTc                                                        >  1:2793366/1‑93 (MQ=255)
                             tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTcc                                                       >  1:2198266/1‑93 (MQ=255)
                             tgaCTATAATAGCCTGCATTAAACGTTTCTGATGATAAATGCAGTGTCGGCTGCGGCGCATCGGACACAGGTCACCCCCTGTGTAGTAACTcc                                                       >  2:3105084/1‑93 (MQ=255)
                              gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCa                                                      >  1:3744984/1‑93 (MQ=255)
                              gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCa                                                      >  1:234726/1‑93 (MQ=255)
                               aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                     >  2:356396/1‑93 (MQ=255)
                                cTATAATAGCCTGCGTTACACGTTTTTGATGATAAATGCAGTGTCGAATGCGGCGCATCCGACACAGTTAACCCCCTGTGAAGTAACTCCATa                                                    <  2:2916547/93‑1 (MQ=255)
                                cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                    <  2:1172617/93‑1 (MQ=255)
                                cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                    <  1:817616/93‑1 (MQ=255)
                                   taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                 <  1:3865504/93‑1 (MQ=255)
                                   taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                 >  1:3475180/1‑93 (MQ=255)
                                   taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                 >  1:4277722/1‑93 (MQ=255)
                                    aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATg                                                >  2:4472859/1‑93 (MQ=255)
                                      taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGCAGTAACTCCATAAATGcc                                              >  1:1663606/1‑93 (MQ=255)
                                      taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGcc                                              <  2:3512575/93‑1 (MQ=255)
                                      taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGcc                                              <  1:4490693/93‑1 (MQ=255)
                                       aGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCt                                             <  1:319084/93‑1 (MQ=255)
                                        gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTc                                            >  2:3110518/1‑93 (MQ=255)
                                        gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTc                                            >  1:1082908/1‑93 (MQ=255)
                                          cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGGGCATTCGACACAGCTCCACCCCTGTGAAGTAACTCCATAAATGCCTCat                                          >  2:4146175/1‑93 (MQ=255)
                                          cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCat                                          >  1:389532/1‑93 (MQ=255)
                                           tGCGTTAAACGTTTTTGATGATAAATGAAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         >  2:2741211/1‑93 (MQ=255)
                                            gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                        >  2:4554500/1‑93 (MQ=255)
                                            gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                        >  2:2093630/1‑93 (MQ=255)
                                            gCGTTAAACGTTGTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                        <  2:3191942/93‑1 (MQ=255)
                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                      >  1:2238727/1‑93 (MQ=255)
                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                      <  2:2677971/93‑1 (MQ=255)
                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                      >  2:2231915/1‑93 (MQ=255)
                                              gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATAGcc                                      >  2:732929/1‑93 (MQ=255)
                                               ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                     <  2:717588/93‑1 (MQ=255)
                                                tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                    >  1:3611643/1‑93 (MQ=255)
                                                  aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTgct                                  >  2:2890834/1‑91 (MQ=255)
                                                   aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                 >  1:4672260/1‑93 (MQ=255)
                                                   aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                 >  1:1419341/1‑93 (MQ=255)
                                                   aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                 >  2:3048673/1‑93 (MQ=255)
                                                    cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                <  1:3680915/93‑1 (MQ=255)
                                                     gTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                               >  1:3630530/1‑93 (MQ=255)
                                                      tttttGATGATAAATGCAGGGTCGGATGCGGCGCACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                              <  1:4158185/93‑1 (MQ=255)
                                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                             >  2:1343460/1‑93 (MQ=255)
                                                       ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                             <  1:4505959/93‑1 (MQ=255)
                                                         ttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                           >  1:1490925/1‑93 (MQ=255)
                                                         ttGATGATAAAGGCAGTGTCGGATGCGGCGAATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                           <  1:1305399/93‑1 (MQ=255)
                                                          tgatgaTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                          <  1:463390/93‑1 (MQ=255)
                                                           gatgatAAATGCAGTGTCGGATGCGGTGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCNTCATATCCTGTTAACCCGGCa                         >  1:635305/1‑93 (MQ=255)
                                                           gatgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATCAATGCCTCATATCCTGTTAACCCGGCa                         >  1:1327472/1‑93 (MQ=255)
                                                            atgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaa                        <  1:3771024/93‑1 (MQ=255)
                                                                tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                    >  2:4715245/1‑93 (MQ=255)
                                                                  aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtct                  >  1:2449191/1‑93 (MQ=255)
                                                                  aaTGCAGCGTCGGATGCGTCGGATCCGACACAGTTCACGCCCTGTGAAGTGACTCCATAAATGACTCGTATCCGGTTAACCCGGCAATAAtct                  >  2:3186238/1‑93 (MQ=255)
                                                                   aTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                 >  1:4373590/1‑93 (MQ=255)
                                                                    tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                >  1:3332017/1‑93 (MQ=255)
                                                                     gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg               >  1:2947764/1‑93 (MQ=255)
                                                                     gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg               >  1:2072757/1‑93 (MQ=255)
                                                                       aGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc             >  1:4328528/1‑93 (MQ=255)
                                                                       aGTGGCGGGGGCGGCGCATACGACCAAGATTACCCCCTGTGAAGGAACCCCAAAAATTCCTCATAACCCGTTCACCCGGCAATAATCTccccc             >  2:1876006/1‑93 (MQ=255)
                                                                        gtgtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgccg            >  1:3109135/1‑93 (MQ=255)
                                                                        gtgtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgccg            <  2:4306831/93‑1 (MQ=255)
                                                                          gtgggATGCGGCGCATCCGACGCAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa          <  1:2263315/90‑1 (MQ=255)
                                                                          gtCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa          >  2:1081663/1‑93 (MQ=255)
                                                                           tCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt         >  1:3256362/1‑93 (MQ=255)
                                                                           tCGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt         >  1:3883589/1‑93 (MQ=255)
                                                                            cGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg        >  2:3316711/1‑93 (MQ=255)
                                                                            cGGATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg        >  2:1589985/1‑93 (MQ=255)
                                                                             ggctGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc       <  1:2477385/90‑1 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGcc      <  1:202186/93‑2 (MQ=255)
                                                                              gATGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      >  1:4018969/1‑93 (MQ=255)
                                                                                tGCGGCGCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa    <  2:718240/93‑1 (MQ=255)
                                                                                  cggcgCATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca  >  1:1822515/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 < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: