Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,897,501 C→A 16.6% D124Y (GAT→TAT)  ygcS ← putative MFS sugar transporter; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,897,5010CA16.6% 458.8 / 69.2 179D124Y (GAT→TAT) ygcSputative MFS sugar transporter; membrane protein
Reads supporting (aligned to +/- strand):  ref base C (62/85);  new base A (19/11);  total (83/96)
Fisher's exact test for biased strand distribution p-value = 4.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 186 total aligned reads displayed.

CATAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAA  >  NC_000913/2897402‑2897600
                                                                                                   |                                                                                                   
 aTAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCg                                                                                                    >  2:3787365/1‑100 (MQ=255)
  tAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGc                                                                                                   >  1:1368886/1‑100 (MQ=255)
   aGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGcc                                                                                                  <  2:2934441/100‑1 (MQ=255)
    gCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCt                                                                                                 >  1:3007790/1‑100 (MQ=255)
    gCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCt                                                                                                 >  2:1224843/1‑100 (MQ=255)
    gCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCt                                                                                                 >  2:4196261/1‑100 (MQ=255)
     ccTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTc                                                                                                >  2:2804099/1‑100 (MQ=255)
      cTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTcc                                                                                               <  2:1533218/100‑1 (MQ=255)
       tACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTccc                                                                                              <  2:2433733/100‑1 (MQ=255)
       tACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTccc                                                                                              >  1:1529090/1‑100 (MQ=255)
        aCGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCa                                                                                             <  2:1210641/100‑1 (MQ=255)
        aCGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAACTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTccc                                                                                              <  1:3467128/99‑1 (MQ=255)
         cGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAg                                                                                            <  1:4190948/100‑1 (MQ=255)
          tgTCCACACCACGCTGATTGCGCCCAGCAAAATACCGCGAGAGCGGCGCGGGGGTAACTCAGCCAGCAAGGCGTGACCGACTGAATAATCGCCTcc                                                                                               >  1:2546211/2‑96 (MQ=255)
          ggTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGa                                                                                           >  2:1903380/1‑100 (MQ=255)
           gTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGAc                                                                                          >  2:2856314/1‑100 (MQ=255)
           gTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATAGCCTCCCAGa                                                                                           <  1:922017/99‑1 (MQ=255)
               acacCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAAt                                                                                      >  2:2780403/1‑100 (MQ=255)
               acacCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAAt                                                                                      <  1:4105705/100‑1 (MQ=255)
                caccacGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATg                                                                                     <  2:1084076/100‑1 (MQ=255)
                caccacGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATg                                                                                     <  2:1958049/100‑1 (MQ=255)
                  ccacGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGcc                                                                                   <  2:1510335/100‑1 (MQ=255)
                     cGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAAt                                                                                <  2:1181434/100‑1 (MQ=255)
                     cGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAAt                                                                                >  2:195829/1‑100 (MQ=255)
                     cGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATAGCCTCCCAGACCAATGCCAAt                                                                                <  2:2353053/100‑1 (MQ=255)
                      gCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATc                                                                               >  1:3866406/1‑100 (MQ=255)
                       cTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCa                                                                              >  1:1479122/1‑100 (MQ=255)
                        tGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCaa                                                                             <  2:1735400/100‑1 (MQ=255)
                          aaTGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCaaaa                                                                           <  1:3005250/100‑1 (MQ=255)
                           aTGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAAt                                                                          >  1:1298056/1‑100 (MQ=255)
                           aTGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATAGCCTCCCAGACCAATGCCAATCAAAAt                                                                          >  1:493831/1‑100 (MQ=255)
                              cgcCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCACTCAAAATgcg                                                                       >  1:2803669/1‑100 (MQ=255)
                               gcCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATgcgc                                                                      <  2:3387763/100‑1 (MQ=255)
                                 ccAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATAGCCTCCCAGACCAATGCCAATCAAAATGCGCAg                                                                    >  2:2230416/1‑100 (MQ=255)
                                  cAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCa                                                                     <  1:4034441/98‑1 (MQ=255)
                                    gCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTcc                                                                 >  2:3076378/1‑100 (MQ=255)
                                    gCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTcc                                                                 <  1:1222101/100‑1 (MQ=255)
                                    gCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTcc                                                                 <  2:1469067/100‑1 (MQ=255)
                                    gCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATAGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTcc                                                                 >  2:2997331/1‑100 (MQ=255)
                                     cAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACTTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCa                                                                >  1:1679096/1‑100 (MQ=255)
                                     cAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCa                                                                >  1:2253885/1‑100 (MQ=255)
                                         aTACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATaa                                                            <  2:514897/100‑1 (MQ=255)
                                           aCCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGa                                                          >  1:1929359/1‑100 (MQ=255)
                                            ccGCGATGGCGGCGCGGGGAAAATTAAGCCAGCAAGGTGGGACATACTGAAAAATGGCCTCACAGACCAATGCAACACAAAAAGGGCAATCCAATAAGa                                                          >  1:3836448/1‑99 (MQ=255)
                                             cgcgATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATg                                                        <  1:2707545/100‑1 (MQ=255)
                                              gcggtgggcgcgcGGGGAGAATGCAGCCAGCTAGGTGTGACCGACTGAGTAAGCGCCTCCCAGACCCATGCCAATCAAAATGCGCAGTCCAATAAg                                                           <  1:841720/87‑1 (MQ=255)
                                                  tGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTg                                                   >  1:100655/1‑100 (MQ=255)
                                                   ggcggcGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTgg                                                  <  1:4155285/100‑1 (MQ=255)
                                                    gcggcgCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGc                                                 <  1:3330106/100‑1 (MQ=255)
                                                    gcggcgCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGc                                                 >  1:3992737/1‑100 (MQ=255)
                                                         gcgGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCgtggt                                            >  2:1212833/1‑100 (MQ=255)
                                                          cgGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCgtggtg                                           >  1:3144824/1‑100 (MQ=255)
                                                           ggggAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGtggtgg                                          <  2:3228150/100‑1 (MQ=255)
                                                            gggAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGc                                         >  1:1807521/1‑100 (MQ=255)
                                                             ggAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCa                                        >  2:2905216/1‑100 (MQ=255)
                                                              gAAAATTCAGCCAGCCAGATGTAACCTAATGACTACTCGCCTCACAGACCAATGCCACTCAACATGCGCAGTCCAATAAGATGCTCTTGCgtggttgt                                         <  1:352422/98‑4 (MQ=255)
                                                                  aTTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCaaaaaa                                   >  2:1305640/1‑100 (MQ=255)
                                                                  aTTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCaaaaaa                                   <  2:2474572/100‑1 (MQ=255)
                                                                  aTTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGAGCAGTCCAATAAGATGCTCTGGCGTGGTGGCaaaaaa                                   >  1:2891931/1‑100 (MQ=255)
                                                                   ttCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAAt                                  <  2:471586/100‑1 (MQ=255)
                                                                    tCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAATATGCTCTGGCGTGGTGGCAAAAAAtt                                 <  2:2604786/100‑1 (MQ=255)
                                                                     cagccagcAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTg                                <  2:1368890/100‑1 (MQ=255)
                                                                     cagccagcAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGAGGCAAAAAATTg                                <  1:1099927/100‑1 (MQ=255)
                                                                     cagccagcAAGGTGTGACCTACTGAATAATAGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTg                                >  2:3637203/1‑100 (MQ=255)
                                                                       gccagcAAGGTGTGACCTACTGAATGATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTa                              <  2:2827806/100‑1 (MQ=255)
                                                                       gccagcAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGATGGCAAAAAATTg                                >  1:1036640/1‑98 (MQ=255)
                                                                        ccagcAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTaa                             <  2:3126864/100‑1 (MQ=255)
                                                                         cagcAAGGTGTGACCTACTGAATAATAGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTaaa                            <  1:3875869/100‑1 (MQ=255)
                                                                              aGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGa                       >  2:679239/1‑100 (MQ=255)
                                                                               ggTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGaa                      <  2:265278/100‑1 (MQ=255)
                                                                                gtgtGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGaag                     <  2:2178536/100‑1 (MQ=255)
                                                                                gtgtGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGaag                     <  2:3931288/100‑1 (MQ=255)
                                                                                 tgtgACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGaagc                    <  2:3879667/100‑1 (MQ=255)
                                                                                  gtgACCTACTGACTAATCGCCTCCCAGACCAATGCCAACCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACTaagta                   >  1:1654929/1‑98 (MQ=255)
                                                                                  gtgACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGaag                     <  1:2280048/98‑1 (MQ=255)
                                                                                     aCCTACTGAATAATAGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGaagcaagc                >  2:3381063/1‑100 (MQ=255)
                                                                                      ccTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAAAAAGATGCTCTGGCATGGAAGAAAAAAATTGTAAAAACGAAGCAAGCg               <  1:2218608/100‑1 (MQ=255)
                                                                                      ccTACTGAATAATAGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCg               >  2:396119/1‑100 (MQ=255)
                                                                                      ccTACTGAATAATAGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCg               >  1:3807502/1‑100 (MQ=255)
                                                                                        tACTGAATAATCGCCTTCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTACAAACGAAGCAAGCGTa             <  1:2489970/100‑1 (MQ=255)
                                                                                        tACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGt              <  1:2089082/99‑1 (MQ=255)
                                                                                         aCTGAATAATAGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTaa            >  2:2321318/1‑100 (MQ=255)
                                                                                          cTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAAt           <  2:2546622/100‑1 (MQ=255)
                                                                                           tGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATc          >  1:764346/1‑100 (MQ=255)
                                                                                           tGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATc          >  1:1956375/1‑100 (MQ=255)
                                                                                            gAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCa         >  1:2494953/1‑100 (MQ=255)
                                                                                            gAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCa         >  1:1091474/1‑100 (MQ=255)
                                                                                            gAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCa         >  2:1645983/1‑100 (MQ=255)
                                                                                            gAATAATCGCATAACAGACCAATGCCAATCAAAATGCGCAGTCAAATAAGATGCTCTGGCGTGGTGGAAAAAAATTGTAAAAACGAAGCAAGAGTa             >  2:2003684/1‑96 (MQ=255)
                                                                                            gAATAATAGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGTGTAATCa         >  2:27561/1‑100 (MQ=255)
                                                                                             aataatCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAg        <  2:4507/100‑1 (MQ=255)
                                                                                             aataatAGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATc          <  2:1499089/98‑1 (MQ=255)
                                                                                              ataatCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGc       <  2:1441675/100‑1 (MQ=255)
                                                                                               taatCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAACCGAAGCAAGCGTAATCAGCa      <  1:2609613/100‑1 (MQ=255)
                                                                                                aatGGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCaa     >  1:3898894/1‑100 (MQ=255)
                                                                                                 atCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAAAAAGATGCTCTGGCGTGGTGACAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCaaa    >  1:2329014/1‑100 (MQ=255)
                                                                                                  tCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCaaaa   <  2:1491454/100‑1 (MQ=255)
                                                                                                  tCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCaaaa   >  1:2082090/1‑100 (MQ=255)
                                                                                                   cGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCaaaaa  >  2:3471397/1‑100 (MQ=255)
                                                                                                   cGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCaaaa   >  1:2397864/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
CATAGCCTACGGTCCACACCACGCTGAATGCGCCCAGCAAAATACCGCGATGGCGGCGCGGGGAAAATTCAGCCAGCAAGGTGTGACCTACTGAATAATCGCCTCCCAGACCAATGCCAATCAAAATGCGCAGTCCAATAAGATGCTCTGGCGTGGTGGCAAAAAATTGTAAAAACGAAGCAAGCGTAATCAGCAAAAA  >  NC_000913/2897402‑2897600

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