Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,143,915 G→T 5.8% intergenic (+154/+58) cusF → / ← PP_5734 copper efflux protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,143,9150GT 5.8% 219.5 / 3.7 69intergenic (+154/+58)cusF/PP_5734copper efflux protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/35);  new base T (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TTATCCTGTCTGCCAACGCTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGACCATGCC  >  NC_002947/6143771‑6144063
                                                                                                                                                |                                                                                                                                                    
ttATCCTGTCTGCCAACGCTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCt                                                                                                                                                  <  8:188937/149‑1 (MQ=255)
   tCCTGTCTGCCAACGCTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGaa                                                                                                                                               <  8:233060/149‑1 (MQ=255)
      tGTCTGCCAACGCTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCg                                                                                                                                            >  3:253317/1‑149 (MQ=255)
       gTCTGCCAACGCTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAACGCTTGGAAGCATCTAGGAAGTCAGTGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTAGTGACACGGCCGGGAAGCGa                                                                                                                                           >  7:237885/1‑149 (MQ=255)
        tCTGCCAACTGTTGATTCGCCCATGACCTCCTGAGGCCTCCTATCAACTGGAAAAGCTGGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAATCGAt                                                                                                                                          <  7:233768/149‑1 (MQ=255)
        tCTGCCAACGCTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGAt                                                                                                                                          <  4:253317/149‑1 (MQ=255)
          cgCCAACGCTTGATTCGCCCATCACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGAtgt                                                                                                                                        <  2:173436/148‑1 (MQ=255)
            ccAACGGTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGAtgt                                                                                                                                        <  2:171635/147‑1 (MQ=255)
            ccAACGCTTGATTGGCCCAGGTCCTCCTGTAGCCTCCTATCTACTGGAAAAGCTTGGAAGGATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTgg                                                                                                                                      <  8:237885/149‑1 (MQ=255)
             cAACGCTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTgg                                                                                                                                      >  4:117840/1‑148 (MQ=255)
                  cTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGgcg                                                                                                                                                    <  6:125795/129‑1 (MQ=255)
                  cTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGgcg                                                                                                                                                    >  5:125795/1‑129 (MQ=255)
                    tGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCg                                                                                                                              >  7:94029/1‑149 (MQ=255)
                       ttCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGCCACATCGCTTCAGCgccgtgt                                                                                                                                      <  4:96774/138‑7 (MQ=255)
                       ttCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGCCACATCGCTTCAGCgccgtgt                                                                                                                                      >  3:96774/1‑132 (MQ=255)
                          gCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTaa                                                                                                                         <  1:313314/148‑1 (MQ=255)
                            ccATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAggc                                                                                                                      >  7:262628/1‑149 (MQ=255)
                               tGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAggcggc                                                                                                                   <  8:94029/149‑1 (MQ=255)
                               tGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAggcggc                                                                                                                   <  7:183933/149‑1 (MQ=255)
                                  cctcctGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGCCACATCGCTTCAGCGCCGTGTCAg                                                                                                                                   <  6:226441/130‑1 (MQ=255)
                                  cctcctGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGCCACATCGCTTCAGCGCCGTGTCAg                                                                                                                                   >  5:226441/1‑130 (MQ=255)
                                      ctGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAACTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTTGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGagcagc                                                                                                            <  8:325880/149‑1 (MQ=255)
                                           gCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGagcagcagc                                                                                                         >  1:153203/1‑147 (MQ=255)
                                           gCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCga                                                                                                       <  3:255270/149‑1 (MQ=255)
                                                ctAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGaa                                                                                                    <  3:117840/147‑1 (MQ=255)
                                                  aGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGAATCAGATCATGGTAAACCCGGCCGCTGGAGCCATTTCTGAAGTACTCACACGCCGCTGAAGAGATTTGGCAGTTCCTCCTAAGGCGGCGAGAAGCAGCGAAAAAAg                                                                                                 >  4:354196/1‑148 (MQ=255)
                                                   gCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGagcag                                                                                                             >  5:126440/1‑135 (MQ=255)
                                                   gCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGagcag                                                                                                             <  6:126440/135‑1 (MQ=255)
                                                          gAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGAtgtg                                                                                                                                       >  7:35316/1‑102 (MQ=255)
                                                          gAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGAtgtg                                                                                                                                       <  8:35316/102‑1 (MQ=255)
                                                          gAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCa                                                                                                                                    >  7:78060/1‑105 (MQ=255)
                                                          gAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCa                                                                                                                                    <  8:78060/105‑1 (MQ=255)
                                                          gAAAAGCTTCGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCa                                                                                          <  8:218856/147‑1 (MQ=255)
                                                             aaGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGag                                                                                                                <  5:265987/122‑1 (MQ=255)
                                                             aaGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGag                                                                                                                >  6:265987/1‑122 (MQ=255)
                                                             aaGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCAcccc                                                                                      <  5:233043/148‑1 (MQ=255)
                                                             aaGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCAcccc                                                                                      >  6:233043/1‑148 (MQ=255)
                                                                       gCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTttcttc                                                                           <  8:262628/149‑1 (MQ=255)
                                                                         aTCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTcttctt                                                                         <  5:252402/149‑1 (MQ=255)
                                                                           cTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGt                                                                                            >  6:173406/1‑128 (MQ=255)
                                                                             aGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAg                                                                     >  5:156760/1‑149 (MQ=255)
                                                                                     aGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACt                                                               <  1:22054/147‑1 (MQ=255)
                                                                                     aGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACt                                                               >  2:22054/1‑147 (MQ=255)
                                                                                     aGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGa                                                             <  2:153203/149‑1 (MQ=255)
                                                                                     aGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCAATTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTg                                                              <  5:315606/148‑1 (MQ=255)
                                                                                           aTTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCAc                                                                                         >  1:238459/1‑115 (MQ=255)
                                                                                           aTTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCAc                                                                                         <  2:238459/115‑1 (MQ=255)
                                                                                           aTTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATagag                                                       <  4:61476/149‑1 (MQ=255)
                                                                                            ttCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATagag                                                       >  4:357639/1‑148 (MQ=255)
                                                                                              cAGATCATGGTAACCCCGGCCTCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGTCGAGCATCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTTCAGGCGACTGAATAgagga                                                     <  3:354196/148‑1 (MQ=255)
                                                                                               aGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAgaggag                                                    >  1:98014/1‑148 (MQ=255)
                                                                                               aGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAgaggag                                                    >  2:4960/1‑148 (MQ=255)
                                                                                                       gTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCgaga                                           <  6:156760/149‑1 (MQ=255)
                                                                                                           ccccGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATaa                                        >  8:218989/1‑148 (MQ=255)
                                                                                                            cccGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATaa                                        >  6:130352/1‑147 (MQ=255)
                                                                                                               ggCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGa                                                                 >  8:269307/1‑119 (MQ=255)
                                                                                                               ggCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGa                                                                 <  7:269307/119‑1 (MQ=255)
                                                                                                                          ccATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATc                        <  3:357639/149‑1 (MQ=255)
                                                                                                                           cATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGAt                                          <  6:58711/130‑1 (MQ=255)
                                                                                                                           cATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGAt                                          >  5:58711/1‑130 (MQ=255)
                                                                                                                            aTTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCt                      >  8:224941/1‑149 (MQ=255)
                                                                                                                                 tGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCAACTGAATAGAGGAGAAATCGAGAt                                          >  5:191517/1‑124 (MQ=255)
                                                                                                                                 tGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCAACTGAATAGAGGAGAAATCGAGAt                                          <  6:191517/124‑1 (MQ=255)
                                                                                                                                   aaGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGAtt                  <  7:218989/146‑1 (MQ=255)
                                                                                                                                        cTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGgcgcac           >  8:232841/1‑146 (MQ=255)
                                                                                                                                         tGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGa         >  1:238950/1‑149 (MQ=255)
                                                                                                                                           acacGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAg                                                                     <  8:149059/87‑1 (MQ=255)
                                                                                                                                           acacGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAg                                                                     >  7:149059/1‑87 (MQ=255)
                                                                                                                                           acacGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGAc        <  2:98014/148‑1 (MQ=255)
                                                                                                                                             acGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGACCAt     >  5:54649/1‑149 (MQ=255)
                                                                                                                                             acGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGACCAt     >  5:60073/1‑149 (MQ=255)
                                                                                                                                                gcgcTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGACCATGcc  >  5:105872/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTATCCTGTCTGCCAACGCTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGACCATGCC  >  NC_002947/6143771‑6144063

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