Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,273,229 A→T 15.3% intergenic (+27/+17) chaC → / ← ychN cation transport regulator/putative sulfur relay protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,273,2290AT15.3% 217.0 / 32.9 85intergenic (+27/+17)chaC/ychNcation transport regulator/putative sulfur relay protein
Reads supporting (aligned to +/- strand):  ref base A (39/33);  new base T (3/10);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 6.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

TTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAG  >  NC_000913/1273130‑1273326
                                                                                                   |                                                                                                 
tttGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCa                                                                                                   >  1:149085/1‑100 (MQ=255)
  tGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTCAACTTCCCGCATAGTGGGGCGTCAGa                                                                                                 >  1:56521/1‑100 (MQ=255)
  tGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGa                                                                                                 <  1:1127179/100‑1 (MQ=255)
    gtGGTATCGTTAAGAAAACTGCTGGCGGAGAATTTTGCGGATGGTGTGTGACGGCCGGGATTCGCCTGAGTAAACTTTCGGCATAGTGGGGCTGCAGACg                                                                                               >  2:228786/2‑100 (MQ=255)
        tATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGcccc                                                                                           >  1:1354763/1‑100 (MQ=255)
         aTCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGccccc                                                                                          <  1:44214/100‑1 (MQ=255)
          tCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccac                                                                                         <  1:1813574/100‑3 (MQ=255)
          tCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccac                                                                                         <  1:1154476/100‑3 (MQ=255)
          tCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccac                                                                                         <  2:956974/100‑3 (MQ=255)
           cGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccact                                                                                        <  2:1776823/100‑4 (MQ=255)
           cGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTc                                                                                        >  1:552186/1‑100 (MQ=255)
           cGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTc                                                                                        >  1:1190023/1‑100 (MQ=255)
            ggTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCa                                                                                       >  1:372190/1‑100 (MQ=255)
             gTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccactat                                                                                      <  2:1766371/100‑6 (MQ=255)
              tAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCaaa                                                                                     >  1:1402399/1‑100 (MQ=255)
              tAAAAAAACTGCTGGCGGAGAAGTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccactatg                                                                                     <  1:2100328/100‑7 (MQ=255)
               aaaaaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccactatgc                                                                                    <  1:980601/100‑8 (MQ=255)
                 aaaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccactatgcgg                                                                                  <  1:1377454/100‑10 (MQ=39)
                 aaaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccactatgcgg                                                                                  <  1:1027936/100‑10 (MQ=39)
                 aaaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACAt                                                                                  <  2:2137542/100‑1 (MQ=255)
                  aaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccactatgcggg                                                                                 <  2:539820/100‑11 (MQ=38)
                     aCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTaaa                                                                              <  2:898983/100‑1 (MQ=255)
                      ctgctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTAGGGCGTCAGACGCCCCCTCAAACATTaaaa                                                                             >  1:1184856/1‑100 (MQ=255)
                           ggCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGa                                                                        <  1:766537/100‑1 (MQ=255)
                            gCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAg                                                                       <  2:342144/100‑1 (MQ=255)
                            gCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAg                                                                       <  2:358625/100‑1 (MQ=255)
                             cGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGc                                                                      >  1:1415177/1‑100 (MQ=255)
                               ccgatctttCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCAc                                                                    <  1:759290/94‑1 (MQ=255)
                                   aTTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTa                                                                >  1:1799281/1‑100 (MQ=255)
                                   aTTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTa                                                                <  1:1391071/100‑1 (MQ=255)
                                   aTTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTa                                                                >  2:246886/1‑100 (MQ=255)
                                   aTTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGAGTCAGAAGCCCCCTCAAACATTAAAATGTGAGCACTTTa                                                                >  2:653481/1‑100 (MQ=255)
                                    ttttCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTAt                                                               <  2:1579437/100‑1 (MQ=255)
                                     tttCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATc                                                              >  2:1374610/1‑100 (MQ=255)
                                      ttCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCAcagatcgg                                                             >  2:759288/1‑93 (MQ=255)
                                      ttCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCg                                                             >  2:895219/1‑100 (MQ=255)
                                        ccGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGc                                                           <  2:552184/100‑1 (MQ=255)
                                        ccGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGc                                                           <  1:1374618/100‑1 (MQ=255)
                                         cGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCt                                                          >  2:931382/1‑100 (MQ=255)
                                         cGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCt                                                          >  1:1596344/1‑100 (MQ=255)
                                          ggATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTg                                                         <  1:1085094/100‑1 (MQ=255)
                                          ggATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTg                                                         <  1:1132409/100‑1 (MQ=255)
                                          ggATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTg                                                         <  2:329068/100‑1 (MQ=255)
                                            aTGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGAc                                                       <  1:1524551/100‑1 (MQ=255)
                                            aTGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGAc                                                       <  2:613078/100‑1 (MQ=255)
                                             tGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACa                                                      <  1:993443/100‑1 (MQ=255)
                                               gtgtgtTACGTACGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGCGAGCACTTTATCGGCTGACAGc                                                    >  2:1830944/1‑100 (MQ=255)
                                                  tgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTc                                                 >  1:963307/1‑100 (MQ=255)
                                                   gtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTcc                                                >  2:842549/1‑100 (MQ=255)
                                                   gtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTcc                                                >  2:1745085/1‑100 (MQ=255)
                                                   gtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTcc                                                >  1:930089/1‑100 (MQ=255)
                                                   gtTACGTCCGGGAGTCGCCTGAGTAAACTGCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTcc                                                <  2:2004321/100‑1 (MQ=255)
                                                    ttACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGAGTCAGAAGCCCCCTCAAACATTAAAATGTGGGCACTTTATCGGCTGACAGCGTCCa                                               >  2:1520155/1‑100 (MQ=255)
                                                         tCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTgcg                                          >  1:2217854/1‑100 (MQ=255)
                                                          ccGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTgcgc                                         <  2:762094/100‑1 (MQ=255)
                                                          ccGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTgcgc                                         <  2:93642/100‑1 (MQ=255)
                                                           cGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTgcg                                          <  2:182956/98‑1 (MQ=255)
                                                            gggATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCa                                       >  2:1329042/1‑100 (MQ=255)
                                                                ttCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTc                                   >  1:2053611/1‑100 (MQ=255)
                                                                 tCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTcc                                  <  1:2263594/100‑1 (MQ=255)
                                                                  cGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCca                                 >  1:1141611/1‑100 (MQ=255)
                                                                   gCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCcac                                <  2:1549178/100‑1 (MQ=255)
                                                                       gAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAaga                            >  1:2273557/1‑100 (MQ=255)
                                                                           aaaCTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCGGTTCCACAAGAGTAc                        >  1:2354905/1‑100 (MQ=255)
                                                                           aaaCTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTAc                        >  1:811172/1‑100 (MQ=255)
                                                                           aaaCTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTAc                        >  2:574746/1‑100 (MQ=255)
                                                                              cTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCCCAAGAGTACCGa                     >  1:1262956/1‑100 (MQ=255)
                                                                              cTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGa                     >  2:207222/1‑100 (MQ=255)
                                                                              cTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGa                     >  2:323183/1‑100 (MQ=255)
                                                                                tCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGAtt                   <  1:2082123/100‑1 (MQ=255)
                                                                                 cccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACGAGAGTACCGAttt                  >  2:2104455/1‑100 (MQ=255)
                                                                                  cctcagaGTGGCGCGTCAGTAGTCCCCTCAATCACTAAACTCTGAACACTTTATCGGCTGTCAGCGTCCTTTTCGCCAGTTCCACAAGAGTACCGAttt                  <  2:1732470/93‑1 (MQ=255)
                                                                                  ccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTc                 >  1:664485/1‑100 (MQ=255)
                                                                                   cGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTcc                >  1:2050541/1‑100 (MQ=255)
                                                                                      ttgagggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcc             >  1:1050727/6‑100 (MQ=255)
                                                                                      aTAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcc             <  1:161002/100‑1 (MQ=255)
                                                                                       tgagggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAccc            >  1:2332541/5‑100 (MQ=255)
                                                                                        aGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcccc           <  1:937755/100‑1 (MQ=255)
                                                                                        aGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcccc           <  1:940515/100‑1 (MQ=255)
                                                                                         agggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCa          >  2:1011234/3‑100 (MQ=255)
                                                                                           ggggCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatc        <  2:253949/100‑1 (MQ=255)
                                                                                           ggggCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatc        <  2:1027559/100‑1 (MQ=255)
                                                                                           ggggCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatc        <  2:963304/100‑1 (MQ=255)
                                                                                             ggCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTGCCACAAGAGTACCGATTTCCACCCCatcaa      >  1:704271/1‑100 (MQ=255)
                                                                                              gCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatcaat     <  1:2179295/100‑1 (MQ=255)
                                                                                                 tCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCag  >  1:52754/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
TTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAG  >  NC_000913/1273130‑1273326

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