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.G89.6% 287.4 / 30.9 96coding (531/909 nt)citEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base . (6/4);  new base G (51/34);  total (58/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  NC_000913/649861‑650054
                                                                                                      |                                                                                                 
gttgttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGttccttcc                                                                                                    <  1:583835/100‑1 (MQ=255)
   gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggg                                                                                               <  1:1249179/100‑1 (MQ=255)
   gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggg                                                                                               <  2:40550/100‑1 (MQ=255)
     tAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggggg                                                                                            >  2:2900953/1‑99 (MQ=255)
     tAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggggg                                                                                            <  1:2744847/100‑2 (MQ=255)
      aGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggggga                                                                                           >  1:2842646/1‑98 (MQ=255)
       gCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACATCGTGCGAAAAGCAGTTCAGTTTCTTTTGGGGAGc                                                                                         >  1:2846397/1‑100 (MQ=255)
       gCGTCGGAATAGACGGTATCGAACGACTTAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGTACAGTAGTTCAGTTCCTTCCgggggag                                                                                          <  2:388026/100‑4 (MQ=255)
         gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggggagcg                                                                                        <  1:1523511/100‑6 (MQ=255)
            ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTc                                                                                     >  1:390492/1‑100 (MQ=255)
            ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTc                                                                                     >  1:2217696/1‑100 (MQ=255)
            ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTc                                                                                     <  1:2080509/100‑1 (MQ=255)
               aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCAAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGt                                                                                  >  1:2132552/1‑100 (MQ=255)
                 aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTAc                                                                                >  2:445443/1‑100 (MQ=255)
                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACg                                                                               >  1:70469/1‑100 (MQ=255)
                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCATGCGAACAGCAGCTCAGTTCCTTCCGGGGGAGCGTTCTGTACg                                                                               >  2:920459/1‑100 (MQ=255)
                   aCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGc                                                                              >  2:2141800/1‑100 (MQ=255)
                   aCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAACAGTTCAGTTCATTCCAGGAGAGCGTTCTGTACGc                                                                              >  1:647384/1‑100 (MQ=255)
                    cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGc                                                                              <  2:1867814/99‑1 (MQ=255)
                    cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCa                                                                             >  1:1289279/1‑100 (MQ=255)
                    cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCa                                                                             >  2:2903585/1‑100 (MQ=255)
                    cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCCAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCa                                                                             >  2:593695/1‑100 (MQ=255)
                     ggTATCGAACGTCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAg                                                                            <  1:1939187/100‑1 (MQ=255)
                       tATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGt                                                                          >  2:2848943/1‑100 (MQ=255)
                         tCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTg                                                                        >  1:1374365/1‑100 (MQ=255)
                          cGAACGCCTGAATACCCGCAAAGCGAGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGttac                                                                       >  2:1140062/1‑98 (MQ=255)
                           ggaCGCCTGAATACCCGCAGAGCGCGCGGCCTGTAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgcg                                                                      >  1:2058214/3‑100 (MQ=255)
                           gAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAATAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgcg                                                                      <  2:895651/100‑1 (MQ=255)
                           gAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgcg                                                                      >  1:2229706/1‑100 (MQ=255)
                            aaCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgcgc                                                                     >  2:2519442/1‑100 (MQ=255)
                             aCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGca                                                                    <  2:1401823/100‑1 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGTGGGAGCGTTCTGTACGCAGGTTGCGca                                                                    <  2:1012304/99‑1 (MQ=255)
                                ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACAt                                                                 >  2:766267/1‑100 (MQ=255)
                                ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGCAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGCGCGTTCTGTACGCAGGTTGCGCACAt                                                                 >  2:652401/1‑100 (MQ=255)
                                catGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAAGAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCTTTCTGTACGCAGGTTGCGCACAt                                                                 <  2:931591/98‑1 (MQ=255)
                                 cTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACATTGTGTTAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACAt                                                                 <  1:805333/99‑1 (MQ=255)
                                 cTGAATACCCGCAGAGCGCGCGGCCTGAAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATa                                                                >  2:2091350/1‑100 (MQ=255)
                                  tGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAg                                                               <  2:2229683/100‑1 (MQ=255)
                                   gAATACCCGTTGTGCGCGAGGCCTGCAAAAGGGAACAGCTTGCGAACAGCAGTTCAGTTCCGTCCGGGGGAGCGTTCTGTACGCAGGTTGCGca                                                                    <  1:1238241/94‑1 (MQ=255)
                                      tACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtct                                                            <  2:2821898/99‑1 (MQ=255)
                                       aCCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          >  2:1468215/1‑100 (MQ=255)
                                       aCCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          <  1:1016854/100‑1 (MQ=255)
                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCAGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  1:757457/1‑100 (MQ=255)
                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  2:1206117/1‑100 (MQ=255)
                                        cccGCAGAGCGCGCGACCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTGCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          >  2:1604471/1‑99 (MQ=255)
                                         ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          >  1:2030853/1‑98 (MQ=255)
                                         ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTAAGCAGGTTGCGCACATAGTCTTCTg                                                        >  1:2540435/1‑100 (MQ=255)
                                         ccGCAGAGCGCGCGGACTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCATTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  2:1096838/1‑99 (MQ=255)
                                          cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                       >  2:1450459/1‑100 (MQ=255)
                                          cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                       <  1:1450566/100‑1 (MQ=255)
                                           gCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCa                                                      >  2:2113546/1‑100 (MQ=255)
                                              gagCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCATGTTGTGCACATAGTCTTCAGCACCg                                                   <  2:1205683/100‑1 (MQ=255)
                                              gagCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGa                                                  >  2:740197/1‑100 (MQ=255)
                                                 cgcgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGGGCACATAGTCTTCTGCACCGAgg                                                <  2:2163862/100‑1 (MQ=255)
                                                    gcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCg                                             <  2:1657768/100‑1 (MQ=255)
                                                    gcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCg                                             >  2:1269168/1‑100 (MQ=255)
                                                     cgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGAACATAGTCTTCTGCACCGAGGGCGa                                            >  1:2598806/1‑100 (MQ=255)
                                                      ggCCTGCAAAATGTAACAGCGTGGGAACAGTAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCATGTTGCGCACATAGTCTTCTGCACCGAGGGc                                              <  2:2060758/97‑1 (MQ=255)
                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        <  2:1347036/100‑1 (MQ=255)
                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCACACATAGTCTTCTGCACCGAGGGCGATAcc                                        <  1:2422445/100‑1 (MQ=255)
                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGTAGCGTTCTGTACGAAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        >  1:103630/1‑100 (MQ=255)
                                                          tGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCg                                       <  2:2522851/100‑1 (MQ=255)
                                                           gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                      >  1:2580583/1‑100 (MQ=255)
                                                           gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                      <  1:127130/100‑1 (MQ=255)
                                                           gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGACTTCTGCACCGAGGGCGATACCGa                                      >  1:1659898/1‑100 (MQ=255)
                                                            cAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGATTGCAAACATAGTCTTCTGCACCGATGGCGATACCGAt                                     >  1:889355/1‑100 (MQ=255)
                                                             aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    <  2:298572/100‑1 (MQ=255)
                                                              aaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                   >  2:2021197/1‑100 (MQ=255)
                                                              aaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                   <  2:70469/100‑1 (MQ=255)
                                                              aaaTGGAACAGCGTGCGAACAGCAGTTAAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                   >  2:2325844/1‑100 (MQ=255)
                                                                 tGGAGCAGCGTGCGAACGGCAGTTCAGTTCCGTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaaa                                 <  2:2491958/99‑1 (MQ=255)
                                                                  tgAACAGCGTGCGAACAGCAGTTAAGTTCTTGCAGGGGGAGCGTTCTGTGTGCAGGTTGCGTTCATAGTCTTCTTCATCGAGGGCGTTACCGATCAAAc                                >  2:1988280/2‑99 (MQ=255)
                                                                  ggAACATCGTGCGAACAGCAGTTCAGTTCCTTCCGTGGGAGCGTTCTGTACGCAGGTTGCGAAGATAGTCTTCTGCACCGAGGGCGATACCGATCAAAc                                <  1:19692/99‑1 (MQ=255)
                                                                  ggAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACg                               >  2:1831480/1‑100 (MQ=255)
                                                                  ggAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCAATACCGATCAAACg                               <  2:1996979/100‑1 (MQ=255)
                                                                        gCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGa                         >  1:957461/1‑100 (MQ=255)
                                                                         cGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGaa                        >  1:478750/1‑100 (MQ=255)
                                                                           tGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGc                      >  1:2472760/1‑100 (MQ=255)
                                                                            gCGAACAGCAGTTCAGTTCCTTTCGGGGGAGTTTTTTGTACGCAGGTTGTGCACATAGTCTGCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     <  1:427463/100‑1 (MQ=255)
                                                                            gCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     >  2:16781/1‑100 (MQ=255)
                                                                            gCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     >  1:2214643/1‑100 (MQ=255)
                                                                             cGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGt                    >  1:1618089/1‑100 (MQ=255)
                                                                             cGAACAGCAGTTCAGTCGCTTCCGGGAGAGCGTTCTGTACGCAGATTGCGTACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGTAAGCGt                    <  2:836258/100‑1 (MQ=255)
                                                                              gAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTg                   <  1:2938920/100‑1 (MQ=255)
                                                                                 cagcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTGCGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATTCTGATCAATCGCTCGGAAGCGTGAg                 <  2:1412706/99‑1 (MQ=255)
                                                                                  agcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCg               <  1:2894858/100‑1 (MQ=255)
                                                                                   gcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGa              >  1:1936857/1‑100 (MQ=255)
                                                                                    cagttcagttCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAt             >  1:1831194/1‑100 (MQ=255)
                                                                                     agttcagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAtt            <  2:1877143/100‑1 (MQ=255)
                                                                                        tcagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTcc         <  2:195879/100‑1 (MQ=255)
                                                                                         cagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCGCATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCa        >  1:2908768/1‑100 (MQ=255)
                                                                                         cagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCa        <  2:2977929/100‑1 (MQ=255)
                                                                                          agttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCAc       >  1:1698515/1‑100 (MQ=255)
                                                                                          agttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCAc       <  1:993858/100‑1 (MQ=255)
                                                                                          agttcCTTCCGGGGGAGCGTTCTGGACGCAGGTTGCGCACATAGTGTTCTGCACAGAGGGGAATACCGATCGAAAGCTCGGAAGCGGGAGCGACTTcaa        >  2:2487105/1‑97 (MQ=255)
                                                                                           gttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACt      >  2:2407259/1‑100 (MQ=255)
                                                                                              ccttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcg   >  2:2887314/1‑100 (MQ=255)
                                                                                              ccttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcg   <  1:2783311/100‑1 (MQ=255)
                                                                                               cttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc  >  1:773137/1‑100 (MQ=255)
                                                                                               cttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGCAAGCGTGAGCGATTTCCACTgcgc  >  2:1304927/1‑100 (MQ=255)
                                                                                                      |                                                                                                 
GTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  NC_000913/649861‑650054

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