Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,143,921 A→C 5.8% intergenic (+160/+52) 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,9210AC 5.8% 222.0 / 4.0 69intergenic (+160/+52)cusF/PP_5734copper efflux protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/33);  new base C (2/2);  total (34/35)
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.72e-01

TCCTGTCTGCCAACGCTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGACCATGCCTGATCC  >  NC_002947/6143774‑6144069
                                                                                                                                                   |                                                                                                                                                    
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)
                 tGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCg                                                                                                                                    >  7:94029/1‑149 (MQ=255)
                    ttCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGCCACATCGCTTCAGCgccgtgt                                                                                                                                            >  3:96774/1‑132 (MQ=255)
                    ttCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGCCACATCGCTTCAGCgccgtgt                                                                                                                                            <  4:96774/138‑7 (MQ=255)
                       gCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTaa                                                                                                                               <  1:313314/148‑1 (MQ=255)
                         ccATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAggc                                                                                                                            >  7:262628/1‑149 (MQ=255)
                            tGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAggcggc                                                                                                                         <  7:183933/149‑1 (MQ=255)
                            tGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAggcggc                                                                                                                         <  8:94029/149‑1 (MQ=255)
                               cctcctGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGCCACATCGCTTCAGCGCCGTGTCAg                                                                                                                                         >  5:226441/1‑130 (MQ=255)
                               cctcctGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGCCACATCGCTTCAGCGCCGTGTCAg                                                                                                                                         <  6:226441/130‑1 (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                                                                                                                   <  6:126440/135‑1 (MQ=255)
                                                gCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGagcag                                                                                                                   >  5:126440/1‑135 (MQ=255)
                                                       gAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGAtgtg                                                                                                                                             <  8:35316/102‑1 (MQ=255)
                                                       gAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGAtgtg                                                                                                                                             >  7:35316/1‑102 (MQ=255)
                                                       gAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCa                                                                                                                                          <  8:78060/105‑1 (MQ=255)
                                                       gAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCa                                                                                                                                          >  7:78060/1‑105 (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                                                                       <  7:269307/119‑1 (MQ=255)
                                                                                                            ggCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGa                                                                       >  8:269307/1‑119 (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)
                                                                                                                                                  gAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGACCATGCCTGATc   >  3:318273/1‑149 (MQ=255)
                                                                                                                                                   aaGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGACCATGCCTGATcc  >  6:17877/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCTGTCTGCCAACGCTTGATTCGCCCATGACCTCCTGACGCCTCCTAGCAACTGGAAAAGCTTGGAAGCATCTAGGAAGTCAGGGCGATTCAGATCATGGTAACCCCGGCCGCTGGAGCCATTTCTGAAGTACTGACACGGCGCTGAAGCGATGTGGCAGTTCCGCCTAAGGCGGCGAGCAGCAGCGAGAAAAGATCGTCACCCCAAGTTTTCTTCTTGCAGGCGACTGAATAGAGGAGAAATCGAGATAACAGCTCGCAAGGCATCCTGATTGGCGCCCGACCATGCCTGATCC  >  NC_002947/6143774‑6144069

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: