Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,458,650 +G intergenic (‑161/+20) b2342 ← / ← b2343 putative acyltransferase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,458,6461.G99.0% 308.9 / ‑3.6 97intergenic (‑157/+24)b2342/b2343putative acyltransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (49/47);  total (50/47)
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.86e-01

Only 100 of 107 total aligned reads displayed.

ACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCC  >  NC_000913/2458554‑2458733
                                                                                             |                                                                                       
aCATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACa                                                                                          >  1:3111282/1‑93 (MQ=255)
aCATTTAGCCAAGGGGAGAAAAGGGAAAGGGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACa                                                                                          <  1:1672548/93‑1 (MQ=255)
 cATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAg                                                                                        >  2:2802591/1‑93 (MQ=255)
 cATTTAGCCAAGTGGAGAAAAGGGAAAGTGGAAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAg                                                                                        <  2:4291839/93‑1 (MQ=255)
  aTTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAgg                                                                                       <  1:464940/93‑1 (MQ=255)
   tttAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAggg                                                                                      <  1:3127469/93‑1 (MQ=255)
   ttgcgCCAGGTGGAAAAAGGGGAGCGGGGGAGCTAATACAGAGAGAGAAAAAAAGAAAAAAAGAGAACACGCCCCTGCCAGCAACTGACAggg                                                                                      <  2:4173941/89‑1 (MQ=255)
      aGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAgggggc                                                                                   >  2:4391179/1‑91 (MQ=255)
          aaGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGcgtac                                                                                >  1:4364526/1‑93 (MQ=255)
          aaGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGcgtac                                                                                >  1:620428/1‑93 (MQ=255)
          aaGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGcgtac                                                                                >  2:194717/1‑93 (MQ=255)
          aaGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGcgtac                                                                                >  2:1892193/1‑93 (MQ=255)
           aGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGcgtacg                                                                               >  2:2446599/1‑93 (MQ=255)
            gTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGcgtacgt                                                                              <  1:181275/93‑1 (MQ=255)
            gTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGcgtacgt                                                                              >  1:2517806/1‑93 (MQ=255)
            gTGAAGAAAGGGAAAAGGGGGAGCTATAACAAAGAAAAAAAAAAAAAAAAAAAAAAAAAACGCCCCTGCCAGCAACTGACAGGGGGcgtacgt                                                                              <  1:863243/93‑1 (MQ=255)
                agaAAAGAGAAAGAGAGAGCTATGACACAGAGAGAAAAGAAGAAGAAAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTAccc                                                                          <  2:1958471/93‑1 (MQ=255)
                 gaAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCg                                                                         <  1:2096129/93‑1 (MQ=255)
                  aaaaGGGAAAGTGGGAGCTATGACAAAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCga                                                                        <  1:2778976/93‑1 (MQ=255)
                    aaGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCgaga                                                                      >  1:2521285/1‑93 (MQ=255)
                     aGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGagag                                                                     <  2:3245649/93‑1 (MQ=255)
                      gggAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAgg                                                                    <  1:1602748/93‑1 (MQ=255)
                       ggAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGa                                                                   <  1:1603031/93‑1 (MQ=255)
                        gAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAGCGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGaa                                                                  <  1:483373/93‑1 (MQ=255)
                         aaaGTGGGAGCTATGCCACGGAGAGAAAGGAAGCAGAGAAAGGAAAACGCCCCTGCCAGCGACTGACAGGGGGCGTACGTACCCGAGCGGact                                                                 <  1:4224990/93‑3 (MQ=255)
                         aaaGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAAt                                                                 <  2:3763175/93‑1 (MQ=255)
                         aaaGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAAt                                                                 >  1:98074/1‑93 (MQ=255)
                         aaaGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAAt                                                                 >  1:846947/1‑93 (MQ=255)
                         aaaGTGGGAGCTATGACACAGAGAAAAAAGAAGAAGAGAAAAGAAAACTCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAAt                                                                 <  1:4092439/93‑1 (MQ=255)
                           aGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTa                                                               >  1:4430078/1‑93 (MQ=255)
                             tGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAAc                                                             <  1:1143836/93‑1 (MQ=255)
                              gggAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACg                                                            >  1:1217339/1‑93 (MQ=255)
                              gggAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACg                                                            >  2:3905524/1‑93 (MQ=255)
                                gAGCTATGACACAGAGAGAAAAGAAAAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCa                                                          <  2:4332583/93‑1 (MQ=255)
                                  gCTATGACACAGAGAGAAATGAAGAAGAGAAAAGAATACGCCCGTGCCAGCAACTGACAGGGGGCGTCCGTCCCCGAGAGGAATTAACGCAGa                                                        <  2:3292338/93‑1 (MQ=255)
                                  gCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGa                                                        <  1:3158620/93‑1 (MQ=255)
                                     aTGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGAccc                                                     <  2:4112819/93‑1 (MQ=255)
                                       gACACAGAGAGAAAAGAAGAAAAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAg                                                   <  1:3576472/93‑1 (MQ=255)
                                         cacaGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTTCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCt                                                 >  1:463570/1‑93 (MQ=255)
                                          acaGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTTCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTg                                                >  2:4211523/1‑93 (MQ=255)
                                          acaGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTg                                                <  1:2953287/93‑1 (MQ=255)
                                             gagagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaa                                             >  1:556054/1‑93 (MQ=255)
                                               gagaAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGaaaa                                           >  1:4246512/1‑93 (MQ=255)
                                                   aaaGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAg                                       >  2:3251357/1‑93 (MQ=255)
                                                    aagaagaagAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGc                                      >  2:1198224/1‑93 (MQ=255)
                                                     agaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCa                                     >  2:4170519/1‑93 (MQ=255)
                                                     agaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCa                                     >  2:1820666/1‑93 (MQ=255)
                                                      gaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAt                                    >  2:2034128/1‑93 (MQ=255)
                                                      gaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAtt                                   >  1:2343491/1‑93 (MQ=255)
                                                      aaagaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAt                                    <  1:2045127/92‑1 (MQ=255)
                                                        agaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCAttt                                  >  2:1355093/1‑93 (MQ=255)
                                                         gaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTc                                 <  1:4697878/93‑1 (MQ=255)
                                                         gaagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTc                                 >  2:1879324/1‑93 (MQ=255)
                                                          aagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCa                                <  2:3164620/93‑1 (MQ=255)
                                                          aagaGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCa                                <  1:1931453/93‑1 (MQ=255)
                                                           agagAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAg                               >  1:1249386/1‑93 (MQ=255)
                                                            gagaAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGc                              >  2:1605891/1‑93 (MQ=255)
                                                            gagaAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGc                              <  1:3994913/93‑1 (MQ=255)
                                                             agaaaagaaaACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCt                             <  2:1195121/93‑1 (MQ=255)
                                                             agaaaagaaaACGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCt                             >  2:303089/1‑93 (MQ=255)
                                                              gaaaagaaaaCGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCATCtt                            <  2:4646822/93‑1 (MQ=255)
                                                              aaaaagaaaaCGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCtt                            <  1:966759/92‑1 (MQ=255)
                                                                aaagaaaaCGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCg                          >  1:4151779/1‑93 (MQ=255)
                                                                aaagaaaaCGCCCCTGCCAGCAACTGACAGGGGGCGTACATACCCGAGAGGAATTAACGCAGACCCAGCTGGCAAATCAGCATTTCAGCTTCg                          >  1:4235074/1‑93 (MQ=255)
                                                                 aagaaaaCGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGc                         <  1:4544679/93‑1 (MQ=255)
                                                                  agaaaaCGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCa                        <  2:3111282/93‑1 (MQ=255)
                                                                   gaaaaCGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAg                       <  2:2465993/93‑1 (MQ=255)
                                                                    aaacCGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAgg                      <  1:2312266/93‑1 (MQ=255)
                                                                     aatcGCCCGTGCCAGCAGCTGGGAGGGGGTGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGc                     <  1:3727970/90‑1 (MQ=255)
                                                                     aaccGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGc                     <  1:3727969/90‑1 (MQ=255)
                                                                       aCGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCaa                   >  2:4501044/1‑93 (MQ=255)
                                                                        cGCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAAATAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCaaa                  <  1:4235602/93‑1 (MQ=255)
                                                                        cGCCCCTGCCAGCAACTCACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCaaa                  >  2:4235887/1‑93 (MQ=255)
                                                                         gCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCaaaa                 >  1:441941/1‑93 (MQ=255)
                                                                         gCCCCTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCaaaa                 <  1:2179917/93‑1 (MQ=255)
                                                                           cccTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGt               <  1:4131202/93‑1 (MQ=255)
                                                                           cccTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGt               <  1:3565113/93‑1 (MQ=255)
                                                                           cccTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGt               >  2:1261333/1‑93 (MQ=255)
                                                                            ccTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGta              >  2:3650085/1‑92 (MQ=255)
                                                                            ccTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTg              >  2:2657578/1‑93 (MQ=255)
                                                                             cTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGa             <  1:2315181/93‑1 (MQ=255)
                                                                             cTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGa             >  2:4404689/1‑93 (MQ=255)
                                                                             cTGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGa             <  1:3944364/93‑1 (MQ=255)
                                                                              tGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGaa            <  2:2259173/93‑1 (MQ=255)
                                                                              tGCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGaa            >  1:408768/1‑93 (MQ=255)
                                                                               gCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGaaa           <  1:243335/93‑1 (MQ=255)
                                                                               gCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGaaa           <  2:1782544/93‑1 (MQ=255)
                                                                               gCCAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGaaa           <  1:618886/93‑1 (MQ=255)
                                                                                ccAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAAt          <  2:2864114/93‑1 (MQ=255)
                                                                                 cAGCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATc         >  2:1005001/1‑93 (MQ=255)
                                                                                   gCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATcta       >  1:3104972/1‑91 (MQ=255)
                                                                                   gCAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGa       >  2:3692629/1‑93 (MQ=255)
                                                                                    cAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGat      <  2:4199211/93‑1 (MQ=255)
                                                                                    cAACTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGat      <  1:830458/93‑1 (MQ=255)
                                                                                     aaCTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGata     >  2:3883707/1‑93 (MQ=255)
                                                                                      aCTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGatat    >  2:1586900/1‑93 (MQ=255)
                                                                                       cTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATc   >  2:1851921/1‑93 (MQ=255)
                                                                                       cTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATc   >  2:1851920/1‑93 (MQ=255)
                                                                                       cTGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATc   >  1:1461287/1‑93 (MQ=255)
                                                                                        tGACAGGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATcc  <  1:3226690/93‑1 (MQ=255)
                                                                                             |                                                                                       
ACATTTAGCCAAGTGGAGAAAAGGGAAAGTGGGAGCTATGACACAGAGAGAAAAGAAGAAGAGAAAAGAAAACGCCCCTGCCAGCAACTGACAGGGGCGTACGTACCCGAGAGGAATTAACGCAGACCCAGCTGGAAAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCC  >  NC_000913/2458554‑2458733

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