Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,636,465 T→C intergenic (‑296/+14) b1551 ← / ← cspI hypothetical protein/cold shock‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,636,4650TC85.0% 506.4 / 79.8 215intergenic (‑296/+14)b1551/cspIhypothetical protein/cold shock‑like protein
Reads supporting (aligned to +/- strand):  ref base T (16/16);  new base C (103/79);  total (119/96)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 216 total aligned reads displayed.

TGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAG  >  NC_000913/1636373‑1636555
                                                                                            |                                                                                          
    ccGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAata                                                                                        >  1:5174438/1‑93 (MQ=255)
    ccGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAata                                                                                        <  2:2940870/93‑1 (MQ=255)
     cGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAatat                                                                                       >  2:1262254/1‑93 (MQ=255)
     cGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAatat                                                                                       >  2:1706384/1‑93 (MQ=255)
      gAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATAtt                                                                                      >  2:694651/1‑93 (MQ=255)
      gAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATAtt                                                                                      >  2:5687404/1‑93 (MQ=255)
          ggATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCt                                                                                  <  1:4574287/93‑1 (MQ=255)
           gATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTa                                                                                 <  2:6301623/93‑1 (MQ=255)
            aTTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTAc                                                                                >  2:6156921/1‑93 (MQ=255)
            aTTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTAc                                                                                >  2:308586/1‑93 (MQ=255)
            aTTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTAc                                                                                <  2:5852068/93‑1 (MQ=255)
             ttCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTAcc                                                                               >  1:1494835/1‑93 (MQ=255)
              tCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCt                                                                              >  1:1002047/1‑93 (MQ=255)
               cATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTc                                                                             <  1:2595475/93‑1 (MQ=255)
               cATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTc                                                                             <  2:596383/93‑1 (MQ=255)
                aTGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCa                                                                            <  2:2118802/93‑1 (MQ=255)
                 tGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCaa                                                                           >  1:4642667/1‑93 (MQ=255)
                  gTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCaaa                                                                          >  1:7478716/1‑93 (MQ=255)
                   tCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAg                                                                         <  2:2443207/93‑1 (MQ=255)
                    cAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAgc                                                                        <  1:4114547/93‑1 (MQ=255)
                    cAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAgc                                                                        >  2:5638860/1‑93 (MQ=255)
                         aCATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCAc                                                                   >  2:6774173/1‑93 (MQ=255)
                         aCATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCAc                                                                   <  2:372037/93‑1 (MQ=255)
                         aCATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCAc                                                                   >  2:152622/1‑93 (MQ=255)
                           aTCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTa                                                                 >  1:5929728/1‑93 (MQ=255)
                           aTCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTa                                                                 >  1:787601/1‑93 (MQ=255)
                             cAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACa                                                               <  1:6553952/93‑1 (MQ=255)
                             cAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACa                                                               >  1:2786888/1‑93 (MQ=255)
                              aaaaGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACAt                                                              <  1:5394918/93‑1 (MQ=255)
                                aaGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGa                                                            >  2:673653/1‑93 (MQ=255)
                                aaGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGa                                                            >  1:4995473/1‑93 (MQ=255)
                                 aGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGaa                                                           >  1:3630244/1‑93 (MQ=255)
                                 aGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGaa                                                           <  2:4243968/93‑1 (MQ=255)
                                  gatgatGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAAc                                                          >  1:6621226/1‑93 (MQ=255)
                                  gatgatGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAAc                                                          >  1:5342212/1‑93 (MQ=255)
                                   atgatgTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACa                                                         >  2:4513678/1‑93 (MQ=255)
                                    tgatgTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAg                                                        >  2:4691442/1‑93 (MQ=255)
                                     gatgTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGc                                                       <  2:2223447/93‑1 (MQ=255)
                                     gatgTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGc                                                       <  2:5126215/93‑1 (MQ=255)
                                      atgTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCg                                                      <  2:3982047/93‑1 (MQ=255)
                                       tgTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCgg                                                     <  2:6572867/93‑1 (MQ=255)
                                          tCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAg                                                  >  2:6355898/1‑93 (MQ=255)
                                          tCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAg                                                  >  2:1877073/1‑93 (MQ=255)
                                          tCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAg                                                  <  1:6760588/93‑1 (MQ=255)
                                          tCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAg                                                  >  2:4797169/1‑93 (MQ=255)
                                             ataCTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGAc                                               >  2:4085653/1‑93 (MQ=255)
                                               aCTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCt                                             >  1:1111837/1‑93 (MQ=255)
                                                cTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCtt                                            <  2:7176507/93‑1 (MQ=255)
                                                cTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCtt                                            <  1:691398/93‑1 (MQ=255)
                                                 tGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCttt                                           <  2:1083168/93‑1 (MQ=255)
                                                 tGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCttt                                           <  2:5439007/93‑1 (MQ=255)
                                                 tGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCttt                                           >  2:4985994/1‑93 (MQ=255)
                                                  gggAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTa                                          <  1:3304438/93‑1 (MQ=255)
                                                  gggAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTa                                          <  1:3608791/93‑1 (MQ=255)
                                                   ggAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGTGCCACTACATGAACAGCGGCAGGACCTTTAg                                         <  1:4492178/93‑1 (MQ=255)
                                                     aaCCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGt                                       >  1:3415055/1‑93 (MQ=255)
                                                      aCCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTc                                      >  2:3350340/1‑93 (MQ=255)
                                                      aCCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTc                                      >  2:6910295/1‑93 (MQ=255)
                                                      aCCATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTc                                      >  2:3563055/1‑93 (MQ=255)
                                                        cATGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCg                                    >  1:932359/1‑93 (MQ=255)
                                                         aTGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGt                                   <  2:2997125/93‑1 (MQ=255)
                                                          tGACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGtt                                  <  2:2017864/93‑1 (MQ=255)
                                                           gACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTtc                                 <  1:1634273/93‑1 (MQ=255)
                                                           gACAACACGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTtc                                 <  1:1345175/93‑1 (MQ=255)
                                                           gACAACACGGGCATCTAAAGTGAATATGGCTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTtc                                 <  2:867881/93‑1 (MQ=255)
                                                              aacaCGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCa                              >  1:3078977/1‑93 (MQ=255)
                                                              aacaCGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCa                              <  1:5647105/93‑1 (MQ=255)
                                                               acacGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCaa                             >  2:2408384/1‑93 (MQ=255)
                                                               acacGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCaa                             <  2:3856914/93‑1 (MQ=255)
                                                               acacGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCaa                             <  2:156625/93‑1 (MQ=255)
                                                                cacGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAAt                            <  1:5982482/93‑1 (MQ=255)
                                                                 acGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATa                           >  2:7099387/1‑93 (MQ=255)
                                                                  cGGGCATCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATAc                          >  1:2498575/1‑93 (MQ=255)
                                                                    ggCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCa                        >  1:3934212/1‑93 (MQ=255)
                                                                       aTCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAAt                     >  1:6476821/1‑93 (MQ=255)
                                                                       aTCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAAt                     >  1:2068702/1‑93 (MQ=255)
                                                                       aTCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAAt                     >  1:7416097/1‑93 (MQ=255)
                                                                        tCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAAtt                    >  1:4318896/1‑93 (MQ=255)
                                                                        tCTAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAAtt                    >  2:1981139/1‑93 (MQ=255)
                                                                          tAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCa                  >  2:345459/1‑93 (MQ=255)
                                                                          tAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCa                  >  1:4779346/1‑93 (MQ=255)
                                                                          tAAAGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCa                  >  1:6258647/1‑93 (MQ=255)
                                                                           aaaGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACATCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCaa                 >  1:1188316/1‑93 (MQ=255)
                                                                           aaaGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCaa                 <  2:3779336/93‑1 (MQ=255)
                                                                            aaGTGAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAAc                <  1:3219742/93‑1 (MQ=255)
                                                                                gAATATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTcc            >  1:4245393/1‑93 (MQ=255)
                                                                                 aaTATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCt           <  2:5317469/93‑1 (MQ=255)
                                                                                 aaTATGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCt           <  1:5571272/93‑1 (MQ=255)
                                                                                  atatGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTg          >  1:4175203/1‑93 (MQ=255)
                                                                                   tatGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGa         <  1:5853895/93‑1 (MQ=255)
                                                                                   tatGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGa         >  1:2430254/1‑93 (MQ=255)
                                                                                   tatGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGa         <  1:2254139/93‑1 (MQ=255)
                                                                                   tatGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGa         >  2:4920865/1‑93 (MQ=255)
                                                                                   tatGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGa         <  1:7215423/93‑1 (MQ=255)
                                                                                    atGGTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGAt        <  1:402213/93‑1 (MQ=255)
                                                                                     tGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGAtt       <  2:2529105/93‑1 (MQ=255)
                                                                                       gTTTGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATtct     >  2:1740268/1‑93 (MQ=255)
                                                                                        tttGCAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATtctc    <  1:1762971/93‑1 (MQ=255)
                                                                                         ttGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCa   >  2:436520/1‑93 (MQ=255)
                                                                                          tGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAg  >  2:737099/1‑93 (MQ=255)
                                                                                            |                                                                                          
TGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAG  >  NC_000913/1636373‑1636555

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