Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,893,920 C→T G107S (GGC→AGC)  yieG ← putative membrane/transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,893,9200CT89.1% 1016.1 / 66.5 331G107S (GGC→AGC) yieGputative membrane/transport protein
Reads supporting (aligned to +/- strand):  ref base C (17/19);  new base T (157/137);  total (175/156)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 333 total aligned reads displayed.

TATCATCCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCA  >  NC_000913/3893828‑3894009
                                                                                            |                                                                                          
  tcatcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCa                                                                                          >  1:6253365/1‑93 (MQ=255)
  tcatcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCa                                                                                          <  2:3933156/93‑1 (MQ=255)
   catcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCAt                                                                                         <  1:5532518/93‑1 (MQ=255)
    atcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCTCCCATCCCGACCTGCCACGGCAAGCCCATc                                                                                        >  2:7436493/1‑93 (MQ=255)
      ccAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGc                                                                                      <  1:3668093/93‑1 (MQ=255)
         gTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTg                                                                                   <  2:853448/93‑1 (MQ=255)
          tAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGt                                                                                  >  1:1791032/1‑93 (MQ=255)
           aGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTa                                                                                 >  1:6447362/1‑93 (MQ=255)
             cGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACa                                                                               <  2:4198649/93‑1 (MQ=255)
              gAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACaa                                                                              >  1:2746412/1‑93 (MQ=255)
               aaCGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAAc                                                                             >  2:6365408/1‑93 (MQ=255)
                aCGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACg                                                                            >  1:2527883/1‑93 (MQ=255)
                 cgcgGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGa                                                                           >  1:1181574/1‑93 (MQ=255)
                  gcgGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGAc                                                                          >  2:290293/1‑93 (MQ=255)
                  gcgGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGAc                                                                          <  2:6510187/93‑1 (MQ=255)
                   cgGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACa                                                                         >  2:456624/1‑93 (MQ=255)
                    ggAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACaa                                                                        >  2:3986983/1‑93 (MQ=255)
                     gAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACaaa                                                                       >  2:5769289/1‑93 (MQ=255)
                      aaaaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACaaaa                                                                      <  2:3453095/93‑1 (MQ=255)
                      aaaaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACaaaa                                                                      <  1:5222575/93‑1 (MQ=255)
                      aaaaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACaaaa                                                                      >  1:3752116/1‑93 (MQ=255)
                      aaaaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACaaaa                                                                      <  2:5021399/93‑1 (MQ=255)
                       aaaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAg                                                                     >  1:4865056/1‑93 (MQ=255)
                       aaaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAg                                                                     >  1:7417550/1‑93 (MQ=255)
                        aaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGc                                                                    <  2:1708922/93‑1 (MQ=255)
                        aaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGc                                                                    <  1:5330655/93‑1 (MQ=255)
                         aTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCg                                                                   >  1:4116894/1‑93 (MQ=255)
                         aTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCg                                                                   >  1:4447924/1‑93 (MQ=255)
                          tCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCga                                                                  >  2:4933212/1‑93 (MQ=255)
                          tCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCga                                                                  >  2:5701133/1‑93 (MQ=255)
                            gTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCgaag                                                                >  2:2986009/1‑93 (MQ=255)
                            gTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCgaag                                                                >  1:3058956/1‑93 (MQ=255)
                              cAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCgaagaa                                                              >  1:2947635/1‑93 (MQ=255)
                              cAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCgaagaa                                                              <  2:6314028/93‑1 (MQ=255)
                                gTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACg                                                            <  2:92917/93‑1 (MQ=255)
                                 tAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGc                                                           <  1:6603181/93‑1 (MQ=255)
                                 tAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGc                                                           <  2:2543306/93‑1 (MQ=255)
                                  aaCAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCa                                                          <  2:1965364/93‑1 (MQ=255)
                                   aNAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCan                                                         <  2:6035478/93‑2 (MQ=255)
                                    cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCAtt                                                        <  2:364072/93‑1 (MQ=255)
                                    cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCAtt                                                        <  1:3518382/93‑1 (MQ=255)
                                    cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCAtt                                                        <  2:5514968/93‑1 (MQ=255)
                                     agcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTc                                                       >  2:1085462/1‑93 (MQ=255)
                                      gcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCa                                                      >  1:3967102/1‑93 (MQ=255)
                                      gcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCa                                                      <  2:2591240/93‑1 (MQ=255)
                                      gcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCa                                                      >  2:1832478/1‑93 (MQ=255)
                                        agACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAgg                                                   <  1:2681430/93‑1 (MQ=255)
                                         gACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGc                                                  >  2:4511821/1‑93 (MQ=255)
                                         gACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGc                                                  >  1:7077272/1‑93 (MQ=255)
                                          aCCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGcc                                                 >  1:5697596/1‑93 (MQ=255)
                                           ccTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGccc                                                >  1:3232110/1‑93 (MQ=255)
                                            cTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCa                                               <  2:5232993/93‑1 (MQ=255)
                                            cTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCa                                               >  2:416178/1‑93 (MQ=255)
                                              aTCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATa                                             <  1:1328096/93‑1 (MQ=255)
                                               tCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAg                                            <  2:6294830/93‑1 (MQ=255)
                                                 gcgcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCg                                          <  2:6103244/93‑1 (MQ=255)
                                                  cgcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCgg                                         >  2:5294256/1‑93 (MQ=255)
                                                   gcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCggg                                        <  1:2327702/93‑1 (MQ=255)
                                                   gcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCggg                                        >  1:6387993/1‑93 (MQ=255)
                                                    ccccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGt                                       <  1:5054506/93‑1 (MQ=255)
                                                      ccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGc                                     <  2:2017907/93‑1 (MQ=255)
                                                      ccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGc                                     >  2:807210/1‑93 (MQ=255)
                                                        agaagaTTGCGCCCATCCCGACCTGCCCCGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCa                                   >  1:4930138/1‑93 (MQ=255)
                                                        agaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCa                                   >  1:2150320/1‑93 (MQ=255)
                                                         gaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGGCCCATAGCGGGTGCCa                                   >  1:6155073/1‑93 (MQ=255)
                                                         gaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAg                                  >  1:5770933/1‑93 (MQ=255)
                                                         gaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAg                                  >  1:6332393/1‑93 (MQ=255)
                                                         gaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAg                                  >  1:2415806/1‑93 (MQ=255)
                                                         gaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAg                                  <  1:7238912/93‑1 (MQ=255)
                                                          aagaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGt                                 <  1:2618276/93‑1 (MQ=255)
                                                           agaTTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTg                                >  1:2179519/1‑93 (MQ=255)
                                                            gaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGc                               >  1:803962/1‑93 (MQ=255)
                                                               tGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAAc                            <  2:6875080/93‑1 (MQ=255)
                                                                 cgcCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTg                          <  2:7338160/93‑1 (MQ=255)
                                                                  gcCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTgg                         <  1:1127451/93‑1 (MQ=255)
                                                                  gcCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTgg                         <  1:2644509/93‑1 (MQ=255)
                                                                  gcCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTgg                         <  2:4847411/93‑1 (MQ=255)
                                                                  gcCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTgg                         <  2:4675250/93‑1 (MQ=255)
                                                                    ccATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCa                       <  1:1542190/93‑1 (MQ=255)
                                                                    ccATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCa                       <  2:5400241/93‑1 (MQ=255)
                                                                     cATGCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAg                      >  1:6976913/1‑93 (MQ=255)
                                                                     cATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAg                      >  2:3259832/1‑93 (MQ=255)
                                                                      aTCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAgg                     <  2:5877140/93‑1 (MQ=255)
                                                                      aTCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAgg                     <  2:2397060/93‑1 (MQ=255)
                                                                       tCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGt                    >  2:2745282/1‑93 (MQ=255)
                                                                        cccGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGtt                   <  2:4389892/93‑1 (MQ=255)
                                                                         ccGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTa                  <  2:1408073/93‑1 (MQ=255)
                                                                          cGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAg                 <  2:1551506/93‑1 (MQ=255)
                                                                          cGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAg                 >  1:4647599/1‑93 (MQ=255)
                                                                            aCCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCa               <  2:3039977/93‑1 (MQ=255)
                                                                            aCCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCa               >  1:2826257/1‑93 (MQ=255)
                                                                              cTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCaaa             <  1:7146325/93‑1 (MQ=255)
                                                                              cTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCaaa             <  2:1243469/93‑1 (MQ=255)
                                                                               tGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAAc            <  1:6045582/93‑1 (MQ=255)
                                                                                  cACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGt         <  1:3610880/93‑1 (MQ=255)
                                                                                  cACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGt         <  2:6693539/93‑1 (MQ=255)
                                                                                      gCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCca     >  2:7116893/1‑93 (MQ=255)
                                                                                      gCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGGGCCAGTGCAACTGGCAGGGTAGCAAACAGTCCca     >  1:6923136/1‑93 (MQ=255)
                                                                                        aaGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcatc   >  1:4316769/1‑93 (MQ=255)
                                                                                         aGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGGTAGCAAACAAACCcctca  >  1:6662888/1‑93 (MQ=255)
                                                                                            |                                                                                          
TATCATCCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCA  >  NC_000913/3893828‑3894009

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