Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 649,964:1 +G 100% coding (527/909 nt) citE ← citrate lyase, citryl‑ACP lyase (beta) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913649,9601.G65.8% 260.8 / 155.6 114coding (531/909 nt)citEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base . (23/16);  new base G (37/37);  total (60/54)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

GTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTA  >  NC_000913/649861‑650059
                                                                                                    |                                                                                                   
gttgttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGttccttcc                                                                                                      >  2:272524/1‑100 (MQ=255)
gttgttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGttccttcc                                                                                                      >  2:3055316/1‑100 (MQ=255)
  tgttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑gg                                                                                                   <  2:1076555/100‑1 (MQ=255)
   gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑ggg                                                                                                  <  2:119495/100‑1 (MQ=255)
   gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑ggg                                                                                                  <  2:2448999/100‑1 (MQ=255)
   gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑ggg                                                                                                  <  2:2610476/100‑1 (MQ=255)
    ttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑gggg                                                                                                 <  2:757516/100‑1 (MQ=255)
    ttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTGCAGTTCCTTCC‑gggg                                                                                                 >  1:3144529/1‑100 (MQ=255)
     tAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑ggggg                                                                                                <  2:3343586/100‑2 (MQ=255)
     tAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGa                                                                                                >  1:3295023/1‑100 (MQ=255)
      aGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAg                                                                                               >  1:1730368/1‑100 (MQ=255)
       gCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑gggggag                                                                                              <  2:32384/100‑4 (MQ=255)
        cGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCg                                                                                             >  1:3345990/1‑100 (MQ=255)
          tCGGAATAGAGGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑gggggagcg                                                                                            <  1:2883099/99‑6 (MQ=255)
           cGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGTGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGtt                                                                                           >  2:2599648/1‑100 (MQ=255)
            ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGTGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCAGGGGAGCGTTc                                                                                          <  2:2831474/100‑1 (MQ=255)
            ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAACGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCt                                                                                         >  2:2395082/1‑100 (MQ=255)
              aaTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTg                                                                                        >  1:1750829/1‑100 (MQ=255)
               aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTa                                                                                      <  1:2726597/100‑1 (MQ=255)
               aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCGGCAAAATGGAGCAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGt                                                                                       <  1:2109659/100‑1 (MQ=255)
                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCTGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACg                                                                                    >  1:876562/1‑100 (MQ=255)
                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACg                                                                                    >  2:1408854/1‑100 (MQ=255)
                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGc                                                                                   >  2:2999140/1‑100 (MQ=255)
                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGc                                                                                   >  1:2805047/1‑100 (MQ=255)
                   aCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGc                                                                                   >  1:2881055/1‑100 (MQ=255)
                   aCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGc                                                                                   >  1:2882026/1‑100 (MQ=255)
                     ggTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAgg                                                                                >  1:1657102/1‑100 (MQ=255)
                           gAACGCCTGAATACCCGCAGAGTGCGCGGCCTGCAAAATGGAACCCCGTGCGCACAGCAGCTCAGTTCCTTCA‑GGGGAGCGTTCTGTACGCAGGTTgcgc                                                                          >  1:823986/1‑100 (MQ=255)
                            aaCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGca                                                                         <  2:1720479/100‑1 (MQ=255)
                            aaCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCGTCCGGGGGAGCGTTCTGTACGCAGGTTgcgc                                                                          <  1:1054523/100‑1 (MQ=255)
                             aCGCCTGAATCCCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGca                                                                         <  1:973194/100‑1 (MQ=255)
                             aCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGca                                                                         <  2:1195788/100‑1 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGcaca                                                                       >  1:1294685/1‑100 (MQ=255)
                                ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACAt                                                                      >  2:1927149/1‑100 (MQ=255)
                                ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACAt                                                                      >  1:1637411/1‑100 (MQ=255)
                                ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATa                                                                     >  2:3177051/1‑100 (MQ=255)
                                  tGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAg                                                                    >  2:1368806/1‑100 (MQ=255)
                                   gAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGt                                                                   <  2:2861946/100‑1 (MQ=255)
                                   gAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtc                                                                  <  1:466027/100‑1 (MQ=255)
                                     aTACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtctt                                                                >  2:481518/1‑100 (MQ=255)
                                       aCCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                               <  1:1437342/100‑1 (MQ=255)
                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                              >  2:317476/1‑100 (MQ=255)
                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                              >  2:975665/1‑100 (MQ=255)
                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                              <  1:2671974/100‑1 (MQ=255)
                                         ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGAGAGCGTTCTGTACGCAGGTTGCGCACATAGTCCTCTGc                                                            >  2:1982374/1‑100 (MQ=255)
                                          cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                            >  2:1547454/1‑100 (MQ=255)
                                           gCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCCGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTCCTGCa                                                           >  2:124630/1‑100 (MQ=255)
                                           gCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCa                                                           <  1:963341/100‑1 (MQ=255)
                                           gCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTATACGCAGGTTGCGCACATAGTCTTCTGCa                                                           >  2:2772581/1‑100 (MQ=255)
                                           gCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCAc                                                          <  2:1657109/100‑1 (MQ=255)
                                             agagCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCAcc                                                         <  2:2155103/100‑1 (MQ=255)
                                                gcgcgcGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGa                                                       <  1:471228/99‑1 (MQ=255)
                                                gcgcgcGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAg                                                      <  2:872334/100‑1 (MQ=255)
                                                 cgcgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTGAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAgg                                                     <  1:1626262/100‑1 (MQ=255)
                                                 cgcgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGAAGGTTGCGCACATAGTCTTCTGCACCGAggg                                                    >  1:2644502/1‑100 (MQ=255)
                                                  gcgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAggg                                                    >  2:1134779/1‑100 (MQ=255)
                                                  gcgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAggg                                                    <  1:2752853/100‑1 (MQ=255)
                                                  gcgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAggg                                                    >  2:1992238/1‑100 (MQ=255)
                                                   cgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGc                                                   >  2:2138958/1‑100 (MQ=255)
                                                   cgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGc                                                   <  1:1186015/100‑1 (MQ=255)
                                                   cgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGc                                                   <  2:2881078/100‑1 (MQ=255)
                                                       gCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATa                                               >  1:1753931/1‑100 (MQ=255)
                                                       gCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAc                                              <  2:379154/100‑1 (MQ=255)
                                                        ccTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                             <  2:2629034/100‑1 (MQ=255)
                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                             >  2:3011961/1‑100 (MQ=255)
                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                             <  1:1502452/100‑1 (MQ=255)
                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                             <  1:2593409/100‑1 (MQ=255)
                                                           gCAAAATGGGACAGCGTGCCGACAGTAGTTCAGCTCCGTCCGGGGGAGCCTGCTGTACGCTGGTTGCGCACATAGTCTTCGGCACCGAggg                                                    <  1:2254527/91‑1 (MQ=255)
                                                           gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                           >  2:364157/1‑100 (MQ=255)
                                                           gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                           >  2:463483/1‑100 (MQ=255)
                                                            cAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGAt                                          <  2:524419/100‑1 (MQ=255)
                                                              aaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                        <  2:956818/100‑1 (MQ=255)
                                                              aaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaa                                       <  1:2082995/100‑1 (MQ=255)
                                                               aaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaa                                       <  2:1556764/100‑1 (MQ=255)
                                                               aaTGGAACAGCGGGCGAACAGCAGGTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaaa                                      <  1:493243/100‑1 (MQ=255)
                                                                aTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaaa                                      <  2:340410/100‑1 (MQ=255)
                                                                aTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCCATACCCACCAAAc                                     >  1:3155719/1‑100 (MQ=255)
                                                                       aGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGa                              >  1:323319/1‑100 (MQ=255)
                                                                        gCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGa                              <  2:1044985/100‑1 (MQ=255)
                                                                         cGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGaa                             <  1:2457784/100‑1 (MQ=255)
                                                                          gTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAg                            <  2:1287544/100‑1 (MQ=255)
                                                                          gTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGc                           >  2:3293745/1‑100 (MQ=255)
                                                                           tGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGCAGc                           >  2:3171162/1‑100 (MQ=255)
                                                                           tGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGc                           <  2:2191763/100‑1 (MQ=255)
                                                                            gCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                          <  2:3075669/100‑1 (MQ=255)
                                                                             cGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGt                         >  1:2130439/1‑100 (MQ=255)
                                                                                 cagcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGc                     >  1:198458/1‑100 (MQ=255)
                                                                                   gcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGa                   >  2:2451626/1‑100 (MQ=255)
                                                                                   gcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGa                   >  2:680805/1‑100 (MQ=255)
                                                                                   gcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGa                   >  2:752226/1‑100 (MQ=255)
                                                                                      gttcagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAttt                >  1:1260771/1‑100 (MQ=255)
                                                                                      gttcagttcCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTc               <  2:1730379/100‑1 (MQ=255)
                                                                                      gttcagttcCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTc               >  2:1816104/1‑100 (MQ=255)
                                                                                       ttcagttcCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTcc              <  2:3113161/100‑1 (MQ=255)
                                                                                         cagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCa             >  1:1419782/1‑100 (MQ=255)
                                                                                         cagttcCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCAc            <  1:627860/100‑1 (MQ=255)
                                                                                         cagttcCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCAc            >  1:982573/1‑100 (MQ=255)
                                                                                                 tcc‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCggg    <  1:3155453/100‑1 (MQ=255)
                                                                                                  cc‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCggg    >  2:3074104/1‑99 (MQ=255)
                                                                                                   c‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTa  >  2:2672879/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
GTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCGGGTA  >  NC_000913/649861‑650059

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