Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,832,048 A→G R39R (CGT→CGC)  yicI ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,0480AG86.1% 734.0 / 98.6 281R39R (CGT→CGC) yicIhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/21);  new base G (114/127);  total (132/149)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 293 total aligned reads displayed.

TCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAAT  >  NC_000913/3831956‑3832137
                                                                                            |                                                                                         
tCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCa                                                                                           <  1:1226308/93‑1 (MQ=255)
   aTCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCgg                                                                                        <  2:166713/93‑1 (MQ=255)
    tCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCggg                                                                                       <  2:2247463/93‑1 (MQ=255)
    tCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCggg                                                                                       <  1:6226098/93‑1 (MQ=255)
    tCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCggg                                                                                       <  2:130560/93‑1 (MQ=255)
     ccGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCgggg                                                                                      <  2:812855/93‑1 (MQ=255)
     ccGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCgggg                                                                                      >  1:4192453/1‑93 (MQ=255)
       gCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCgggggg                                                                                    >  1:2812380/1‑93 (MQ=255)
       gCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCgggggg                                                                                    <  2:615561/93‑1 (MQ=255)
       ccacaCCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCgggggg                                                                                    <  2:2654407/92‑1 (MQ=255)
        cacaCCGACAATCCCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGgc                                                                                   <  1:1758679/93‑1 (MQ=255)
         acacCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGgca                                                                                  <  1:85/93‑1 (MQ=255)
          cacCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGgcag                                                                                 >  1:6559417/1‑93 (MQ=255)
            ccGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGgcagca                                                                               >  1:502898/1‑93 (MQ=255)
             cGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCAt                                                                              >  2:5586367/1‑93 (MQ=255)
             cGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCAt                                                                              <  2:431599/93‑1 (MQ=255)
             cGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCAt                                                                              >  2:5754683/1‑93 (MQ=255)
                  aTACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCCGCATAGacc                                                                         >  1:6106845/1‑93 (MQ=255)
                  aTACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGacc                                                                         <  1:4869015/93‑1 (MQ=255)
                    aCCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGaccac                                                                       >  1:6472713/1‑93 (MQ=255)
                    aCCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGaccac                                                                       >  2:475363/1‑93 (MQ=255)
                    aCCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGaccac                                                                       >  2:2788655/1‑93 (MQ=255)
                      cTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGAccacca                                                                     <  1:721484/93‑1 (MQ=255)
                         ccTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCAttt                                                                  <  1:216519/93‑1 (MQ=255)
                            gtgGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCAt                                                               >  2:1509192/1‑93 (MQ=255)
                            gtgGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCAt                                                               <  2:298929/93‑1 (MQ=255)
                             tgGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCAtt                                                              >  1:5618101/1‑93 (MQ=255)
                              ggggAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTa                                                             >  1:3555754/1‑93 (MQ=255)
                              ggggAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTa                                                             >  1:3027460/1‑93 (MQ=255)
                                ggAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATc                                                           <  1:197295/93‑1 (MQ=255)
                                 gagaAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATcc                                                          <  2:2023447/93‑1 (MQ=255)
                                  agaAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTTCGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCct                                                         <  2:934349/93‑1 (MQ=255)
                                  agaAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCct                                                         <  2:2968068/93‑1 (MQ=255)
                                   gaAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCACAGACCACCATTTCATTATCctg                                                        >  2:3799583/1‑93 (MQ=255)
                                    aaaaaGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCctgc                                                       <  1:626335/93‑1 (MQ=255)
                                     aaaaGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCctgct                                                      <  1:3971027/93‑1 (MQ=255)
                                      aaaGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCctgctg                                                     >  1:1807415/1‑93 (MQ=255)
                                        aGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGtt                                                   >  1:1093287/1‑93 (MQ=255)
                                        aGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGtt                                                   >  2:5863663/1‑93 (MQ=255)
                                          cgcAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCa                                                 >  1:4502435/1‑93 (MQ=255)
                                           gcAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCaa                                                >  2:6071232/1‑93 (MQ=255)
                                            cAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAAc                                               >  2:1845105/1‑93 (MQ=255)
                                             aaCGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAAcc                                              <  2:583138/93‑1 (MQ=255)
                                             aaCGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAAcc                                              <  1:4694375/93‑1 (MQ=255)
                                             aaCGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAAcc                                              <  2:2072896/93‑1 (MQ=255)
                                              aCGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCt                                             <  1:4134120/93‑1 (MQ=255)
                                              aCGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCt                                             >  2:77080/1‑93 (MQ=255)
                                                gTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCg                                           >  1:4797634/1‑93 (MQ=255)
                                                 taaataaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGa                                          <  1:234932/93‑1 (MQ=255)
                                                  aaataaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGaa                                         <  1:5368242/93‑1 (MQ=255)
                                                     taaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAacac                                      >  1:2916323/1‑93 (MQ=255)
                                                      aaaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACAcc                                     <  2:3998523/93‑1 (MQ=255)
                                                       aaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCt                                    >  1:5150497/1‑93 (MQ=255)
                                                       aaGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCt                                    >  1:3681553/1‑93 (MQ=255)
                                                        aGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTg                                   >  2:1077520/1‑93 (MQ=255)
                                                        aGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTg                                   >  2:4121787/1‑93 (MQ=255)
                                                        aGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTg                                   >  1:1545507/1‑93 (MQ=255)
                                                         ggCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGa                                  <  1:329096/93‑1 (MQ=255)
                                                          gCGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGaa                                 <  2:3249478/93‑1 (MQ=255)
                                                           cGTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAg                                <  1:5124795/93‑1 (MQ=255)
                                                            gTATCAAGCTGGCAGGTGCGTTCACGCACATCACGGGGGGCAGCATAGACCCCCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGc                               <  2:3279086/93‑1 (MQ=255)
                                                            gTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGc                               >  1:3895236/1‑93 (MQ=255)
                                                            gTATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGc                               <  2:2577928/93‑1 (MQ=255)
                                                             tATCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCg                              >  2:5179288/1‑93 (MQ=255)
                                                              aTCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCgg                             <  2:1029812/93‑1 (MQ=255)
                                                              aTCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCgg                             <  2:6342137/93‑1 (MQ=255)
                                                               tCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCggg                            >  2:3299790/1‑93 (MQ=255)
                                                               tCAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCggg                            >  2:5895518/1‑93 (MQ=255)
                                                                cAAGCTGCCAGGTACGTTCCCGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGt                           >  1:6450922/1‑93 (MQ=255)
                                                                cAAGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGt                           <  2:2460165/93‑1 (MQ=255)
                                                                 aaGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTg                          >  2:5929773/1‑93 (MQ=255)
                                                                 aaGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTg                          >  1:575910/1‑93 (MQ=255)
                                                                  aGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGtca                         <  2:6421235/93‑3 (MQ=255)
                                                                  aGCTGCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGa                         >  1:1595541/1‑93 (MQ=255)
                                                                  aGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGa                         >  1:6017976/1‑93 (MQ=255)
                                                                      gCCAGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCa                     >  1:2602437/1‑93 (MQ=255)
                                                                         aGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGACGCGGGCGAATCAAAt                  >  1:630631/1‑93 (MQ=255)
                                                                         aGGTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAAt                  >  1:1704725/1‑93 (MQ=255)
                                                                          ggTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAAtt                 >  2:1157667/1‑93 (MQ=255)
                                                                          ggTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAAtt                 <  1:3462858/93‑1 (MQ=255)
                                                                           gTACGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTg                <  1:3746557/93‑1 (MQ=255)
                                                                            tACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGa               <  1:212670/93‑1 (MQ=255)
                                                                              cGTTCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAgg             >  2:2118676/1‑93 (MQ=255)
                                                                                ttCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGcc           >  1:5824879/1‑93 (MQ=255)
                                                                                ttCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGca           >  2:1416392/1‑92 (MQ=255)
                                                                                gtCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGcc           <  2:5547854/92‑1 (MQ=255)
                                                                                 tCACGCACATCGCGGGGGTCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCa          <  1:3077752/93‑1 (MQ=255)
                                                                                 tCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCa          <  2:4510603/93‑1 (MQ=255)
                                                                                 tCACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCa          <  1:4007814/93‑1 (MQ=255)
                                                                                  cACGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAg         <  2:354418/93‑1 (MQ=255)
                                                                                    cGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGt       <  2:1958779/93‑1 (MQ=255)
                                                                                    cGCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGt       <  2:4933928/93‑1 (MQ=255)
                                                                                     gCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGtt      >  2:3513379/1‑93 (MQ=255)
                                                                                     gCACATCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGtt      <  2:3127528/93‑1 (MQ=255)
                                                                                       acaTCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGGGAATCAAATTGAGGCCAGGTTGa    >  1:4557178/1‑93 (MQ=255)
                                                                                       acaTCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGa    >  2:5694816/1‑93 (MQ=255)
                                                                                        caTCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGaa   <  1:4400595/93‑1 (MQ=255)
                                                                                        caTCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGGTGaa   >  2:3431223/1‑93 (MQ=255)
                                                                                         aTCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAAt  <  2:5204702/93‑1 (MQ=255)
                                                                                         aTCGCGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAAt  <  1:6374935/93‑1 (MQ=255)
                                                                                            |                                                                                         
TCAATCCGCACACCGACAATACCTTCCTGTGGGGAGAAAAAGCGCAACGTAAATAAAGGCGTATCAAGCTGCCAGGTACGTTCACGCACATCACGGGGGGCAGCATAGACCACCATTTCATTATCCTGCTGTTCAACCTCGAACACCTGAAGCGGGTGAATCAAATTGAGGCCAGGTTGAAT  >  NC_000913/3831956‑3832137

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