Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,522,256 G→C G591A (GGC→GCC)  mrcA → penicillin‑binding protein 1A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,522,2560GC83.5% 847.6 / 146.2 309G591A (GGC→GCC) mrcApenicillin‑binding protein 1A
Reads supporting (aligned to +/- strand):  ref base G (17/34);  new base C (108/150);  total (125/184)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

Only 100 of 316 total aligned reads displayed.

ACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCG  >  NC_000913/3522164‑3522347
                                                                                            |                                                                                           
   aaaCATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAAcgccgg                                                                                          >  1:3918476/1‑93 (MQ=255)
    aaCATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACgccggc                                                                                         >  1:1366613/1‑93 (MQ=255)
      cATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCtt                                                                                       <  1:7358822/93‑1 (MQ=255)
        ttGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGACTACGCGGTCATGGCGAACGCCGGCTTcc                                                                                     <  1:6391889/93‑1 (MQ=255)
         tGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCt                                                                                    <  2:206282/93‑1 (MQ=255)
         tGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCt                                                                                    <  2:6257222/93‑1 (MQ=255)
           tCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCtgg                                                                                  <  2:7003286/93‑1 (MQ=255)
             cacaCCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCtggtg                                                                                <  1:5923564/93‑1 (MQ=255)
               cacCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGa                                                                              <  2:443832/93‑1 (MQ=255)
               cacCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGa                                                                              >  2:6501671/1‑93 (MQ=255)
               cacCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGa                                                                              <  2:1088376/93‑1 (MQ=255)
               cacCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGa                                                                              >  2:3266546/1‑93 (MQ=255)
                acCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGAc                                                                             <  2:6607607/93‑1 (MQ=255)
                 ccGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGAcc                                                                            <  1:2816134/93‑1 (MQ=255)
                  cGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGAccc                                                                           >  1:656189/1‑93 (MQ=255)
                  cGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGAccc                                                                           <  1:3248847/93‑1 (MQ=255)
                  cGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGAccc                                                                           <  2:5417743/93‑1 (MQ=255)
                    aaTCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGt                                                                         <  1:4634887/93‑1 (MQ=255)
                    aaTCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGt                                                                         >  1:6812051/1‑93 (MQ=255)
                    aaTCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGt                                                                         >  2:5047086/1‑93 (MQ=255)
                     aTCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTg                                                                        <  1:1045969/93‑1 (MQ=255)
                     aTCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTg                                                                        <  1:2547469/93‑1 (MQ=255)
                      tCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTgg                                                                       >  1:2716077/1‑93 (MQ=255)
                        gctggcgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGtt                                                                     <  1:3696358/93‑1 (MQ=255)
                         ctggcgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGttt                                                                    <  2:6432148/93‑1 (MQ=255)
                          tggcgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTa                                                                   <  1:2275591/93‑1 (MQ=255)
                            gcgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATc                                                                 <  2:1306165/93‑1 (MQ=255)
                            gcgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATc                                                                 >  2:2234066/1‑93 (MQ=255)
                            gcgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATc                                                                 <  2:231533/93‑1 (MQ=255)
                            gcgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATc                                                                 <  2:3552069/93‑1 (MQ=255)
                             cgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCa                                                                <  1:5442628/93‑1 (MQ=255)
                             cgctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCa                                                                <  1:1326624/93‑1 (MQ=255)
                              gctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAg                                                               >  1:1964850/1‑93 (MQ=255)
                               ctggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGc                                                              >  2:73412/1‑93 (MQ=255)
                                tggGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCa                                                             <  1:6027808/93‑1 (MQ=255)
                                  ggTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCaaa                                                           >  1:6055162/1‑93 (MQ=255)
                                   gTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCaaaa                                                          <  1:29350/93‑1 (MQ=255)
                                    ttCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAAt                                                         <  1:5696363/93‑1 (MQ=255)
                                    ttCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAAt                                                         <  1:5494634/93‑1 (MQ=255)
                                    ttCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAAt                                                         >  1:5626008/1‑93 (MQ=255)
                                       aGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGa                                                      >  1:2158325/1‑93 (MQ=255)
                                       aGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGa                                                      <  2:2397451/93‑1 (MQ=255)
                                        gCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGaa                                                     <  1:5383855/93‑1 (MQ=255)
                                        gCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGaa                                                     >  1:5256900/1‑93 (MQ=255)
                                         cGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGaaa                                                    <  1:4081914/93‑1 (MQ=255)
                                          gTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGaaaa                                                   <  2:5139323/93‑1 (MQ=255)
                                           tCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAAc                                                  <  1:4259310/93‑1 (MQ=255)
                                           tCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAAc                                                  <  2:1133964/93‑1 (MQ=255)
                                              ttCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGAt                                               >  1:5045223/1‑93 (MQ=255)
                                               tCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATc                                              <  2:528200/93‑1 (MQ=255)
                                                cACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCa                                             <  2:6720346/93‑1 (MQ=255)
                                                 aCCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAg                                            >  2:3725706/1‑93 (MQ=255)
                                                  ccccAACGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAgg                                           <  2:4340629/93‑1 (MQ=255)
                                                   cccAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAggg                                          <  1:6654862/93‑1 (MQ=255)
                                                    ccAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGggc                                         >  1:4913826/1‑93 (MQ=255)
                                                    ccAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGggc                                         >  2:1781087/1‑93 (MQ=255)
                                                     cAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGggcg                                        <  2:802277/93‑1 (MQ=255)
                                                     cAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGggcg                                        >  2:1665662/1‑93 (MQ=255)
                                                     cAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGggcg                                        <  1:1248992/93‑1 (MQ=255)
                                                     cAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGggcg                                        <  1:1636757/93‑1 (MQ=255)
                                                      aaTGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGggcgg                                       >  2:4647261/1‑93 (MQ=255)
                                                      aaTGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGggcgg                                       >  2:6032790/1‑93 (MQ=255)
                                                       aTGCAGATGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGggcggc                                      >  2:2684351/1‑93 (MQ=255)
                                                        tGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGgcggcg                                     >  2:2132197/1‑93 (MQ=255)
                                                         gCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGt                                    <  1:2973123/93‑1 (MQ=255)
                                                          cAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGGaagcc                                                         >  2:3947169/1‑68 (MQ=255)
                                                          cAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGGaagcc                                                         <  1:3947169/71‑4 (MQ=255)
                                                          cAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTg                                   <  1:2448362/93‑1 (MQ=255)
                                                           aGGTGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGa                                  <  2:707902/93‑1 (MQ=255)
                                                              tGGCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGAttt                               >  2:3224256/1‑93 (MQ=255)
                                                               ggCGCGCGGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGAtttt                              <  2:4521362/93‑1 (MQ=255)
                                                                   cgcgGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGaa                          <  2:7154817/93‑1 (MQ=255)
                                                                   cgcgGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGaa                          <  1:1572756/93‑1 (MQ=255)
                                                                     cgGCTACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGc                        <  1:6965577/93‑1 (MQ=255)
                                                                      ggCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCg                       >  2:520165/1‑93 (MQ=255)
                                                                         tACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGaaa                    >  2:4915894/1‑93 (MQ=255)
                                                                         tACGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGaaa                    >  2:5582629/1‑93 (MQ=255)
                                                                          aCGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAAc                   <  2:1900968/93‑1 (MQ=255)
                                                                          aCGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAAc                   <  1:2663404/93‑1 (MQ=255)
                                                                          aCGCGGTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAAc                   <  1:7089167/93‑1 (MQ=255)
                                                                           cgcgGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAAcc                  <  1:843032/93‑1 (MQ=255)
                                                                              ggTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGaa               <  2:7314034/93‑1 (MQ=255)
                                                                               gTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGaaa              <  2:1459306/93‑1 (MQ=255)
                                                                               gTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGaaa              <  1:3644413/93‑1 (MQ=255)
                                                                               gTCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGaaa              >  2:467996/1‑93 (MQ=255)
                                                                                tCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAg             <  2:7468557/93‑1 (MQ=255)
                                                                                tCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAg             >  2:4214341/1‑93 (MQ=255)
                                                                                tCATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAg             >  1:4006508/1‑93 (MQ=255)
                                                                                 cATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGt            <  2:4814768/93‑1 (MQ=255)
                                                                                 cATGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGt            <  2:2451507/93‑1 (MQ=255)
                                                                                  aTGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTa           <  1:2323047/93‑1 (MQ=255)
                                                                                  aTGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTa           <  1:5789629/93‑1 (MQ=255)
                                                                                   tGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAg          >  2:5877201/1‑93 (MQ=255)
                                                                                   tGGCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAg          <  2:7526430/93‑1 (MQ=255)
                                                                                    nncGACCGCCGGCTTCCTGGTGGCCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGcc        <  1:2235296/91‑1 (MQ=255)
                                                                                     gCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGcc        <  1:4861132/93‑1 (MQ=255)
                                                                                     gCGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGcc        <  2:5750736/93‑1 (MQ=255)
                                                                                      cGAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCt       <  2:2975499/93‑1 (MQ=255)
                                                                                       gAACGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTg      >  2:6824585/1‑93 (MQ=255)
                                                                                         aCGCCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGcc    <  1:4456201/93‑1 (MQ=255)
                                                                                            |                                                                                           
ACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCCTGCCCG  >  NC_000913/3522164‑3522347

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