Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 458,009 T→C intergenic (+85/‑103) clpX → / → lon ATP‑dependent specificity component of clpP serine protease/ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913458,0090TC87.0% 667.2 / 84.9 246intergenic (+85/‑103)clpX/lonATP‑dependent specificity component of clpP serine protease/ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (14/18);  new base C (112/102);  total (126/120)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 250 total aligned reads displayed.

GAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTA  >  NC_000913/457919‑458099
                                                                                          |                                                                                          
gAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCgg                                                                                          >  1:6210918/1‑93 (MQ=255)
gAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCgg                                                                                          <  1:3318702/93‑1 (MQ=255)
 aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCggg                                                                                         >  2:1110042/1‑93 (MQ=255)
 aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCggg                                                                                         <  1:1813696/93‑1 (MQ=255)
 aataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCggg                                                                                         <  2:6489674/93‑1 (MQ=255)
  ataatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCgggg                                                                                        >  1:2997889/1‑93 (MQ=255)
    aatTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGa                                                                                      <  1:3335085/93‑1 (MQ=255)
        aaCCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACa                                                                                  >  1:1246004/1‑93 (MQ=255)
         aCCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACAt                                                                                 <  2:4155528/93‑1 (MQ=255)
           cATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATcc                                                                               <  1:1613863/93‑1 (MQ=255)
           cATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATcc                                                                               >  2:4057/1‑93 (MQ=255)
            aTTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATccc                                                                              <  2:3177491/93‑1 (MQ=255)
             ttCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATcccc                                                                             >  1:2026694/1‑93 (MQ=255)
              tCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCa                                                                            <  2:1687970/93‑1 (MQ=255)
              tCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCa                                                                            <  2:2717873/93‑1 (MQ=255)
                ccATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCata                                                                          <  2:1048153/93‑1 (MQ=255)
                 cATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCatat                                                                         >  1:5287265/1‑93 (MQ=255)
                 cATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCatat                                                                         >  2:5669130/1‑93 (MQ=255)
                   tACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATAc                                                                       >  2:2136648/1‑93 (MQ=255)
                    aCAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACt                                                                      <  2:1189994/93‑1 (MQ=255)
                    aCAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACt                                                                      <  1:2860570/93‑1 (MQ=255)
                       aTTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGAc                                                                   >  2:3229702/1‑93 (MQ=255)
                       aTTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGAc                                                                   >  1:6420193/1‑93 (MQ=255)
                        ttAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACg                                                                  >  2:5845396/1‑93 (MQ=255)
                           gTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTAc                                                               <  2:3005035/93‑1 (MQ=255)
                            ttAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACa                                                              >  1:4866217/1‑93 (MQ=255)
                            ttAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACa                                                              <  1:6277035/93‑1 (MQ=255)
                              aaCCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATg                                                            <  1:6462430/93‑1 (MQ=255)
                              aaCCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATg                                                            >  2:5018338/1‑93 (MQ=255)
                                 cAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTa                                                         <  2:1789650/93‑1 (MQ=255)
                                  aaaaaGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTaa                                                        <  1:1685663/93‑1 (MQ=255)
                                   aaaaGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAAt                                                       >  2:448315/1‑93 (MQ=255)
                                    aaaGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATa                                                      >  1:1196933/1‑93 (MQ=255)
                                    aaaGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATa                                                      <  1:5539186/93‑1 (MQ=255)
                                       ggggggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGAt                                                   >  1:3465642/1‑93 (MQ=255)
                                        gggggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATg                                                  <  1:1349275/93‑1 (MQ=255)
                                          gggAGTTTATCCCCCCTTTAATTTTTCCTCGATTCTCGGCGTTGAATGCGGGGGAACCATCCCCATATACTGACGTACATGTTAATAGATGGc                                                <  1:3137176/93‑1 (MQ=255)
                                            gATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGt                                              >  2:2101161/1‑93 (MQ=255)
                                            gATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGt                                              >  1:4195906/1‑93 (MQ=255)
                                              ttttATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                                            >  2:2327061/1‑93 (MQ=255)
                                               tttATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGaa                                           <  2:469425/93‑1 (MQ=255)
                                                ttATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAg                                          <  1:4295891/93‑1 (MQ=255)
                                                ttATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAg                                          >  2:6000708/1‑93 (MQ=255)
                                                 tATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGc                                         <  1:6162552/93‑1 (MQ=255)
                                                  aTCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGca                                        <  2:3456437/93‑1 (MQ=255)
                                                    ctcCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGcaca                                      <  1:3946394/93‑1 (MQ=255)
                                                    ctcCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGcaca                                      <  1:3598354/93‑1 (MQ=255)
                                                     tcCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAg                                     >  1:3726084/1‑93 (MQ=255)
                                                     tcCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAg                                     >  1:4600770/1‑93 (MQ=255)
                                                     tcCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGCGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAg                                     <  1:4803717/93‑1 (MQ=255)
                                                      ccccTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGt                                    >  2:6315474/1‑93 (MQ=255)
                                                      ccccTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGt                                    <  1:5035601/93‑1 (MQ=255)
                                                       cccTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTc                                   <  1:1167728/93‑1 (MQ=255)
                                                        ccTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGTGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCg                                  >  2:2627385/1‑93 (MQ=255)
                                                        ccTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCg                                  <  2:4303671/93‑1 (MQ=255)
                                                        ccTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCg                                  >  1:6520287/1‑93 (MQ=255)
                                                         cTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgt                                 >  1:2557855/1‑93 (MQ=255)
                                                         cTTTAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgt                                 <  2:2861371/93‑1 (MQ=255)
                                                          tttAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtt                                >  2:160265/1‑92 (MQ=255)
                                                          tttAATTTTTCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtg                                >  2:542454/1‑93 (MQ=255)
                                                           ttAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATGCTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtgt                               >  1:3675861/1‑93 (MQ=255)
                                                              aTTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCAt                            >  1:5141838/1‑93 (MQ=255)
                                                               tttttCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATc                           >  1:3513507/1‑93 (MQ=255)
                                                               tttttCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATc                           <  1:4054134/93‑1 (MQ=255)
                                                                ttttCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCt                          <  2:1837479/93‑1 (MQ=255)
                                                                 tttCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTg                         <  2:2388551/93‑1 (MQ=255)
                                                                   tCCTCTATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGAt                       <  1:777581/93‑1 (MQ=255)
                                                                      tctATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAc                    >  1:5704738/1‑93 (MQ=255)
                                                                       ctATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTAcc                   <  2:818444/93‑1 (MQ=255)
                                                                        tATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCt                  >  1:3825478/1‑93 (MQ=255)
                                                                        tATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCt                  >  2:537212/1‑93 (MQ=255)
                                                                        tATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCt                  <  2:5032231/93‑1 (MQ=255)
                                                                        tATTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCt                  <  1:6068288/93‑1 (MQ=255)
                                                                         aTTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTg                 <  1:3075298/93‑1 (MQ=255)
                                                                         aTTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTg                 >  1:5902340/1‑93 (MQ=255)
                                                                         aTTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTg                 >  1:2816278/1‑93 (MQ=255)
                                                                         aTTCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTg                 >  2:3744557/1‑93 (MQ=255)
                                                                          ttCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTgg                >  1:4729746/1‑93 (MQ=255)
                                                                          ttCTCGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTgg                >  1:5212760/1‑93 (MQ=255)
                                                                           tctcGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGc               <  1:804008/93‑1 (MQ=255)
                                                                             tcGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCgg             >  1:5425086/1‑93 (MQ=255)
                                                                             tcGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCgg             >  2:5208374/1‑93 (MQ=255)
                                                                             tcGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCgg             >  1:2683888/1‑93 (MQ=255)
                                                                             tcGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCgg             >  1:811553/1‑93 (MQ=255)
                                                                              cGGCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGa            <  2:1490561/93‑1 (MQ=255)
                                                                               ggCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGaa           <  2:6580683/93‑1 (MQ=255)
                                                                                gCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGaaa          >  1:5324675/1‑93 (MQ=255)
                                                                                gCGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGaaa          <  2:3035263/93‑1 (MQ=255)
                                                                                 cGTTGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAAt         >  2:4793315/1‑93 (MQ=255)
                                                                                  gTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAAtt        >  2:5983156/1‑93 (MQ=255)
                                                                                   ttGAATGCGGGGGAAAGATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTa       <  2:36035/93‑1 (MQ=255)
                                                                                   ttGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTa       >  1:5221848/1‑93 (MQ=255)
                                                                                    tGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTaa      <  1:242592/93‑1 (MQ=255)
                                                                                    tGAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTaa      <  2:3432217/93‑1 (MQ=255)
                                                                                     gAATGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTaaa     <  1:2254529/93‑1 (MQ=255)
                                                                                      aaTGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAAc    >  2:5877255/1‑93 (MQ=255)
                                                                                      aaTGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAAc    <  2:2921530/93‑1 (MQ=255)
                                                                                       aTGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACt   <  2:5142897/93‑1 (MQ=255)
                                                                                       aTGCGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACt   >  2:5515557/1‑93 (MQ=255)
                                                                                        tgtgGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTa  <  1:2958758/93‑1 (MQ=255)
                                                                                          |                                                                                          
GAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTA  >  NC_000913/457919‑458099

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