Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,357,746 T→C 12.1% intergenic (+92/‑162) PP_3831 → / → PP_3832 hypothetical protein/CsrA‑like carbon storage regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,357,7460TC12.1% 164.8 / 12.4 68intergenic (+92/‑162)PP_3831/PP_3832hypothetical protein/CsrA‑like carbon storage regulator
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base C (5/3);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

AAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCGGACTTATACATC  >  NC_002947/4357602‑4357892
                                                                                                                                                |                                                                                                                                                  
aaGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACg                                                                                                                                                <  1:218294/149‑1 (MQ=255)
 aGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGATGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACg                                                                                                                                                <  2:107865/148‑1 (MQ=255)
    ggCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCTACAAACCGTTGACGGttt                                                                                                                                            >  5:23048/1‑149 (MQ=255)
      ccTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTCAACGGTTTgt                                                                                                                                          <  1:296443/149‑2 (MQ=255)
       cTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCt                                                                                                                                         <  1:212207/149‑1 (MQ=255)
           cTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCGTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTgg                                                                                                                                     <  6:384857/149‑1 (MQ=255)
           cTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTgg                                                                                                                                     >  5:266385/1‑149 (MQ=255)
              ttgttTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCAc                                                                                                                                  <  7:147162/149‑1 (MQ=255)
              ttgttTCTTCAGGCACTTGAGGAACAGGAGGACCATAGAGCTTGTCTTTTCAGCCCATGATGCGAATAAATTTCCCATTAGCCGCCAACCCGCCTATTCTTGCCAGCGAGCACCTTGGCCGACATACCTTTGACGGTTTGCTTTGGc                                                                                                                                    <  8:95154/147‑1 (MQ=255)
                gttTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCt                                                                                                                                >  3:323301/1‑149 (MQ=255)
                     ttcAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAAc                                                                                                                           >  5:419559/1‑149 (MQ=255)
                         ggCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGaa                                                                                                                       <  6:266385/149‑1 (MQ=255)
                           cACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGa                                                                                                                     <  8:280984/149‑1 (MQ=255)
                           cACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGa                                                                                                                     <  6:419559/149‑1 (MQ=255)
                            aCTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGa                                                                                                                     >  4:172287/1‑148 (MQ=255)
                              ttGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGacac                                                                                                                  >  7:232746/1‑149 (MQ=255)
                              ttGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGacac                                                                                                                  >  3:395006/1‑149 (MQ=255)
                                gAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACg                                                                                                                 <  8:348245/148‑1 (MQ=255)
                                        ggaggaCCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAAt                                                                                                        <  4:395006/149‑1 (MQ=255)
                                           ggaCCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTtg                                                                                                     >  5:359344/1‑149 (MQ=255)
                                            gaCCATTGAGCTTGTCTATTGAGACCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGTCACCTGCAACTg                                                                                                                         <  6:166446/128‑1 (MQ=255)
                                            gaCCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTtgt                                                                                                    <  4:323301/149‑1 (MQ=255)
                                            gaCCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTCAACGGTTTgtcg                                                                                                                                        >  3:248523/1‑110 (MQ=255)
                                            gaCCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTCAACGGTTTgtcg                                                                                                                                        <  4:248523/113‑4 (MQ=255)
                                            gaCAATTGAGCTTGTCTATTCATCCCATGATGCGAAGAAATTTCCCATTAGACGCCAACCCCACTATTCTTGCCAGCGAGCAACTTGGCAGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTg                                                                                                                         >  5:166446/1‑128 (MQ=255)
                                                aTTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTtgtctgt                                                                                                >  4:67025/1‑149 (MQ=255)
                                                aTTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTCGACGGTTTGCTTTgg                                                                                                                                     >  1:195744/1‑112 (MQ=255)
                                                aTTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTCGACGGTTTGCTTTgg                                                                                                                                     <  2:195744/112‑1 (MQ=255)
                                                 ttGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTtgtctgtc                                                                                               <  8:232746/149‑1 (MQ=255)
                                                     gCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCggg                                                                                            >  4:137888/1‑148 (MQ=255)
                                                     gCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGACAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCggg                                                                                            <  4:313799/148‑1 (MQ=255)
                                                         gTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAAt                                                                                        >  7:41248/1‑148 (MQ=255)
                                                          tCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTa                                                                                      <  8:41248/149‑1 (MQ=255)
                                                            tATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGaa                                                                                         <  2:278922/144‑1 (MQ=255)
                                                            tATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGaa                                                                                         >  1:278922/1‑144 (MQ=255)
                                                             aTTCAGCACATGATGAGAAGAAATTTCCCATTAGCCGCCAAACCCCCTATTCTTGCCAGCTAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCGCATGCAACTGAAGACACGAAAAAGAATTTGTGTGTGTGGTAATTACTa                                                                                   >  6:344692/1‑149 (MQ=255)
                                                                 aGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTgg                                                                                                                                     >  7:19145/1‑95 (MQ=255)
                                                                 aGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTgg                                                                                                                                     <  8:19145/95‑1 (MQ=255)
                                                                  gCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTATTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATccgc                                                                              >  6:371808/1‑149 (MQ=255)
                                                                        gatgCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGGCAGCGAGCACCTTGGCCGACAAACCGTTGTCGGTTTGCTTTGGCACCTGGAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGaaa                                                                         <  7:82933/148‑1 (MQ=255)
                                                                         atgCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGTCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGc                                                                       <  7:214735/149‑1 (MQ=255)
                                                                            cGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGt                                                                     <  3:309254/148‑1 (MQ=255)
                                                                              aagaaATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCtt                                                                  <  1:158823/149‑1 (MQ=255)
                                                                                  aaTTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGATTTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAAt                                                              <  6:23048/149‑1 (MQ=255)
                                                                                     ttCCCATTACCCGCCAACCCCCCTAGTCATTCCAGCGAGCACCTTGGCCGACAAACCGTAGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAACGCGTCTTTAAt                                                              >  1:22707/1‑146 (MQ=255)
                                                                                         cATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAgaagg                                                       <  3:67025/149‑1 (MQ=255)
                                                                                         cATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAgaag                                                        >  1:346070/1‑148 (MQ=255)
                                                                                         cATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAgaag                                                        >  4:130569/1‑148 (MQ=255)
                                                                                         cATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCAGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAgaagg                                                       <  6:192495/149‑1 (MQ=255)
                                                                                          aTTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGATTTACTATCCGCCGAAAGCGTCTTTAATGAgaagg                                                       >  5:221280/1‑148 (MQ=255)
                                                                                              gccgccAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAagag                                                  >  4:122956/1‑149 (MQ=255)
                                                                                               ccgccAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAgg                                                 >  6:188335/1‑149 (MQ=255)
                                                                                               ccgccAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAgg                                                 <  2:346070/149‑1 (MQ=255)
                                                                                                cgccAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAggg                                                <  3:122956/149‑1 (MQ=255)
                                                                                                       cccccTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCt                                         >  3:8805/1‑149 (MQ=255)
                                                                                                         cccTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTcc                                       >  3:21625/1‑149 (MQ=255)
                                                                                                                      aGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGc                           >  6:23248/1‑148 (MQ=255)
                                                                                                                       gCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGATTTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCaa                         <  6:221280/149‑1 (MQ=255)
                                                                                                                           gCACCTTGGCCTACAAACCGTTGACGGATTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAg                     <  5:344692/149‑1 (MQ=255)
                                                                                                                                  ggCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGc               <  3:172287/148‑1 (MQ=255)
                                                                                                                                  ggCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGc               <  4:8805/148‑1 (MQ=255)
                                                                                                                                  ggCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCg              >  1:384865/1‑149 (MQ=255)
                                                                                                                                   caaagcAAACCGTCAACGGTTTGTTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCgg             >  7:255117/6‑149 (MQ=255)
                                                                                                                                   caaagcAAACCGTCAACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCg              >  6:154513/6‑148 (MQ=255)
                                                                                                                                    ccGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCGGa            >  6:413837/1‑149 (MQ=255)
                                                                                                                                     cGACAAACCGTGGACGGTTTGCTTTGGCACCTTCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCGGAc           <  8:237554/149‑1 (MQ=255)
                                                                                                                                      agcAAACCGTCAACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGATAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCGGACt          >  6:120361/3‑149 (MQ=255)
                                                                                                                                              gTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCGGACTTATACATc  <  6:359344/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAGGAGGACCATTGAGCTTGTCTATTCAGCCCATGATGCGAAGAAATTTCCCATTAGCCGCCAACCCCCCTATTCTTGCCAGCGAGCACCTTGGCCGACAAACCGTTGACGGTTTGCTTTGGCACCTGCAACTGAAGACACGCAACAGAATTTGTCTGTCGGGGAATTACTATCCGCCGAAAGCGTCTTTAATGAGAAGGAAGAGGGACAAGCTCCCGCCGCCGCGGCAACTAGCCCAGCGGACTTATACATC  >  NC_002947/4357602‑4357892

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