Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,971,398 T→C 6.0% L579L (CTT→CTC)  PP_2599 → chemotaxis sensory transducer family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,971,3980TC 6.0% 195.4 / 2.5 66L579L (CTT→CTC) PP_2599chemotaxis sensory transducer family protein
Reads supporting (aligned to +/- strand):  ref base T (30/32);  new base C (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

ACGACTTCGCTTTCGTCGCCCAGGCCAACATCGGCGGCGAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGGACGACCCGCGTAACATCTGGGGGTTCGACATGCATCTGCCAGGC  >  NC_002947/2971250‑2971545
                                                                                                                                                    |                                                                                                                                                   
acgacTTCGCTTTCGTCGCCCAGGCCAACATCGGCGGCGAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCtt                                                                                                                                                     <  5:84557/149‑1 (MQ=255)
     ttCGCTTTCGTCGCCCAGGCCAACATCGGCGGCGAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCgg                                                                                                                                                >  1:81146/1‑149 (MQ=255)
       cGCTTTCGTCGCCCAGGCCAACATCGGCGGCGAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTc                                                                                                                                                  <  6:99416/145‑1 (MQ=255)
       cGCTTTCGTCGCCCAGGCCAACATCGGCGGCGAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTc                                                                                                                                                  >  5:99416/1‑145 (MQ=255)
         cTTTCGTCGCCCAGGCCAACATCGGCGGCGAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTg                                                                                                                                            <  2:170860/149‑1 (MQ=255)
                   ccAGGCCAACATCGGCGGCGAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACAAGGCGTTCCGCCTCGTcac                                                                                                                                                <  3:63547/135‑3 (MQ=255)
                   ccAGGCCAACATCGGCGGCGAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACAAGGCGTTCCGCCTCGTcac                                                                                                                                                >  4:63547/1‑133 (MQ=255)
                          aaCATCGGCGGCGAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTAcc                                                                                                                           >  4:179404/1‑149 (MQ=255)
                                  cggcgAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGAc                                                                                                                        <  6:19556/144‑1 (MQ=255)
                                  cggcgAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGAc                                                                                                                        >  5:19556/1‑144 (MQ=255)
                                  cggcgAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAgagg                                                                                                                   <  3:104432/149‑1 (MQ=255)
                                   ggcgAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAgagga                                                                                                                  <  7:1663/149‑1 (MQ=255)
                                        gTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTAcc                                                                                                                           <  3:272484/135‑1 (MQ=255)
                                        gTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTAcc                                                                                                                           >  4:272484/1‑135 (MQ=255)
                                            aTGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGc                                                                                                          <  3:179404/148‑1 (MQ=255)
                                            aTGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCt                                                                                                         >  6:133960/1‑149 (MQ=255)
                                             tGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTg                                                                                                        >  4:61067/1‑149 (MQ=255)
                                             tGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTg                                                                                                        <  8:221183/149‑1 (MQ=255)
                                                aGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCACCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGcc                                                                                                     <  8:291550/149‑1 (MQ=255)
                                                       cGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGcc                                                                                              >  8:175884/1‑149 (MQ=255)
                                                       cGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGcc                                                                                              >  1:146991/1‑149 (MQ=255)
                                                         atatCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGcccc                                                                                            <  8:203934/149‑1 (MQ=255)
                                                         atatCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACAAGGCGTTCCGCCTCGTcacg                                                                                                                                               <  3:277937/98‑4 (MQ=255)
                                                         atatCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACAAGGCGTTCCGCCTCGTcacg                                                                                                                                               >  4:277937/1‑95 (MQ=255)
                                                              cgtcgcTTGCCGATCACCTGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAg                                                                                                                                      <  1:122744/99‑1 (MQ=38)
                                                              cgtcgcTTGCCGATCACCTGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAg                                                                                                                                      >  2:122744/4‑102 (MQ=38)
                                                              cAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCa                                                                                       >  4:52724/1‑149 (MQ=255)
                                                               gtcgcTTGCCGATCACCTGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCAc                                                                                      >  2:186122/3‑149 (MQ=255)
                                                                 cgcTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTg                                                                                    <  7:175884/149‑1 (MQ=255)
                                                                      gCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTg                                                                                <  3:52724/148‑1 (MQ=255)
                                                                         ggTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCt                                                                            >  8:266289/1‑149 (MQ=255)
                                                                         ggTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCt                                                                            >  3:2034/1‑149 (MQ=255)
                                                                           tGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTg                                                                                <  1:117136/143‑1 (MQ=255)
                                                                           tGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTg                                                                                >  2:117136/1‑143 (MQ=255)
                                                                           tGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTgg                                                                          <  2:81146/149‑1 (MQ=255)
                                                                            gATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTggg                                                                         >  5:291309/1‑149 (MQ=255)
                                                                            gATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTggg                                                                         >  4:269726/1‑149 (MQ=255)
                                                                             acctgCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTggg                                                                         <  1:186122/144‑1 (MQ=255)
                                                                             aTCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTggg                                                                         <  4:2034/148‑1 (MQ=255)
                                                                               cGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCa                                                                      >  8:160340/1‑149 (MQ=255)
                                                                                ggCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGcg                                                                                                         >  5:308515/1‑112 (MQ=255)
                                                                                ggCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCt                                                                                                         <  6:308515/113‑1 (MQ=255)
                                                                                ggCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTgg                                                                          <  5:240526/144‑1 (MQ=255)
                                                                                ggCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGAAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTgg                                                                          >  6:240526/1‑144 (MQ=255)
                                                                                 gCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCAcc                                                                    <  5:244128/149‑1 (MQ=255)
                                                                                    aGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGc                                                                 <  7:160340/149‑1 (MQ=255)
                                                                                        aCGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCAcgc                                                             <  7:266289/149‑1 (MQ=255)
                                                                                           ctggctgCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCAcgccgc                                                          <  5:133960/149‑1 (MQ=255)
                                                                                               ctgCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGccc                                                                                             >  6:190331/1‑110 (MQ=255)
                                                                                               ctgCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGccc                                                                                             <  5:190331/110‑1 (MQ=255)
                                                                                                      aTTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGa                                                                                                                         <  2:286778/75‑1 (MQ=255)
                                                                                                      aTTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGa                                                                                                                         >  1:286778/1‑75 (MQ=255)
                                                                                                      aTTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGgtggtg                                                  >  5:57119/1‑146 (MQ=255)
                                                                                                      aTTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGaa                                                <  6:57119/148‑1 (MQ=255)
                                                                                                         tCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGga                                             <  4:24400/148‑1 (MQ=255)
                                                                                                            aTAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGgacgac                                         >  1:267877/1‑149 (MQ=255)
                                                                                                               aCCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAg                                               <  8:125397/140‑1 (MQ=255)
                                                                                                               aCCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGGGGAGAAg                                               >  7:125397/1‑140 (MQ=255)
                                                                                                                      gggCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGGACGACCCGCGTAAc                                >  2:38915/1‑148 (MQ=255)
                                                                                                                               ccGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGGACGACCCGCGTAACATCTgg                          <  3:61067/145‑1 (MQ=255)
                                                                                                                               ccGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGAACTGCGTGCCCTGGGGCACCCCCACGCCGCCGGTGGTGAAGGACGACCCGCGTAACATCTggggct                      >  4:36932/1‑147 (MQ=255)
                                                                                                                                 gTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGGACGACCCGCGTAACATCTGGGGGTTc                    >  3:250604/1‑149 (MQ=255)
                                                                                                                                   gACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGa                                                                                          <  7:217252/77‑1 (MQ=255)
                                                                                                                                   gACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGa                                                                                          >  8:217252/1‑77 (MQ=255)
                                                                                                                                    aCGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGGACGACCCGCGTAACATCTGGGGGTTCGAc                 <  8:171949/149‑1 (MQ=255)
                                                                                                                                             aaCGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGGACGACCCGCGTAACATCTGGGGGTTCGACATGCATCTg        <  1:245145/149‑1 (MQ=255)
                                                                                                                                                gCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGGACGACCCGCGTAACATCTGGGGGTTCGACATGCATCTGc       >  7:153498/1‑147 (MQ=255)
                                                                                                                                                   ttGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGGACGACCCGCGTAACATCTGGTGGTTCGACATGCATCTGCCAGGc  >  4:46715/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGACTTCGCTTTCGTCGCCCAGGCCAACATCGGCGGCGAGTTCATGCAGGACACCGATATCCAGCGCTTGCAGGTGATCGGCAAGCGACGCTGGCTGCGCCATTTCGATAACCGCCTGGGCATCTCCCGTGACGAGGCGGAACGCCTTGTCGGCGTGCCGCAGCCTGAACTACCGACAGAGGAGCCATTGCTGGCCGACCGCCCCGAGCACTGCGTGCCCTGGGGCACCCGCACGCCGCCGGTGGTGAAGGACGACCCGCGTAACATCTGGGGGTTCGACATGCATCTGCCAGGC  >  NC_002947/2971250‑2971545

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: ‑