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,9030AG100.0% 655.0 / NA 203K138R (AAA→AGA) b1757putative thiosulfate sulfur transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (97/106);  total (97/106)

Only 100 of 208 total aligned reads displayed.

TCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT  >  NC_000913/1837813‑1837992
                                                                                          |                                                                                         
tCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGAc                                                                                         <  2:3550707/93‑1 (MQ=255)
tCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGAc                                                                                         <  2:2631218/93‑1 (MQ=255)
 cGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACt                                                                                        <  1:851787/93‑1 (MQ=255)
  gACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTg                                                                                       <  1:7441740/93‑1 (MQ=255)
   aCGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGc                                                                                      <  1:5516089/93‑1 (MQ=255)
   aCGCCGAAAAAACTCGACCTCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGGGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGc                                                                                      <  2:7835056/93‑1 (MQ=255)
      ccGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGc                                                                                   >  1:510719/1‑93 (MQ=255)
      ccGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGc                                                                                   >  1:3225355/1‑93 (MQ=255)
       cGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCa                                                                                  <  2:2420840/93‑1 (MQ=255)
        gTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAt                                                                                 >  2:210683/1‑93 (MQ=255)
        gTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAt                                                                                 >  1:2441096/1‑93 (MQ=255)
         tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                >  1:4013412/1‑93 (MQ=255)
         tCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                >  2:5251811/1‑93 (MQ=255)
         tCAAAACGCGACTGCAAAAAGCAGGTTCAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtt                                                                                <  2:6717035/93‑1 (MQ=255)
          cAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAttt                                                                               >  1:4278699/1‑93 (MQ=255)
           aaaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtttt                                                                              >  1:4906878/1‑93 (MQ=255)
           aaaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCAtttt                                                                              >  1:2077060/1‑93 (MQ=255)
            aaaCGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTg                                                                             <  2:6963051/93‑1 (MQ=255)
                gcgACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagca                                                                         >  2:597838/1‑93 (MQ=255)
                gcgACTGCAAAAAGCAGGTGTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagca                                                                         <  1:3154936/93‑1 (MQ=255)
                 cgACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagcag                                                                        >  1:7115454/1‑93 (MQ=255)
                 cgACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGagcag                                                                        <  2:7463181/93‑1 (MQ=255)
                   aCTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCt                                                                      <  2:7906341/93‑1 (MQ=255)
                    cTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTg                                                                     >  1:228296/1‑93 (MQ=255)
                     tGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTgg                                                                    <  1:1155441/93‑1 (MQ=255)
                       cAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGtt                                                                  <  2:7049311/93‑1 (MQ=255)
                       cAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGtt                                                                  <  1:2957241/93‑1 (MQ=255)
                        aaaaaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGttt                                                                 <  1:4388903/93‑1 (MQ=255)
                        aaaaaGCAGGTTTAACGCATATCTCCAGCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGttt                                                                 <  2:4592514/93‑1 (MQ=255)
                         aaaaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTa                                                                <  1:7903133/93‑1 (MQ=255)
                           aaGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATc                                                              <  2:676679/93‑1 (MQ=255)
                                ggTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCaa                                                         <  2:7768582/93‑1 (MQ=255)
                                 gTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAAt                                                        <  2:2305458/93‑1 (MQ=255)
                                 gTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAAt                                                        >  1:2285880/1‑93 (MQ=255)
                                 gTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAAt                                                        >  1:3676805/1‑93 (MQ=255)
                                  tttAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATg                                                       >  1:3918986/1‑93 (MQ=255)
                                  tttAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATg                                                       <  2:2028897/93‑1 (MQ=255)
                                    tAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGc                                                     <  2:7491331/93‑1 (MQ=255)
                                     aaCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCt                                                    >  1:5819864/1‑93 (MQ=255)
                                      aCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTg                                                   <  1:2352470/93‑1 (MQ=255)
                                      aCGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTg                                                   >  1:2668184/1‑93 (MQ=255)
                                       cGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGc                                                  <  1:2753885/93‑1 (MQ=255)
                                          atatCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACg                                               >  2:3328328/1‑93 (MQ=255)
                                          atatCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACg                                               >  2:4417055/1‑93 (MQ=255)
                                           tatCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGa                                              >  2:5395004/1‑93 (MQ=255)
                                            atCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAc                                             <  1:7022295/93‑1 (MQ=255)
                                            atCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAc                                             >  1:6756197/1‑93 (MQ=255)
                                            atCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAc                                             >  2:5763938/1‑93 (MQ=255)
                                            atCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAc                                             <  2:2974589/93‑1 (MQ=255)
                                            atCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTCAGCAGCTGGTTTATCCGCAATGGCTGCACGAc                                             >  1:3035873/1‑93 (MQ=255)
                                             tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            >  1:6941560/1‑93 (MQ=255)
                                             tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            >  1:4877055/1‑93 (MQ=255)
                                             tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            >  1:6065614/1‑93 (MQ=255)
                                             tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            <  2:3205086/93‑1 (MQ=255)
                                             tctcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGAcc                                            >  1:6781044/1‑93 (MQ=255)
                                               tcCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTg                                          <  2:5150748/93‑1 (MQ=255)
                                                ccATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGc                                         >  1:1382967/1‑93 (MQ=255)
                                                 cATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGca                                        <  1:811477/93‑1 (MQ=255)
                                                 cATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGca                                        >  2:6301071/1‑93 (MQ=255)
                                                  aTCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       >  2:1906085/1‑93 (MQ=255)
                                                  aTCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaa                                       >  1:1895762/1‑93 (MQ=255)
                                                   tCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaac                                      >  2:7520360/1‑93 (MQ=255)
                                                    ccTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaaca                                     <  1:5348147/93‑1 (MQ=255)
                                                    ccTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaaca                                     >  1:4721802/1‑93 (MQ=255)
                                                     cTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaacaa                                    >  2:6874146/1‑93 (MQ=255)
                                                     cTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGcaacaa                                    <  2:6534644/93‑1 (MQ=255)
                                                      tGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAg                                   >  1:1533980/1‑93 (MQ=255)
                                                        aGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGt                                 <  1:564638/93‑1 (MQ=255)
                                                           gACGCGCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTaaa                              <  1:3476577/93‑1 (MQ=255)
                                                             cgcgCTAAGCGAACCTTCGCGGCTGCAAAGACTGACGAATTTTGAGCAGATGGTTTATCCGAAATGGCTGAACGACCTGCAACAAGGTAAaga                            <  1:1261102/93‑1 (MQ=255)
                                                             cgcgCTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAaga                            <  1:3399222/93‑1 (MQ=255)
                                                              gcgcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAagag                           <  1:4299211/93‑1 (MQ=255)
                                                              gcgcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAagag                           <  1:5132903/93‑1 (MQ=255)
                                                               cgcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAgg                          <  2:3122520/93‑1 (MQ=255)
                                                               cgcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAgg                          <  1:4553116/93‑1 (MQ=255)
                                                                gcTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGt                         >  2:2861857/1‑93 (MQ=255)
                                                                 cTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGtt                        <  1:1601688/93‑1 (MQ=255)
                                                                 cTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGtt                        >  2:7661595/1‑93 (MQ=255)
                                                                 cTAAGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGtt                        >  2:4794325/1‑93 (MQ=255)
                                                                   aaGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTAc                      >  2:1512884/1‑93 (MQ=255)
                                                                   aaGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTAc                      <  1:3859258/93‑1 (MQ=255)
                                                                   aaGCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTAc                      <  1:4122302/93‑1 (MQ=255)
                                                                     gCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACgg                    <  1:3107736/93‑1 (MQ=255)
                                                                     gCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACgg                    >  1:6689164/1‑93 (MQ=255)
                                                                     gCGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGATACgg                    >  2:4830887/1‑93 (MQ=255)
                                                                      cGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGc                   <  1:5723829/93‑1 (MQ=255)
                                                                      cGAACCTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGc                   <  2:7224607/93‑1 (MQ=255)
                                                                           cTTCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAc              <  2:5355510/93‑1 (MQ=255)
                                                                            ttCCCGTCTGCAAAGAGTGCTGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  1:6401809/93‑1 (MQ=255)
                                                                            ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  1:7814728/93‑1 (MQ=255)
                                                                            ttCCCGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAAcc             <  2:683907/93‑1 (MQ=255)
                                                                              cccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGACACCTg           >  1:4169793/1‑93 (MQ=255)
                                                                              cccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTg           <  2:3895980/93‑1 (MQ=255)
                                                                              cccGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTg           <  2:1180106/93‑1 (MQ=255)
                                                                                cGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGcc         <  2:6078836/93‑1 (MQ=255)
                                                                                cGTCTGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGcc         <  2:5700955/93‑1 (MQ=255)
                                                                                    tGCAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTg     >  2:2192072/1‑93 (MQ=255)
                                                                                     tcAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGa    <  2:6444804/92‑1 (MQ=255)
                                                                                      cAAAGACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGAc   >  1:7183611/1‑93 (MQ=255)
                                                                                       aaagaCTCCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACt  <  2:1653292/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCGACGCCGTCAAAACGCGACTGCAAAAAGCAGGTTTAACGCATATCTCCATCCTGAGTGACGCGCTAAGCGAACCTTCCCGTCTGCAAAAACTGCCGCATTTTGAGCAGCTGGTTTATCCGCAATGGCTGCACGACCTGCAACAAGGTAAAGAGGTTACGGCGAAACCTGCCGGTGACT  >  NC_000913/1837813‑1837992

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