Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,486,588 A→G D348D (GAT→GAC yjgR ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,486,5880AG82.4% 645.6 / 143.5 307D348D (GAT→GACyjgRhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (120/133);  total (147/160)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 307 total aligned reads displayed.

TCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAC  >  NC_000913/4486497‑4486677
                                                                                           |                                                                                          
   cACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCaa                                                                                        >  2:193226/1‑93 (MQ=255)
    aCGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCaaa                                                                                       >  1:4920352/1‑93 (MQ=255)
     cGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAAc                                                                                      >  1:1625454/1‑93 (MQ=255)
       cTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGc                                                                                    >  1:3623662/1‑93 (MQ=255)
       cTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGc                                                                                    <  1:5475482/93‑1 (MQ=255)
        tCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGcc                                                                                   >  1:4389114/1‑93 (MQ=255)
        tCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGcc                                                                                   <  2:3523468/93‑1 (MQ=255)
        tCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGcc                                                                                   >  2:550992/1‑93 (MQ=255)
        tCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGcc                                                                                   <  1:1359794/93‑1 (MQ=255)
        tCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGcc                                                                                   <  2:5336744/93‑1 (MQ=255)
           accacAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGa                                                                                <  1:2132448/93‑1 (MQ=255)
           accacAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGa                                                                                >  1:1234201/1‑93 (MQ=255)
            ccacAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGAt                                                                               <  1:5756148/93‑1 (MQ=255)
             cacaGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGAtt                                                                              <  2:1125683/93‑1 (MQ=255)
              acaGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTg                                                                             <  1:4137951/93‑1 (MQ=255)
              acaGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTg                                                                             <  2:5671833/93‑1 (MQ=255)
               caGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTgg                                                                            >  1:2395531/1‑93 (MQ=255)
                 gAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGcc                                                                          >  1:2802767/1‑93 (MQ=255)
                    gggCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGc                                                                      >  1:5636832/1‑93 (MQ=255)
                    gggCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGc                                                                      >  1:2573120/1‑93 (MQ=255)
                     ggCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCa                                                                     <  1:3880764/93‑1 (MQ=255)
                     ggCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTACCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCa                                                                     >  1:7370730/1‑93 (MQ=255)
                       cTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATg                                                                   >  1:5607305/1‑93 (MQ=255)
                       cTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATg                                                                   <  2:4831738/93‑1 (MQ=255)
                        ttCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATgg                                                                  <  2:4979028/93‑1 (MQ=255)
                        ttCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATgg                                                                  >  2:1174011/1‑93 (MQ=255)
                        ttCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATgg                                                                  >  1:151461/1‑93 (MQ=255)
                          ccTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGtt                                                                >  1:4733719/1‑93 (MQ=255)
                           cTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGttt                                                               >  1:3003195/1‑93 (MQ=255)
                            tttCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTg                                                              <  2:2626611/93‑1 (MQ=255)
                             ttCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTgc                                                             <  1:7102187/93‑1 (MQ=255)
                             ttCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTgc                                                             >  1:3400131/1‑93 (MQ=255)
                              tCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTgcg                                                            >  2:6068452/1‑93 (MQ=255)
                              tCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTgcg                                                            <  1:59537/93‑1 (MQ=255)
                                  aTCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGcgccgc                                                        <  1:1189022/93‑1 (MQ=255)
                                  aTCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGcgccgc                                                        >  2:2016088/1‑93 (MQ=255)
                                  aTCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGcgccgc                                                        <  1:1953144/93‑1 (MQ=255)
                                   tCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCt                                                       <  1:3718033/93‑1 (MQ=255)
                                    ccAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGGCATGGTTTGCGCCGCt                                                       >  1:2362898/1‑93 (MQ=255)
                                    ccAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTg                                                      <  2:470835/93‑1 (MQ=255)
                                      aGAAAGGAGATCACCCCCCCGCCGTTCCCCATTTCCTGAATCGCCTTTCCGGGTCAAAACCCGGAATGGCCCGCCTGGTTTGCGCCGCTGCCt                                                   <  1:2272594/93‑1 (MQ=255)
                                      aGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGcc                                                    >  2:3789288/1‑93 (MQ=255)
                                      aGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGcc                                                    <  1:1257202/93‑1 (MQ=255)
                                      aGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGcc                                                    >  2:6008597/1‑93 (MQ=255)
                                       gAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCt                                                   <  2:1829108/93‑1 (MQ=255)
                                       gAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCt                                                   <  2:647526/93‑1 (MQ=255)
                                       gAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCt                                                   >  2:6952878/1‑93 (MQ=255)
                                        aaaaGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCtt                                                  <  1:1067992/93‑1 (MQ=255)
                                         aaaGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTc                                                 <  2:6345928/93‑1 (MQ=255)
                                           agagATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCAc                                               >  2:1574984/1‑93 (MQ=255)
                                            gagaTCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACt                                              >  1:7215065/1‑93 (MQ=255)
                                            gagaTCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACt                                              <  1:6894067/93‑1 (MQ=255)
                                                tCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCtt                                          >  1:1028581/1‑93 (MQ=255)
                                                 cAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCttt                                         <  1:3790295/93‑1 (MQ=255)
                                                  aaCGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTc                                        >  2:2199997/1‑93 (MQ=255)
                                                  aaCGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTc                                        <  2:4437579/93‑1 (MQ=255)
                                                  aaCGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTc                                        >  2:4029860/1‑93 (MQ=255)
                                                   aCGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCt                                       >  1:405859/1‑93 (MQ=255)
                                                   aCGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCt                                       >  2:3875885/1‑93 (MQ=255)
                                                   aCGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCt                                       >  2:4324396/1‑93 (MQ=255)
                                                   aCGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCt                                       <  2:366604/93‑1 (MQ=255)
                                                    cGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTg                                      <  1:6526375/93‑1 (MQ=255)
                                                    cGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTg                                      <  2:1690066/93‑1 (MQ=255)
                                                     gCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGa                                     >  1:924373/1‑93 (MQ=255)
                                                      ccTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGAt                                    >  2:2797595/1‑93 (MQ=255)
                                                       cTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATc                                   >  1:6723953/1‑93 (MQ=255)
                                                       cTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATc                                   <  2:2874042/93‑1 (MQ=255)
                                                        tCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCt                                  <  2:809054/93‑1 (MQ=255)
                                                        tCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCt                                  <  2:3213217/93‑1 (MQ=255)
                                                        tCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCt                                  <  1:4257348/93‑1 (MQ=255)
                                                            cGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTg                              >  1:4836592/1‑93 (MQ=255)
                                                            cGGTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTg                              <  2:6153404/93‑1 (MQ=255)
                                                             ggTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTgg                             >  1:5199986/1‑93 (MQ=255)
                                                             ggTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTgg                             <  2:4669260/93‑1 (MQ=255)
                                                              gTGCCCAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGc                            <  2:2484919/93‑1 (MQ=255)
                                                                  ccAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTaa                        >  1:2902449/1‑93 (MQ=255)
                                                                   cAGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTaaa                       >  1:4627911/1‑93 (MQ=255)
                                                                    aGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTaaaa                      <  1:27455/93‑1 (MQ=255)
                                                                    aGTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTaaaa                      >  2:1351388/1‑93 (MQ=255)
                                                                     gTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGGCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAg                     >  1:7509057/1‑93 (MQ=255)
                                                                     gTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAg                     >  2:3353175/1‑93 (MQ=255)
                                                                     gTTCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAg                     >  1:705287/1‑93 (MQ=255)
                                                                      ttCCTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGc                    <  1:6608678/93‑1 (MQ=255)
                                                                        ccTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGccc                  <  2:1232902/93‑1 (MQ=255)
                                                                        ccTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGccc                  <  1:3042762/93‑1 (MQ=255)
                                                                         cTGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCg                 >  2:4246419/1‑93 (MQ=255)
                                                                          tGAATCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGc                >  2:7189498/1‑93 (MQ=255)
                                                                           gAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCa               >  1:5296065/1‑93 (MQ=255)
                                                                             aTCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCaaa             >  1:231874/1‑93 (MQ=255)
                                                                             aTCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCaaa             >  1:6137646/1‑93 (MQ=255)
                                                                              tCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAg            >  1:5695447/1‑93 (MQ=255)
                                                                              tCGCCTTTTCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAg            >  2:6854014/1‑93 (MQ=255)
                                                                                gCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCa          >  1:801388/1‑93 (MQ=255)
                                                                                   ttttCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGt       >  2:7231688/1‑93 (MQ=255)
                                                                                   ttttCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGt       >  2:7482022/1‑93 (MQ=255)
                                                                                   ttttCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGAACTTTTGGCGTAAAAGCCCGCAAAGCATGt       >  2:3450788/1‑93 (MQ=255)
                                                                                     ttCGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTg     <  2:5067872/93‑1 (MQ=255)
                                                                                       cGGTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGaa   <  2:2128691/93‑1 (MQ=255)
                                                                                        ggTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAc  >  2:6093577/1‑93 (MQ=255)
                                                                                        ggTGTCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAc  <  1:1790210/93‑1 (MQ=255)
                                                                                           |                                                                                          
TCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAC  >  NC_000913/4486497‑4486677

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