Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,837,903 A→G K138R (AAA→AGA)  b1757 → putative thiosulfate sulfur transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,837,9030AG98.5% 439.0 / ‑2.9 138K138R (AAA→AGA) b1757putative thiosulfate sulfur transferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (70/65);  total (72/66)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

Only 100 of 147 total aligned reads displayed.

GACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT  >  NC_000913/1837815‑1837992
                                                                                        |                                                                                         
gACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTg                                                                                       <  2:2210022/93‑1 (MQ=255)
  cGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTgac                                                                                     >  1:219306/1‑91 (MQ=255)
  cGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGcc                                                                                     <  1:759706/93‑1 (MQ=255)
   gCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCg                                                                                    >  2:4113614/1‑93 (MQ=255)
     cGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCa                                                                                  >  2:1217985/1‑93 (MQ=255)
     cGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCa                                                                                  <  2:1580423/93‑1 (MQ=255)
      gTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAt                                                                                 >  1:1124746/1‑93 (MQ=255)
       tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                >  2:3114000/1‑93 (MQ=255)
       tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                >  1:4737242/1‑93 (MQ=255)
       tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                <  2:4552723/93‑1 (MQ=255)
        cAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAttt                                                                               <  2:4232139/93‑1 (MQ=255)
         aaaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtttt                                                                              >  2:1607918/1‑93 (MQ=255)
         aaaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtttt                                                                              <  1:4673037/93‑1 (MQ=255)
         aaaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtttt                                                                              >  1:371265/1‑93 (MQ=255)
          aaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTg                                                                             >  2:3389535/1‑93 (MQ=255)
          aaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTg                                                                             >  1:3865767/1‑93 (MQ=255)
          aaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCATCTGCAAAGACTGCCGCATTTTg                                                                             >  2:3389536/1‑93 (MQ=255)
           aaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGa                                                                            <  1:3611523/93‑1 (MQ=255)
              gcgACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagca                                                                         >  2:1973455/1‑93 (MQ=255)
                gACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagcagc                                                                       <  2:1152873/93‑1 (MQ=255)
                 aCTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCt                                                                      <  2:684320/93‑1 (MQ=255)
                  cTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTg                                                                     <  2:3700737/93‑1 (MQ=255)
                     cAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGtt                                                                  <  1:2774074/93‑1 (MQ=255)
                     cAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGtt                                                                  <  2:4566975/93‑1 (MQ=255)
                         aaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAATACTGCCGCATTTTGAGCAGCTGGTTTATc                                                              <  2:3470132/93‑1 (MQ=255)
                         aaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATc                                                              <  2:3329470/93‑1 (MQ=255)
                         aaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATc                                                              <  2:3329469/93‑1 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAGCCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             <  2:1496330/93‑1 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             >  2:4475384/1‑93 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             >  2:3444779/1‑93 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             <  2:1294059/93‑1 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATcc                                                             >  2:2464688/1‑93 (MQ=255)
                          aGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGCTTATcc                                                             >  2:387428/1‑93 (MQ=255)
                           gCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATccc                                                            >  2:691366/1‑92 (MQ=255)
                            cAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTCATCAGCTGGTTTATCCGc                                                           >  1:860704/1‑93 (MQ=255)
                                tttAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATg                                                       >  1:3574473/1‑93 (MQ=255)
                                tttAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATg                                                       <  2:3097490/93‑1 (MQ=255)
                                  tAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTTCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGc                                                     <  1:1880411/93‑1 (MQ=255)
                                  tAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGc                                                     <  1:2101995/93‑1 (MQ=255)
                                  tAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGc                                                     <  2:4273880/93‑1 (MQ=255)
                                   aaCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCt                                                    <  1:3343852/93‑1 (MQ=255)
                                    aCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGCCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTg                                                   >  1:55395/1‑93 (MQ=255)
                                     cGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGc                                                  <  1:3195575/93‑1 (MQ=255)
                                      gCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCa                                                 >  1:670776/1‑93 (MQ=255)
                                       cATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCAc                                                <  2:4669001/93‑1 (MQ=255)
                                        atatCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACg                                               >  1:1765839/1‑93 (MQ=255)
                                         tatCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGa                                              >  1:2001381/1‑93 (MQ=255)
                                           tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCATCTGGTTTATCCGCAATGGCTGCACGAcc                                            <  1:1280176/93‑1 (MQ=255)
                                           tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            >  1:2487699/1‑93 (MQ=255)
                                           tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            >  2:2562956/1‑93 (MQ=255)
                                           tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATCGCTGCACGAcc                                            <  2:871068/93‑1 (MQ=255)
                                            ctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCt                                           <  2:2891034/93‑1 (MQ=255)
                                            ctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCt                                           <  2:1899882/93‑1 (MQ=255)
                                             tcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTg                                          >  2:1194795/1‑93 (MQ=255)
                                             tcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTg                                          <  1:3076384/93‑1 (MQ=255)
                                              ccATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGc                                         <  1:2911887/93‑1 (MQ=255)
                                              ccATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGc                                         >  1:450789/1‑93 (MQ=255)
                                                aTCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       >  1:399185/1‑93 (MQ=255)
                                                aTCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       >  1:2855586/1‑93 (MQ=255)
                                                aTCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       >  1:4518907/1‑93 (MQ=255)
                                                   cTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaacaa                                    >  2:3942064/1‑93 (MQ=255)
                                                   cTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaacaa                                    >  1:2525975/1‑93 (MQ=255)
                                                    tGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAg                                   >  1:4169523/1‑93 (MQ=255)
                                                    tGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAg                                   >  1:1770557/1‑93 (MQ=255)
                                                     gAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAgg                                  >  1:2485609/1‑93 (MQ=255)
                                                     gAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAgg                                  >  2:644122/1‑93 (MQ=255)
                                                      aGTGACGCGCTAAGCGACCCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGt                                 <  2:3857845/93‑1 (MQ=255)
                                                       gTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTa                                <  1:655253/93‑1 (MQ=255)
                                                         gACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTaaa                              >  2:343349/1‑93 (MQ=255)
                                                           cgcgCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAaga                            <  2:860938/93‑1 (MQ=255)
                                                              gcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGAAACAAGGGAAAGAGGt                         >  2:271380/1‑93 (MQ=255)
                                                               cTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGtt                        <  1:2469464/93‑1 (MQ=255)
                                                                tAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTa                       >  2:1080312/1‑93 (MQ=255)
                                                                tAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTa                       >  2:2075205/1‑93 (MQ=255)
                                                                tAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTa                       >  2:4103979/1‑93 (MQ=255)
                                                                 aaGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTAc                      <  1:1213865/93‑1 (MQ=255)
                                                                 aaGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTAc                      <  2:1042288/93‑1 (MQ=255)
                                                                  aGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACg                     <  1:3135208/93‑1 (MQ=255)
                                                                    cGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGc                   <  1:1188207/93‑1 (MQ=255)
                                                                      aaCCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGa                 >  1:2344513/1‑93 (MQ=255)
                                                                          ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  2:2961034/93‑1 (MQ=255)
                                                                          ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  1:2219428/93‑1 (MQ=255)
                                                                          ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  1:870963/93‑1 (MQ=255)
                                                                          ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             >  1:2314153/1‑93 (MQ=255)
                                                                          ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  1:3238789/93‑1 (MQ=255)
                                                                             ccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGc          <  1:3162196/93‑1 (MQ=255)
                                                                             ccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGc          <  1:4181688/93‑1 (MQ=255)
                                                                             ccGTCGGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGc          <  2:3508977/93‑1 (MQ=255)
                                                                              cGTCTGCACAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGcc         <  1:3519086/93‑1 (MQ=255)
                                                                              cGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGGTACGGCGAAAACTGcc         >  2:1478219/1‑93 (MQ=255)
                                                                               gTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCg        >  1:339225/1‑93 (MQ=255)
                                                                                tCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCgg       >  1:3315674/1‑93 (MQ=255)
                                                                                tCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCgg       >  2:44448/1‑93 (MQ=255)
                                                                                tCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCgg       <  2:2085097/93‑1 (MQ=255)
                                                                                 cTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGt      <  1:1506961/93‑1 (MQ=255)
                                                                                  tGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTg     <  2:1497450/93‑1 (MQ=255)
                                                                                    cAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGAc   <  1:1288195/93‑1 (MQ=255)
                                                                                     aaagaCTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACt  <  1:3028321/93‑1 (MQ=255)
                                                                                     aaagaCTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACt  >  2:227592/1‑93 (MQ=255)
                                                                                     aaagaCCGCCGCATGTGGAGCAGCGGGGTTATCCGCAAAGGCTGCACGACCCGCAACAAGGGAAAGAGGTTACGGCGAAACCGGCCGGCGAcg  >  2:3616432/1‑92 (MQ=255)
                                                                                        |                                                                                         
GACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT  >  NC_000913/1837815‑1837992

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: ‑