Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,583,521 T→C 10.9% P737P (CCA→CCG)  parC ← DNA topoisomerase IV subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,583,5210TC10.9% 165.7 / 9.9 64P737P (CCA→CCG) parCDNA topoisomerase IV subunit A
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base C (3/4);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

GGCACCGCAAGGGGCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACGGCTGG  >  NC_002947/5583375‑5583669
                                                                                                                                                  |                                                                                                                                                    
ggCACCGCAAGGGGCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTgg                                                                                                                                                    <  3:202919/149‑1 (MQ=255)
 ggaCCGCAAGGGGCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTgg                                                                                                                                                    >  7:61676/3‑148 (MQ=255)
     cGCAAGGGGCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGAAACCCATCGACGCGCTGGAACCCGCGTGGCCGCt                                                                                                                                               >  3:276242/1‑149 (MQ=255)
            ggCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGcc                                                                                                                                        >  3:123622/1‑149 (MQ=255)
             gCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAACCGGCGTTGAAATTGCCGTGCCAGACATTACAACCGGTTATGTCGGGACTTCCGCCGGCAACCCCTCGACGCGCGGGAACCCGCGTGGCAGCTTGCTGCcg                                                                                                                                       >  6:189129/1‑149 (MQ=255)
                 atcatcCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGcc                                                                                                                                   <  4:18192/149‑1 (MQ=255)
                     tcCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGAcc                                                                                                                               >  5:187759/1‑149 (MQ=255)
                       cgcgAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGTGCGGCAGCAAGCTGCCacgcgggttcc                                                                                                                             <  4:273707/149‑8 (MQ=255)
                       cgcgAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTtgctg                                                                                                                                          >  7:12068/1‑136 (MQ=255)
                       cgcgAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTtgctg                                                                                                                                          <  8:12068/136‑1 (MQ=255)
                        gagAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTTACGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGAAGGAACCCGCGTGGCAGCTTGCTGCCCCGACGCCcgcg                                                                                                                            >  4:57497/3‑149 (MQ=255)
                          gAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCt                                                                                                                          >  5:323671/1‑149 (MQ=255)
                                    tcctaatCCAGGATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTa                                                                                                                <  3:57497/143‑1 (MQ=255)
                                    tCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTg                                                                                                                  >  8:105246/1‑147 (MQ=255)
                                    tCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTg                                                                                                                  >  5:99336/1‑147 (MQ=255)
                                        aCTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGc                                                                                                            <  1:130931/149‑1 (MQ=255)
                                        aCTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGc                                                                                                            <  4:123622/149‑1 (MQ=255)
                                           ccAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTc                                                                                                          >  5:164740/1‑148 (MQ=255)
                                            cAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTc                                                                                                          <  8:61676/147‑1 (MQ=255)
                                             aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   <  5:188557/121‑2 (MQ=255)
                                             aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   >  6:188557/1‑120 (MQ=255)
                                             aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   >  7:137772/1‑120 (MQ=255)
                                             aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   <  8:137772/121‑2 (MQ=255)
                                             aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   <  1:177682/121‑2 (MQ=255)
                                             aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGCGGCAGCAAGCTGCCacgcg                                                                                                                                   >  2:177682/1‑120 (MQ=255)
                                             aGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACATGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGATTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAg                                                                                                       >  7:173212/1‑149 (MQ=255)
                                               cATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGgt                                                                                                     >  4:56408/1‑149 (MQ=255)
                                                   tCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGgtcgt                                                                                                  <  3:56408/148‑1 (MQ=255)
                                                       cccGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGcc                                                                                             <  6:99336/149‑1 (MQ=255)
                                                       cccGAAACGGCGTTGAAATTGCAGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTcc                                                                                                <  2:86931/146‑1 (MQ=255)
                                                        ccGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCt                                                                                            >  6:96887/1‑149 (MQ=255)
                                                        ccGAAACGGCGTTGAAATTGCCGTTCCACACATTACAACCCGTTATGTCGGCACTTCCACCTCCAAACCATCGACGCGCAGGAACCCGCGTGGCAGCTAGCTGCCGCGCCGACCGCGCTCTCCCTTGTAGTGCTCCAGGTCGTCCCCCt                                                                                            >  3:185953/1‑149 (MQ=255)
                                                         cGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCtt                                                                                           <  5:96887/149‑1 (MQ=255)
                                                          gAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTc                                                                                          <  7:219671/149‑1 (MQ=255)
                                                           aaaCGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGCTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCa                                                                                         <  5:326414/149‑1 (MQ=255)
                                                             aCGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCaaa                                                                                       <  6:164740/149‑1 (MQ=255)
                                                                  gTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTc                                                                                          <  6:309768/141‑1 (MQ=255)
                                                                  gTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTc                                                                                          >  5:309768/1‑141 (MQ=255)
                                                                               ttCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGAAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCt                                                                                                           >  3:32335/1‑111 (MQ=255)
                                                                                   aGACATTACAACCCGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTCCGCTTGCCGGCCt                                                                 >  8:253214/1‑149 (MQ=255)
                                                                                    gACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTg                                                                >  7:239912/1‑149 (MQ=255)
                                                                                    cgacctACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTg                                                                <  3:270315/144‑1 (MQ=255)
                                                                                       aTTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCaa                                                             <  8:239912/149‑1 (MQ=255)
                                                                                       aTTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCaa                                                             >  5:209816/1‑149 (MQ=255)
                                                                                           cAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTg                                                                >  4:270315/1‑142 (MQ=255)
                                                                                            aaCCGGTTATGTCGGCACTTCCACCTGCACCCCCTCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACaa                                                        <  4:44608/149‑1 (MQ=255)
                                                                                             aCCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAg                                                       <  2:182060/149‑1 (MQ=255)
                                                                                              cgggTTATGTCGGGACTTCCAGCTGCTACCCATCGACGCGCTGGACCCCCCGTTGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTTCGCTTGCCGGCCTGCAATACAAgg                                                      <  4:185953/147‑1 (MQ=255)
                                                                                                 gTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTc                                                   >  6:200581/1‑149 (MQ=255)
                                                                                                 gTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTc                                                   <  6:187759/149‑1 (MQ=255)
                                                                                                       tCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACTAGGGTCGCGCCt                                             >  5:32506/1‑149 (MQ=255)
                                                                                                       tCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCt                                             <  6:323671/149‑1 (MQ=255)
                                                                                                        cGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCtt                                            <  7:105246/149‑1 (MQ=255)
                                                                                                          gCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCg                                          <  6:209816/149‑1 (MQ=255)
                                                                                                                  aCCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCAc                                  <  5:200581/149‑1 (MQ=255)
                                                                                                                      gCAACCCATCGACGCGCTCGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCAAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGTGTCGCGCCTTCGGGAATCACGGc                               >  6:317560/1‑148 (MQ=255)
                                                                                                                         aCCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTg                                                                >  1:36009/1‑112 (MQ=255)
                                                                                                                         aCCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTg                                                                <  2:36009/112‑1 (MQ=255)
                                                                                                                          cccATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGAt                          <  8:173212/149‑1 (MQ=255)
                                                                                                                            cATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGc                               >  3:53089/1‑142 (MQ=255)
                                                                                                                                          aaCCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAAttctt           >  1:80285/1‑148 (MQ=255)
                                                                                                                                             ccGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACg       >  2:282690/1‑149 (MQ=255)
                                                                                                                                              cgcgTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACgg      >  3:91790/1‑149 (MQ=255)
                                                                                                                                               gcgTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACGGc     >  7:170924/1‑149 (MQ=255)
                                                                                                                                                  tGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACGGCTgg  <  4:53089/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGCACCGCAAGGGGCATATCATCCGCGAAATAACCGTCCGACTCCAGCATTTCTGCCCGAAACGGCGTTGAAATTGCCGTTCCAGACATTACAACCGGTTATGTCGGCACTTCCACCTGCAACCCATCGACGCGCTGGAACCCGCGTGGCAGCTTGCTGCCGCGCCGACCGCGCTCGCCCTTGTAGTGCTCCAGGTCGTCCGCCTTCAAAGACAGGGTACGCTTGCCGGCCTGCAATACAAGGGTCGCGCCTTCGGGAATCACGGCCAGATCGGTAACAAATTCTTCACGGCTGG  >  NC_002947/5583375‑5583669

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