Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,050,464 G→A S282S (TCG→TCA)  hemN → O2‑independent coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,050,4640GA85.5% 893.3 / 106.1 311S282S (TCG→TCA) hemNO2‑independent coproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base G (20/25);  new base A (149/117);  total (169/142)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

CTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTCCA  >  NC_000913/4050373‑4050556
                                                                                           |                                                                                            
ctgctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCgg                                                                                             >  2:6466707/1‑93 (MQ=255)
 tgctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCagg                                                                                            >  1:5354975/1‑93 (MQ=255)
  gctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGc                                                                                           >  1:1099272/1‑93 (MQ=255)
  gctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGc                                                                                           <  1:4227336/93‑1 (MQ=255)
  cctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGc                                                                                           <  1:2325146/92‑1 (MQ=255)
      aGCGCAAAATCAAAGATTCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATc                                                                                       <  1:3658936/93‑1 (MQ=255)
      aGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATc                                                                                       <  1:2643383/93‑1 (MQ=255)
       gcgcAAATTCAAAGATGCGGCCCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCa                                                                                      <  1:776111/93‑1 (MQ=255)
       gcgcAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCa                                                                                      <  1:1258681/93‑1 (MQ=255)
       gcgcAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCa                                                                                      <  1:390055/93‑1 (MQ=255)
        cgaaaaaTCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAg                                                                                     <  1:4916713/90‑1 (MQ=255)
            aaaTCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTa                                                                                 >  1:459185/1‑93 (MQ=255)
            aaaTCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTa                                                                                 <  1:5781051/93‑1 (MQ=255)
              aTCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATc                                                                               >  2:2860431/1‑93 (MQ=255)
              aTCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATc                                                                               <  1:2144156/93‑1 (MQ=255)
              aTCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATc                                                                               <  2:6174430/93‑1 (MQ=255)
               tCAAAGATTCTCACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCg                                                                              <  2:2993744/93‑1 (MQ=255)
                 aaaGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGt                                                                            >  2:265695/1‑93 (MQ=255)
                  aaGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTa                                                                           >  2:4825619/1‑93 (MQ=255)
                  aaGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTa                                                                           >  2:5047788/1‑93 (MQ=255)
                   aGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTAt                                                                          >  2:1862737/1‑93 (MQ=255)
                   aGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTAt                                                                          <  2:5331130/93‑1 (MQ=255)
                   aGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTAt                                                                          >  2:5891080/1‑93 (MQ=255)
                   aGATCCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTAt                                                                          <  1:545686/93‑1 (MQ=255)
                     aTGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATgg                                                                        <  1:2051933/93‑1 (MQ=255)
                      tGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGa                                                                       <  2:1351794/93‑1 (MQ=255)
                        gtGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTGCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATc                                                                     <  1:3198322/92‑1 (MQ=255)
                           aCCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACt                                                                  <  2:4171141/93‑1 (MQ=255)
                            ccTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACtt                                                                 <  1:4431782/93‑1 (MQ=255)
                             cTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACttt                                                                <  1:3847028/93‑1 (MQ=255)
                              tGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTg                                                               >  1:3033025/1‑93 (MQ=255)
                              tGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCAACGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTg                                                               >  2:3569318/1‑93 (MQ=255)
                               gCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGCAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGc                                                              >  2:1258995/1‑93 (MQ=255)
                               gCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGc                                                              >  2:1900308/1‑93 (MQ=255)
                                ccGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTGTCAGTTTATCGGTATGGATCACTTTGcc                                                             >  2:728984/1‑93 (MQ=255)
                                ccGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGcc                                                             <  1:552268/93‑1 (MQ=255)
                                 cGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGccc                                                            >  2:2821336/1‑93 (MQ=255)
                                    gTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTc                                                         >  1:1629929/1‑93 (MQ=255)
                                     tCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                                        <  2:4727540/93‑1 (MQ=255)
                                     tCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                                        >  2:4340195/1‑93 (MQ=255)
                                     tCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                                        <  2:4561491/93‑1 (MQ=255)
                                     tCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                                        <  2:5954908/93‑1 (MQ=255)
                                     tCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                                        >  1:430185/1‑93 (MQ=255)
                                      ccGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCCCTTTGCCCGTCCg                                                       <  1:4723961/93‑1 (MQ=255)
                                      ccGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCg                                                       >  1:1135131/1‑93 (MQ=255)
                                       cGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCgg                                                      <  2:712802/93‑1 (MQ=255)
                                        gcagcaAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGa                                                     <  2:233507/93‑1 (MQ=255)
                                           gcaAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGa                                                  >  2:2092460/1‑93 (MQ=255)
                                              aaaaCTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGa                                               >  1:68411/1‑93 (MQ=255)
                                              aaaaCTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGa                                               >  2:3401760/1‑93 (MQ=255)
                                               aaaCTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAg                                              <  2:3881810/93‑1 (MQ=255)
                                                aaCTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGc                                             <  2:6087212/93‑1 (MQ=255)
                                                 aCTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCt                                            <  1:590696/93‑1 (MQ=255)
                                                 aCTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCt                                            <  1:5298849/93‑1 (MQ=255)
                                                   tCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTgg                                          >  1:5071903/1‑93 (MQ=255)
                                                       tatCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGt                                      >  1:2781347/1‑93 (MQ=255)
                                                        atCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTg                                     <  1:2278735/93‑1 (MQ=255)
                                                        atCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTg                                     >  2:227307/1‑93 (MQ=255)
                                                         tCCTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTgg                                    >  2:2229128/1‑93 (MQ=255)
                                                          ccTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGc                                   >  1:5732445/1‑93 (MQ=255)
                                                           cTGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGcc                                  >  1:1930740/1‑93 (MQ=255)
                                                            tGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGccc                                 <  1:3703182/93‑1 (MQ=255)
                                                            tGCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGccc                                 >  2:1556739/1‑93 (MQ=255)
                                                             gCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCa                                >  1:5220184/1‑93 (MQ=255)
                                                             gCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCa                                >  2:1750027/1‑93 (MQ=255)
                                                             gCAGGAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCa                                >  1:1477109/1‑93 (MQ=255)
                                                              cAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAg                               <  2:5724773/93‑1 (MQ=255)
                                                                 gAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGt                            >  2:5148447/1‑93 (MQ=255)
                                                                 gAAACCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGt                            >  2:4247180/1‑93 (MQ=255)
                                                                    aCCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCGGGCGGTGGCCCAGCGTGaa                         <  1:6461850/93‑1 (MQ=255)
                                                                    aCCATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGaa                         >  2:3717589/1‑93 (MQ=255)
                                                                      cATCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAgg                       >  2:3313185/1‑93 (MQ=255)
                                                                       aTCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGc                      >  2:5524490/1‑93 (MQ=255)
                                                                        tCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCg                     <  2:1720316/93‑1 (MQ=255)
                                                                        tCGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCg                     >  1:2887003/1‑93 (MQ=255)
                                                                         cGCCTTCCTGACGCAATCATGCTATGAGTTTTTCGTTATGGATCACTTTGCCCGTCCGGATGACTAGCTGGCGGTGGCCCAGCGTGAAGGCGt                    <  1:6302414/93‑1 (MQ=255)
                                                                         cGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGt                    >  2:6062826/1‑93 (MQ=255)
                                                                         cGCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGt                    >  1:1744205/1‑93 (MQ=255)
                                                                          gCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtg                   <  1:2223783/93‑1 (MQ=255)
                                                                          gCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtg                   >  1:5600144/1‑93 (MQ=255)
                                                                          gCCTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtg                   <  2:3222814/93‑1 (MQ=255)
                                                                           ccTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtgc                  >  1:311446/1‑93 (MQ=255)
                                                                           ccTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtgc                  >  2:663720/1‑93 (MQ=255)
                                                                            cTTCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtgct                 >  1:1412326/1‑93 (MQ=255)
                                                                             ttCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtgctg                >  2:2549260/1‑93 (MQ=255)
                                                                             ttCCTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtgctg                >  2:1544389/1‑93 (MQ=255)
                                                                                cTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCTGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCAt             >  2:888418/1‑93 (MQ=255)
                                                                                cTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGCAGGCGTGCTGCAt             >  1:6585419/1‑93 (MQ=255)
                                                                                cTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCAt             <  2:6283863/93‑1 (MQ=255)
                                                                                cTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCAt             >  1:3396547/1‑93 (MQ=255)
                                                                                cTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCAt             <  2:6514543/93‑1 (MQ=255)
                                                                                cTGACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCAt             >  1:1074013/1‑93 (MQ=255)
                                                                                  gACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCg           <  2:5745555/93‑1 (MQ=255)
                                                                                  gACGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCg           >  2:4022736/1‑93 (MQ=255)
                                                                                   aCGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGt          >  1:6093263/1‑93 (MQ=255)
                                                                                   aCGCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGt          >  1:1702308/1‑93 (MQ=255)
                                                                                     gCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTaa        <  2:5242653/93‑1 (MQ=255)
                                                                                     gCAATCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTaa        <  1:39709/93‑1 (MQ=255)
                                                                                       aaTCGGGCTATCAGTTTATCGGTATGCATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACt      <  2:3338174/93‑1 (MQ=255)
                                                                                        aTCAGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACtt     >  1:3149347/1‑93 (MQ=255)
                                                                                           |                                                                                            
CTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTCCA  >  NC_000913/4050373‑4050556

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