Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,580,094 C→T D233N (GAT→AAT)  ydeM ← putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,0940CT96.4% 447.2 / 2.6 138D233N (GAT→AAT) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (67/66);  total (69/69)
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.99e-01

Only 100 of 138 total aligned reads displayed.

CCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAT  >  NC_000913/1579998‑1580189
                                                                                                |                                                                                               
 cTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttt                                                                                             <  2:1276373/100‑1 (MQ=255)
 cTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttt                                                                                             <  2:1437126/100‑1 (MQ=255)
  ttGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtttt                                                                                            >  2:334396/1‑100 (MQ=255)
  ttGCAGGACTCCTAGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtttt                                                                                            >  1:2152576/1‑100 (MQ=255)
   tGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttttt                                                                                           >  2:886718/1‑100 (MQ=255)
   tGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttttt                                                                                           <  2:729032/100‑1 (MQ=255)
   tGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtttt                                                                                            >  1:1864606/1‑99 (MQ=255)
     cAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAGAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTaa                                                                                         >  1:1899652/1‑100 (MQ=255)
     cAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTaa                                                                                         >  1:1220854/1‑100 (MQ=255)
     cAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTaa                                                                                         <  2:3245369/100‑1 (MQ=255)
      aGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAc                                                                                        <  2:2887627/100‑1 (MQ=255)
      aGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAc                                                                                        <  2:3204697/100‑1 (MQ=255)
       ggACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAAcc                                                                                       >  1:565866/1‑100 (MQ=255)
        gACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAccc                                                                                      <  2:3016436/100‑1 (MQ=255)
           tCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCAtt                                                                                   <  2:1750229/100‑1 (MQ=255)
           tCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCAtt                                                                                   <  2:1043354/100‑1 (MQ=255)
            ccTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTg                                                                                  >  2:1039296/1‑100 (MQ=255)
             cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                 >  2:1907578/1‑100 (MQ=255)
             cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                 >  1:2557996/1‑100 (MQ=255)
             cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                 <  2:497550/100‑1 (MQ=255)
             cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                 <  1:3412851/100‑1 (MQ=255)
              tGGAAAATACAACTGGTATTCCCATGCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGAGGAAAATTTCACCCACATTGTTTTTAACCCATTGCa                                                                                <  2:2026599/100‑1 (MQ=255)
              tGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCa                                                                                <  2:2023851/100‑1 (MQ=255)
               ggAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCAt                                                                               >  2:1850120/1‑100 (MQ=255)
                 aaaaTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaa                                                                             >  1:2493244/1‑100 (MQ=255)
                    aTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaaaa                                                                          >  1:2588995/1‑100 (MQ=255)
                     tACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAAt                                                                         >  1:1727455/1‑100 (MQ=255)
                     tACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAAt                                                                         <  1:2661396/100‑1 (MQ=255)
                      aCAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATg                                                                        <  2:1054522/100‑1 (MQ=255)
                       cAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATgg                                                                       >  2:310844/1‑100 (MQ=255)
                       cAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATgg                                                                       >  1:853731/1‑100 (MQ=255)
                        aaCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGt                                                                      >  2:1122798/1‑100 (MQ=255)
                         aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                     >  1:1987170/1‑100 (MQ=255)
                         aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                     <  2:310073/100‑1 (MQ=255)
                          cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                    <  2:522116/100‑1 (MQ=255)
                          cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTg                                                                    <  2:1478071/100‑1 (MQ=255)
                           tGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCAAAAAAATAAAAATGGTTg                                                                    <  1:140805/99‑1 (MQ=255)
                                 gCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaaa                                                             >  2:2557047/1‑100 (MQ=255)
                                 gCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaaa                                                             <  1:766013/100‑1 (MQ=255)
                                  cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAAc                                                            <  2:2094558/100‑1 (MQ=255)
                                  cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAAc                                                            >  2:1694507/1‑100 (MQ=255)
                                  cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAAc                                                            >  2:1572826/1‑100 (MQ=255)
                                   ccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACt                                                           >  1:759372/1‑100 (MQ=255)
                                     aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCAGATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg                                                         >  1:2229848/1‑100 (MQ=255)
                                      ttCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGc                                                        <  2:626652/100‑1 (MQ=255)
                                      ttCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGc                                                        <  2:2670179/100‑1 (MQ=255)
                                       tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGcc                                                       <  1:1572804/100‑1 (MQ=255)
                                       tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGcc                                                       <  2:525885/100‑1 (MQ=255)
                                       tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGcc                                                       <  1:3471779/100‑1 (MQ=255)
                                        ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                      >  2:2374196/1‑100 (MQ=255)
                                        ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                      <  2:2107180/100‑1 (MQ=255)
                                        ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTCTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                      <  1:2267159/100‑1 (MQ=255)
                                         cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCAt                                                     >  1:611716/1‑100 (MQ=255)
                                          ccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGACGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATa                                                    >  2:2447039/1‑100 (MQ=255)
                                             aaaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGc                                                 >  2:1353386/1‑100 (MQ=255)
                                                aaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGt                                              <  2:1465359/100‑1 (MQ=255)
                                                  cGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTgg                                            >  2:2392818/1‑100 (MQ=255)
                                                   ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTggg                                           <  2:1987199/100‑1 (MQ=255)
                                                   ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTggg                                           <  1:2047101/100‑1 (MQ=255)
                                                    gCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGa                                          >  2:2747946/1‑100 (MQ=255)
                                                     cTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAg                                         >  2:3250879/1‑100 (MQ=255)
                                                      tGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTGGACATAAACTTGCCATAAGCCGTGGGAgg                                        >  1:3342888/1‑100 (MQ=255)
                                                       gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                       <  2:1905300/100‑1 (MQ=255)
                                                       gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                       >  1:2534938/1‑100 (MQ=255)
                                                       gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                       >  2:1872378/1‑100 (MQ=255)
                                                          aaaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcaca                                    <  2:2229880/100‑1 (MQ=255)
                                                               tttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaa                               >  2:2285177/1‑100 (MQ=255)
                                                                ttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaaa                              <  2:1189612/100‑1 (MQ=255)
                                                                ttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaaa                              <  2:1193639/100‑1 (MQ=255)
                                                                ttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaaa                              <  2:355528/100‑1 (MQ=255)
                                                                  cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                            <  2:1157568/100‑1 (MQ=255)
                                                                  cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                            >  1:3159037/1‑100 (MQ=255)
                                                                  cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                            >  2:2649113/1‑100 (MQ=255)
                                                                  cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                            >  2:2671234/1‑100 (MQ=255)
                                                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                           <  2:1392266/100‑1 (MQ=255)
                                                                     cttGACGGACGAAAATTTCACCCACATCGATACTAACCCCTTGCAAAAAAACGATCGACCTCAACTAACCACAAGCCCTGGAACCCACAGAAAAAACa                           >  1:444334/3‑98 (MQ=255)
                                                                     aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                         <  2:724267/100‑1 (MQ=255)
                                                                     aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                         <  1:2838957/100‑1 (MQ=255)
                                                                      cTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                         <  1:1072596/99‑1 (MQ=255)
                                                                      cTGACGGATGAAAATTTCACCCACATTGTTTTTAACACATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATg                        <  1:1918714/100‑1 (MQ=255)
                                                                      aatgtgggtGAAATTTTCATCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATg                        >  2:1272771/9‑100 (MQ=255)
                                                                        gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGAt                      >  2:1139149/1‑100 (MQ=255)
                                                                        gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGAt                      >  1:954693/1‑100 (MQ=255)
                                                                         aCGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATa                     <  2:2180275/100‑1 (MQ=255)
                                                                          cGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATAc                    >  2:222763/1‑100 (MQ=255)
                                                                          cGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATAc                    <  1:1039283/100‑1 (MQ=255)
                                                                           tgATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACg                   <  1:2750878/99‑1 (MQ=255)
                                                                               gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt               >  2:2844247/1‑100 (MQ=255)
                                                                               gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt               >  2:913238/1‑100 (MQ=255)
                                                                               gAAAATTTCACCCACATTGATTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGAAAAAAAAAAAAAAATCAATGATACGGAAt               <  1:244097/100‑1 (MQ=255)
                                                                                aaaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtg              >  2:1903645/1‑100 (MQ=255)
                                                                                 aaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgt             >  1:1456227/1‑100 (MQ=255)
                                                                                  aaTTTCACCCACATTGTTTTTAACCCAGTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgtg            <  2:390597/100‑1 (MQ=255)
                                                                                     ttCACCCACATTGTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtct        <  1:3321769/100‑1 (MQ=255)
                                                                                       cACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtctc       >  2:2707508/1‑100 (MQ=255)
                                                                                       cACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtctc       <  1:93996/100‑1 (MQ=255)
                                                                                       cACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtctc       >  2:1301289/1‑100 (MQ=255)
                                                                                        aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa      <  2:2944426/100‑1 (MQ=255)
                                                                                           cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa   >  2:3354411/1‑100 (MQ=255)
                                                                                            acaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAt  <  1:1080648/100‑1 (MQ=255)
                                                                                                |                                                                                               
CCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAT  >  NC_000913/1579998‑1580189

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