Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,139,919 T→G 5.2% I16L (ATT→CTT)  hybG ← hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,139,9190TG 5.2% 495.6 / 16.8 135I16L (ATT→CTT) hybGhydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base T (64/64);  new base G (4/3);  total (68/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-02

Only 100 of 142 total aligned reads displayed.

CTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGACACTG  >  NC_000913/3139772‑3140057
                                                                                                                                                   |                                                                                                                                          
cttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCGCTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAat                                                                                                                                            >  4:967380/1‑148 (MQ=255)
cttCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAat                                                                                                                                            <  1:434299/148‑1 (MQ=255)
  tCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAatat                                                                                                                                          >  2:176572/1‑148 (MQ=255)
      cTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTc                                                                                                                                      >  2:49887/1‑148 (MQ=255)
      cTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTc                                                                                                                                      <  1:176572/148‑1 (MQ=255)
            cgATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGAc                                                                                                                                >  2:99252/1‑148 (MQ=255)
               tgatgCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGAcagc                                                                                                                             <  1:49887/148‑1 (MQ=255)
                      cATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAat                                                                                                                                            <  1:579635/126‑1 (MQ=255)
                              aTCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcctggcctgg                                                                                                              >  4:1066983/1‑148 (MQ=255)
                              aTCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAAGATATTCAc                                                                                                                                    >  2:177433/1‑126 (MQ=255)
                                  gACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACg                                                                                                          <  2:107807/148‑1 (MQ=255)
                                              cACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACatt                                                                                              <  3:880696/148‑1 (MQ=255)
                                                cccACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACattat                                                                                            <  2:865867/148‑1 (MQ=255)
                                                  cACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcctggcc                                                                                                                 <  1:1015415/125‑1 (MQ=255)
                                                  cACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcctggcc                                                                                                                 >  2:1015415/1‑125 (MQ=255)
                                                  cACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcc                                                                                                                      <  2:1126036/120‑1 (MQ=255)
                                                   aCTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTaa                                                                                         >  2:38848/1‑148 (MQ=255)
                                                      ggCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTc                                                                                      <  1:259030/148‑1 (MQ=255)
                                                        cccAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCa                                                                                                       >  2:410825/1‑129 (MQ=255)
                                                          cAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGAcagccagc                                                                                                                         >  2:133162/1‑109 (MQ=255)
                                                          cAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGAcagccagc                                                                                                                         <  1:133162/109‑1 (MQ=255)
                                                          cAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGa                                                                                                             >  4:525827/1‑121 (MQ=255)
                                                          cAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGa                                                                                                             <  3:525827/121‑1 (MQ=255)
                                                           aGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGtt                                                                                 <  3:172793/148‑1 (MQ=255)
                                                             tAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCATCTCCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTAt                                                                               <  2:587685/148‑1 (MQ=255)
                                                               gATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcctggc                                                                                                                  >  4:885912/1‑111 (MQ=255)
                                                                    gCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCa                                                                            <  3:1018772/144‑1 (MQ=255)
                                                                     cAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAAGATATTCAcc                                                                                                                                   <  1:61533/88‑1 (MQ=255)
                                                                     cAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAAGATATTCAcc                                                                                                                                   >  2:61533/1‑88 (MQ=255)
                                                                      aGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCa                                                                      <  3:336956/148‑1 (MQ=255)
                                                                           tACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcctggcc                                                                                                                 >  4:15083/1‑100 (MQ=255)
                                                                            aCCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACtt                                                                <  1:1032650/148‑1 (MQ=255)
                                                                              cTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATaca                                                                                                  >  3:139381/1‑112 (MQ=255)
                                                                              cTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTAtt                                                                              >  1:270752/1‑132 (MQ=255)
                                                                              cTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTAtt                                                                              <  2:270752/132‑1 (MQ=255)
                                                                               ttCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACg                                                                                                          <  1:366213/103‑1 (MQ=255)
                                                                                tCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATaca                                                                                                  >  3:767836/1‑110 (MQ=255)
                                                                                tCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATa                                                                                                    >  4:354197/1‑108 (MQ=255)
                                                                                tCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGa                                                            <  1:455871/148‑1 (MQ=255)
                                                                                 cacaAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcctg                                                                                                                    >  1:926226/1‑91 (MQ=255)
                                                                                 cacaAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCAcctg                                                                                                                    <  2:926226/91‑1 (MQ=255)
                                                                                   caAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGt                                                                                  >  2:1118369/1‑123 (MQ=255)
                                                                                   caAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGt                                                                                  <  1:1118369/123‑1 (MQ=255)
                                                                                   caAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCa                                                                            >  1:119727/1‑129 (MQ=255)
                                                                                   caAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCa                                                                            <  4:723155/129‑1 (MQ=255)
                                                                                   caAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAAc                                                                          <  4:1115055/131‑1 (MQ=255)
                                                                                   caAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAAc                                                                          >  3:1115055/1‑131 (MQ=255)
                                                                                   caAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGa                                                         >  2:945217/1‑148 (MQ=255)
                                                                                   caAATCAGAGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAAGATATTCAc                                                                                                                                    >  1:965964/1‑73 (MQ=255)
                                                                                    aaaTCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGAt                                                        >  3:801350/1‑148 (MQ=255)
                                                                                     aaTCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCg                                                                                    >  4:323117/1‑119 (MQ=255)
                                                                                     aaTCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCg                                                                                    <  3:323117/119‑1 (MQ=255)
                                                                                          gggCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGAt                                                        >  2:536045/1‑142 (MQ=255)
                                                                                          gggCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGAt                                                        <  1:536045/142‑1 (MQ=255)
                                                                                           ggCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGaa                                                 <  1:38848/148‑1 (MQ=255)
                                                                                                      caTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATaca                                                                                                  >  3:327146/1‑88 (MQ=255)
                                                                                                      caTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATaca                                                                                                  <  4:327146/88‑1 (MQ=255)
                                                                                                      caTCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTAt                                                                               <  4:211711/107‑1 (MQ=255)
                                                                                                          gcgcTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGaca                                  <  4:801350/148‑1 (MQ=255)
                                                                                                           cgcTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCaa                                                                           <  4:744926/106‑1 (MQ=255)
                                                                                                           cgcTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCaa                                                                           >  3:744926/1‑106 (MQ=255)
                                                                                                               tGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACttt                                                               <  2:989146/114‑1 (MQ=255)
                                                                                                               tGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACttt                                                               >  1:989146/1‑114 (MQ=255)
                                                                                                                     cacaTACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACt                                                                 <  4:888266/106‑1 (MQ=255)
                                                                                                                     cacaTACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACt                                                                 >  3:888266/1‑106 (MQ=255)
                                                                                                                     cacaTACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGc                                             <  1:1049914/126‑1 (MQ=255)
                                                                                                                     cacaTACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGc                                             >  2:1049914/1‑126 (MQ=255)
                                                                                                                     cacaTACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTc                       >  2:749523/1‑148 (MQ=255)
                                                                                                                         tACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCTGTTATTCa                                                                            <  4:982189/91‑1 (MQ=255)
                                                                                                                         tACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCa                                                                            >  3:982189/1‑91 (MQ=255)
                                                                                                                         tACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCa                                                                      >  3:659684/1‑97 (MQ=255)
                                                                                                                         tACTTCAACCAGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATc                                                       >  3:917276/1‑112 (MQ=255)
                                                                                                                            ttCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTAtt                                                                              <  1:927532/86‑1 (MQ=255)
                                                                                                                            ttCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTAtt                                                                              >  2:927532/1‑86 (MQ=255)
                                                                                                                             tCAACCTGCGCAAGCTGGTGAATATCTTCACCGAcag                                                                                                                              <  2:1057553/37‑1 (MQ=255)
                                                                                                                             tCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACg                                                                                                          >  2:345564/1‑57 (MQ=255)
                                                                                                                             tCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCg                                                                                    >  1:992532/1‑79 (MQ=255)
                                                                                                                             tCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCa                                                                            <  1:377683/87‑1 (MQ=255)
                                                                                                                             tCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCa                                                                            >  2:377683/1‑87 (MQ=255)
                                                                                                                             tCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACtt                                                                <  4:774218/99‑1 (MQ=255)
                                                                                                                             tCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTg                                                             >  1:299678/1‑102 (MQ=255)
                                                                                                                             tCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTg                                                             >  3:1048113/1‑102 (MQ=255)
                                                                                                                             tCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTg                                                             <  4:1048113/102‑1 (MQ=255)
                                                                                                                             tCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGa                                                            <  1:755734/103‑1 (MQ=255)
                                                                                                                                 ccTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTg                                                             <  1:907941/98‑1 (MQ=255)
                                                                                                                                  cTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGc                                             <  2:317449/113‑1 (MQ=255)
                                                                                                                                  cTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGc                                             >  1:317449/1‑113 (MQ=255)
                                                                                                                                   tGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTg                                                             <  3:152800/96‑1 (MQ=255)
                                                                                                                                   tGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTg                                                             >  4:152800/1‑96 (MQ=255)
                                                                                                                                    gcgcAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCgg        <  3:991930/148‑1 (MQ=255)
                                                                                                                                     cgcAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGa                                                            <  2:711527/95‑1 (MQ=255)
                                                                                                                                     cgcAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGa                                                            >  1:711527/1‑95 (MQ=255)
                                                                                                                                       cAAGCTGGTGACTATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATAcaca                                                                                                >  3:80454/1‑57 (MQ=255)
                                                                                                                                       cAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATac                                                                                                   <  2:749763/54‑1 (MQ=255)
                                                                                                                                       cAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATAcaca                                                                                                <  4:80454/57‑1 (MQ=255)
                                                                                                                                          gtcggTGAAGATATTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGACACTg  >  3:1113934/4‑148 (MQ=255)
                                                                                                                                             ggTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGacac                                 <  2:296887/114‑1 (MQ=255)
                                                                                                                                               tGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGaca                                  <  1:859037/111‑1 (MQ=255)
                                                                                                                                               tGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGaca                                  >  2:859037/1‑111 (MQ=255)
                                                                                                                                                  atatCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGa                                                            >  1:108106/1‑82 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CTTCATCTTCGTCGATGATGCTCATGGCAAATCCGACGTGTACCAGCACCCACTGGCCCAGTAGATCGGCAGGGTTACCTTCACAAATCAGGGCGATATTCACATCGCGCTTGATACCACATACTTCAACCTGCGCAAGCTGGTGAATATCTTCACCGACAGCCAGCACCTGGCCTGGAACGCCAATACACATTATTAACTCCGGTTATTCAACTTCAATACTTTTGACGATCAGCGAATCGCCGGTATCGACACGCAACCGCTCGCCGTGACAGAGCGGACACTG  >  NC_000913/3139772‑3140057

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