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,9200CT99.7% 1280.0 / ‑6.3 360G107S (GGC→AGC) yieGputative membrane/transport protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (192/166);  total (193/167)
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.99e-01

Only 100 of 370 total aligned reads displayed.

TCATCCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCA  >  NC_000913/3893830‑3894009
                                                                                          |                                                                                         
tcatcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCa                                                                                         >  2:4532942/1‑93 (MQ=255)
tcatcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCa                                                                                         >  2:686017/1‑93 (MQ=255)
 catcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCAt                                                                                        >  1:4474885/1‑93 (MQ=255)
 catcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCAt                                                                                        <  1:3265131/93‑1 (MQ=255)
  atcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATc                                                                                       >  1:488224/1‑93 (MQ=255)
  atcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATc                                                                                       <  2:4844272/93‑1 (MQ=255)
  atcCAGTAGCGAACGAGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATc                                                                                       <  1:5445785/93‑1 (MQ=255)
   tcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCg                                                                                      >  1:8264438/1‑93 (MQ=255)
   tcCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCg                                                                                      >  2:7425760/1‑93 (MQ=255)
      aGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCt                                                                                   <  2:2309538/93‑1 (MQ=255)
      aGTAGCGAACGCGAAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCt                                                                                   <  2:7714054/93‑1 (MQ=255)
        tAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGt                                                                                 <  1:5577323/93‑1 (MQ=255)
          gCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTAc                                                                               <  1:6083068/93‑1 (MQ=255)
          gCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTAc                                                                               <  1:2798643/93‑1 (MQ=255)
          gCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTAc                                                                               >  2:6170205/1‑93 (MQ=255)
           cGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACa                                                                              >  2:7248903/1‑93 (MQ=255)
           cGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACa                                                                              >  1:6740013/1‑93 (MQ=255)
              aCGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAAca                                                                           >  2:7263200/1‑92 (MQ=255)
                gcgGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACaacaac                                                                         >  1:3107125/1‑93 (MQ=255)
                gcgGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGAc                                                                         >  1:6679357/1‑93 (MQ=255)
                     aaaTCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAg                                                                    >  2:7944615/1‑93 (MQ=255)
                     aaaTCGTCAGTAACAGCAGACCTATCACGCCACATAAGATTGAGCCCATCCCGACCTGCCAAGGCAAGCTCATAGCCTGTACAACGACAAAAg                                                                    >  2:3172832/1‑93 (MQ=255)
                         cGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCgaa                                                                >  2:4079057/1‑93 (MQ=255)
                         cGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCgaa                                                                <  1:8455449/93‑1 (MQ=255)
                            cAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCgaagaa                                                             >  2:8325064/1‑93 (MQ=255)
                            cAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCgaagaa                                                             <  1:7935273/93‑1 (MQ=255)
                             aGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAAc                                                            <  1:3640758/93‑1 (MQ=255)
                             aGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAAc                                                            <  1:1246055/93‑1 (MQ=255)
                               tAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGc                                                          <  2:5532526/93‑1 (MQ=255)
                               tAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGc                                                          >  2:6227848/1‑93 (MQ=255)
                                aaCAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCa                                                         <  2:1538087/93‑1 (MQ=255)
                                aaCAGCAGACCTATCCCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGCAGAACGCa                                                         >  2:8709537/1‑93 (MQ=255)
                                  cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCAtt                                                       >  2:2852803/1‑93 (MQ=255)
                                  cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCAtt                                                       >  1:2460680/1‑93 (MQ=255)
                                  cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCAtt                                                       >  2:7489845/1‑93 (MQ=255)
                                  cagcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCAtt                                                       <  1:1545896/93‑1 (MQ=255)
                                   agcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGCACAACGACAAAAGCGAAGAACGCATTc                                                      <  1:1097350/93‑1 (MQ=255)
                                    gcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCa                                                     <  2:565933/93‑1 (MQ=255)
                                    gcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCa                                                     <  2:7649376/93‑1 (MQ=255)
                                    gcagACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCa                                                     <  2:3402143/93‑1 (MQ=255)
                                     cagACCTAACGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAg                                                    <  2:8530038/93‑1 (MQ=255)
                                        aCCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGcc                                                 >  2:198388/1‑93 (MQ=255)
                                         ccTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGccc                                                <  2:7196386/93‑1 (MQ=255)
                                          cTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCa                                               <  2:3010087/93‑1 (MQ=255)
                                          cTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCa                                               <  2:1862896/93‑1 (MQ=255)
                                          cTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCa                                               >  2:8449142/1‑93 (MQ=255)
                                             tCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAg                                            <  1:6800703/93‑1 (MQ=255)
                                              cgcgCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCACATAGc                                           <  2:7795408/93‑1 (MQ=255)
                                              cgcgCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGc                                           <  2:6340978/93‑1 (MQ=255)
                                              cgcgCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGc                                           <  1:8088756/93‑1 (MQ=255)
                                               gcgcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGAATTAAGGCCCATACCg                                          >  2:3283900/1‑93 (MQ=38)
                                               ccgcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCg                                          <  1:2342462/92‑1 (MQ=255)
                                                cgcCCCAGAAGATTGGGCCCATGCCGGGCTGCCACAGCAAGCTCATCGCCTGTACAACAACAGACGCCAAGAACGCATTCACGCCCCTAgtcg                                         >  2:7373877/1‑90 (MQ=255)
                                                cgcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCgg                                         >  2:1378945/1‑93 (MQ=255)
                                                cgcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCgg                                         <  1:8184093/93‑1 (MQ=255)
                                                 gcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCggg                                        >  1:4326595/1‑93 (MQ=255)
                                                 gcCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCggg                                        >  1:7324050/1‑93 (MQ=255)
                                                   cccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACACCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTg                                      >  1:3734653/1‑93 (MQ=255)
                                                   cccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTg                                      >  2:5972904/1‑93 (MQ=255)
                                                   cccAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTg                                      >  1:7372868/1‑93 (MQ=255)
                                                     cAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGcc                                    >  1:5983234/1‑93 (MQ=255)
                                                      agaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCa                                   >  2:3610582/1‑93 (MQ=255)
                                                      agaagaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCa                                   >  1:7881115/1‑93 (MQ=255)
                                                        aagaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGt                                 >  2:2543633/1‑93 (MQ=255)
                                                         agaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTg                                >  2:504708/1‑93 (MQ=255)
                                                          gaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGc                               <  1:651907/93‑1 (MQ=255)
                                                          gaTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGc                               <  2:1027278/93‑1 (MQ=255)
                                                           aTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCa                              >  1:4044625/1‑93 (MQ=255)
                                                           aTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCa                              >  2:5361398/1‑93 (MQ=255)
                                                           aTTGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCa                              <  2:5226941/93‑1 (MQ=255)
                                                             tGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAAc                            <  1:1000408/93‑1 (MQ=255)
                                                             tGCGCCCATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAAc                            <  2:1293916/93‑1 (MQ=255)
                                                                  ccATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCa                       >  2:366097/1‑93 (MQ=255)
                                                                   cATCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAg                      <  1:269295/93‑1 (MQ=255)
                                                                    aTCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAgg                     >  1:1724580/1‑93 (MQ=255)
                                                                     tCCCGACCTGCCACGGCAAGCTCATCGCCTGTATAACGACAAAAGCGAAGAACGCATTTAGGCCAATAGCGGGAGCCAGAGCAACTGGCAGGt                    >  2:4620971/1‑93 (MQ=255)
                                                                     tCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGt                    >  1:8729868/1‑93 (MQ=255)
                                                                     tCCCGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGt                    <  1:3778715/93‑1 (MQ=255)
                                                                     acccGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGt                    <  1:2072159/92‑1 (MQ=255)
                                                                       ccGACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTa                  <  1:4638208/93‑1 (MQ=255)
                                                                        cGACCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAg                 <  1:2001718/93‑1 (MQ=255)
                                                                         gACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGc                >  2:4311821/1‑93 (MQ=255)
                                                                         gACCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGc                >  1:8561460/1‑93 (MQ=255)
                                                                         caCCTGCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGc                <  1:7301709/92‑1 (MQ=255)
                                                                          cccTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCa               <  1:7234699/92‑1 (MQ=255)
                                                                          aCCTGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCa               <  2:4428081/93‑1 (MQ=255)
                                                                            cggCCACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCaaa             >  2:6100575/3‑93 (MQ=255)
                                                                             tGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAAc            >  1:1375573/1‑93 (MQ=255)
                                                                             tGCCACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAAc            >  1:5329760/1‑93 (MQ=255)
                                                                               ccACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAg          >  1:6287696/1‑93 (MQ=255)
                                                                               ccACGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAg          <  2:5001707/93‑1 (MQ=255)
                                                                               ccACGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAg          <  1:7747404/93‑1 (MQ=255)
                                                                                 aCGGCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTc        <  2:6485240/93‑1 (MQ=255)
                                                                                  cGGCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTcc       <  1:3451734/93‑1 (MQ=255)
                                                                                   ggCAAGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTccc      <  1:7024508/93‑1 (MQ=255)
                                                                                   ggCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTccc      >  1:8812897/1‑93 (MQ=255)
                                                                                   ggCAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTccc      <  2:2668392/93‑1 (MQ=255)
                                                                                     cAAGCTCATCGCCTGTACAACAACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcat    <  1:5710488/93‑1 (MQ=255)
                                                                                       aGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcatca  >  2:3768415/1‑93 (MQ=255)
                                                                                       aGCTCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCcatca  >  2:423059/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCATCCAGTAGCGAACGCGGAAAATCGTCAGTAACAGCAGACCTATCGCGCCCCAGAAGATTGCGCCCATCCCGACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCA  >  NC_000913/3893830‑3894009

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