Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 903,232 Δ1 bp coding (459/516 nt) ybjP ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913903,2320T.99.4% 1384.2 / ‑0.7 308coding (459/516 nt)ybjPhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (114/192);  minor base C (2/0);  total (116/192)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 319 total aligned reads displayed.

AATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGA  >  NC_000913/903144‑903319
                                                                                        |                                                                                       
aaTTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTa                                                                                    >  2:9443090/1‑93 (MQ=255)
aaTTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTa                                                                                    <  2:839276/93‑1 (MQ=255)
aaTTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTa                                                                                    >  1:5579979/1‑93 (MQ=255)
 aTTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCTATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTaa                                                                                   <  2:4345718/93‑1 (MQ=255)
 aTTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTaa                                                                                   <  2:3806895/93‑1 (MQ=255)
 aTTTACGGGGTTTATTGGTTGATCAAGCCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTaa                                                                                   <  1:5187834/93‑1 (MQ=255)
       ggggTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGcaca                                                                             <  1:5919213/93‑1 (MQ=255)
         ggTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACAtc                                                                           <  2:3234662/93‑1 (MQ=255)
         ggTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACAtc                                                                           <  1:8268679/93‑1 (MQ=255)
           tttATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACAtcgt                                                                         >  1:521140/1‑93 (MQ=255)
           tttATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACAtcgt                                                                         >  1:4556769/1‑93 (MQ=255)
           tttATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACAtcgt                                                                         <  1:6816934/93‑1 (MQ=255)
           tttATTGGTTAATCAAGGCTTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACAtcgt                                                                         <  2:1462062/93‑1 (MQ=255)
            ttATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGAGTGACGGAGCTTACCCGCCGTGGCATGAACGCTACCGCCATGTAACGCACAtcgtc                                                                        <  2:9920415/93‑1 (MQ=255)
             tATTGTTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACAtcgtcg                                                                       <  2:219416/93‑1 (MQ=255)
             tATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACAtcgtcg                                                                       >  1:8471348/1‑93 (MQ=255)
             tATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACAtcgtcg                                                                       >  2:5314910/1‑93 (MQ=255)
             tATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACAtcgtcg                                                                       >  2:6541719/1‑93 (MQ=255)
              tttGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGa                                                                      <  2:2325332/92‑1 (MQ=255)
              tttGGTTGAGCATGGCTTGAGCGTTTCTCGACGGACTGACGCGCCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGa                                                                      <  1:5314912/92‑1 (MQ=255)
              aTTGGTTGCTCAAGGGGTGAGCGATTCTCGATGGCCCGAGGGGGCGTACCCGCCGTGGCCTGAAGGCTACCGCCAAGTAACGCACATCGTCGa                                                                      <  2:8595125/93‑1 (MQ=255)
              aTTGGTTGATCAAGGCTTTAGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGa                                                                      <  1:118794/93‑1 (MQ=255)
              aTTGGTTGATCAAGGCGTTGGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGa                                                                      <  2:419374/93‑1 (MQ=255)
              aTTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGa                                                                      <  2:6645650/93‑1 (MQ=255)
              aTTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGa                                                                      >  1:6121559/1‑93 (MQ=255)
               ttGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGAt                                                                     <  2:6261862/93‑1 (MQ=255)
                  gTTGATCAGGGCGTTAGCGATTCTCGATGGACTGAGGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATggc                                                                  <  2:5445005/93‑3 (MQ=255)
                    tGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGccc                                                                >  1:1499347/1‑93 (MQ=255)
                     gATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCa                                                               <  2:4834334/93‑1 (MQ=255)
                      aTCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGAGGCATGACCGCAACCGCAAAGCAACTCACATAGTCTATGCCCCaa                                                              >  1:8673274/1‑93 (MQ=255)
                        cAAGGCTTTAGCGATTCACGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAaca                                                            <  1:8958314/93‑1 (MQ=255)
                        cAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAaca                                                            <  2:7900689/93‑1 (MQ=255)
                        cAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAaca                                                            <  2:5116536/93‑1 (MQ=255)
                         aaGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAacac                                                           >  2:9497687/1‑93 (MQ=255)
                          aGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACt                                                          <  2:4543731/93‑1 (MQ=255)
                          aGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACt                                                          >  2:3956833/1‑93 (MQ=255)
                          aGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACt                                                          <  1:8078231/93‑1 (MQ=255)
                           ggCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTg                                                         <  1:9138889/93‑1 (MQ=255)
                             cGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACCCTGGc                                                       >  1:9914729/1‑93 (MQ=255)
                              gTTAGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGcc                                                      <  1:5743789/93‑1 (MQ=255)
                               ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGGCAAGCAACGCACATCGTCGATGGCCCAACACTGGCCt                                                     >  1:226102/1‑93 (MQ=255)
                               tgaGCGATGCTGGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCt                                                     <  2:6063947/91‑1 (MQ=255)
                                    gATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTg                                                >  1:2502144/1‑93 (MQ=255)
                                     aTTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGa                                               <  2:7487761/93‑1 (MQ=255)
                                     aTTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGa                                               >  1:2081857/1‑93 (MQ=255)
                                      ttCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGaa                                              >  1:188957/1‑93 (MQ=255)
                                      ttCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGaa                                              >  2:2611822/1‑93 (MQ=255)
                                       tctcGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAAt                                             <  1:6011419/93‑1 (MQ=255)
                                         tcGATGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCa                                           <  1:8722299/93‑1 (MQ=255)
                                         tcGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCa                                           >  1:4279738/1‑93 (MQ=255)
                                         ccGATGGCCTGACGGGGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCa                                           <  1:7393426/92‑1 (MQ=255)
                                           gATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCAtt                                         <  1:8909482/93‑1 (MQ=255)
                                           gATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCAtt                                         <  2:901414/93‑1 (MQ=255)
                                           gATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCAtt                                         >  2:7339000/1‑93 (MQ=255)
                                            aTGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTa                                        >  1:509011/1‑93 (MQ=255)
                                            aTGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTa                                        >  2:3882130/1‑93 (MQ=255)
                                             tgtaCTGACAGAGCTTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTaa                                       <  2:5952976/90‑1 (MQ=255)
                                              ggACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAac                                      >  2:9222787/1‑93 (MQ=255)
                                               gacggacgGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAaca                                     <  2:1042531/93‑1 (MQ=255)
                                               gACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAaca                                     >  1:1854709/1‑93 (MQ=255)
                                                aCTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTGacac                                    >  2:7733276/1‑93 (MQ=255)
                                                aCTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAacac                                    <  1:3436262/93‑1 (MQ=255)
                                                  tGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACtt                                  >  2:1742016/1‑93 (MQ=255)
                                                   gACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTc                                 <  1:2802422/93‑1 (MQ=255)
                                                    ccGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCa                                <  2:975901/92‑1 (MQ=255)
                                                    aCGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCa                                <  1:9180007/93‑1 (MQ=255)
                                                      ggAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATc                              <  1:8026125/93‑1 (MQ=255)
                                                       gAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATcc                             >  2:1587610/1‑93 (MQ=255)
                                                       gAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATcc                             <  1:6422731/93‑1 (MQ=255)
                                                        aGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCt                            <  2:501208/93‑1 (MQ=255)
                                                        aGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCt                            <  1:669358/93‑1 (MQ=255)
                                                         gCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTg                           <  1:3234239/93‑1 (MQ=255)
                                                          cGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGc                          <  1:4073106/93‑1 (MQ=255)
                                                           gtccccGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGcc                         <  1:3685625/90‑1 (MQ=255)
                                                            tACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCa                        >  2:3684300/1‑93 (MQ=255)
                                                            tACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCa                        >  1:3834122/1‑93 (MQ=255)
                                                             aCCCGCCGTGGCATGAACGCTACCTCCAAGTAACGCACATCGTCGATGGCCTAACACTGGCCTTCCTGAATCATTAACACTTCATACTGCCaa                       >  1:7725308/1‑93 (MQ=255)
                                                              cccGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAc                      <  1:45420/93‑1 (MQ=255)
                                                              cccGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAc                      <  2:5252785/93‑1 (MQ=255)
                                                               ccgccgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAcc                     >  1:4939068/1‑93 (MQ=255)
                                                                cgccgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAccc                    >  1:967355/1‑93 (MQ=255)
                                                                 gccgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCGACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCt                   <  1:8560968/93‑1 (MQ=255)
                                                                 gccgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCt                   <  2:2072111/93‑1 (MQ=255)
                                                                 gccgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCt                   <  2:283423/93‑1 (MQ=255)
                                                                  ccgTGGCATAAACGCTACCGCCAAGTAACGCACATCGCCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTg                  <  1:9026833/93‑1 (MQ=255)
                                                                   cgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTgg                 >  2:1174055/1‑93 (MQ=255)
                                                                     tGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTc               <  2:8811384/93‑1 (MQ=255)
                                                                     tGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTc               <  1:1682857/93‑1 (MQ=255)
                                                                      ggCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCa              <  1:7552705/93‑1 (MQ=255)
                                                                          tGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCt          >  1:2834558/1‑93 (MQ=255)
                                                                            aaCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGt        <  2:2481465/93‑1 (MQ=255)
                                                                            aaCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGt        <  2:7392535/93‑1 (MQ=255)
                                                                            aaCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGt        <  1:6269787/93‑1 (MQ=255)
                                                                             aCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGtt       <  2:698362/93‑1 (MQ=255)
                                                                             aCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGtt       <  1:7042138/93‑1 (MQ=255)
                                                                             aCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGtt       >  1:1033389/1‑93 (MQ=255)
                                                                               gCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTg     <  2:8179706/93‑1 (MQ=255)
                                                                               gCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTg     >  1:4385120/1‑93 (MQ=255)
                                                                                cTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTga    >  2:6835105/1‑93 (MQ=255)
                                                                                cTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTga    >  1:4786832/1‑93 (MQ=255)
                                                                                        |                                                                                       
AATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGA  >  NC_000913/903144‑903319

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