Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707 G→C 5.4% D38E (GAC→GAG parA ← chromosome partition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,7070GC 5.4% 281.0 / 4.0 87D38E (GAC→GAGparAchromosome partition protein
Reads supporting (aligned to +/- strand):  ref base G (43/39);  new base C (3/2);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

GCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTGCGGCGC  >  NC_002947/1561‑1854
                                                                                                                                                   |                                                                                                                                                     
gCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCATACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTc                                                                                                                                                     <  4:91712/149‑1 (MQ=255)
   ggCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGAt                                                                                                                                                  <  1:332343/149‑1 (MQ=255)
        caACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATcagca                                                                                                                                             <  1:260185/149‑1 (MQ=255)
        caACTGGAAGCCTCCGTGCTCGGAGTAGTGCAAAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGAGGCGTTGCCCTGCGGATCGAGGTCGATcagca                                                                                                                                             <  7:250817/149‑1 (MQ=255)
          aCTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATCCGCAgg                                                                                                                                           <  3:248127/149‑1 (MQ=255)
            ctggaGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGac                                                                                                                                         <  1:298014/145‑1 (MQ=255)
              gAAGCCTCTGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGacac                                                                                                                                       >  4:90955/1‑149 (MQ=255)
              gAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGacac                                                                                                                                       >  6:158166/1‑149 (MQ=255)
                   cTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACacgacac                                                                                                                                  >  5:80868/1‑147 (MQ=255)
                          cTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCg                                                                                                                           >  1:363760/1‑149 (MQ=255)
                             ggAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcag                                                                                                                        >  2:221154/1‑149 (MQ=255)
                               agtagtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcagcc                                                                                                                      <  4:245107/149‑1 (MQ=255)
                                       aTAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATccgca                                                                                                                                             >  1:120860/1‑118 (MQ=255)
                                       aTAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATccgca                                                                                                                                             <  2:120860/118‑1 (MQ=255)
                                              gggCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATcagc                                                                                                                                              >  8:180791/1‑110 (MQ=255)
                                              gggCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATcagc                                                                                                                                              <  7:180791/110‑1 (MQ=255)
                                               ggCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACACTTGGTCGCAGCCAGCGAGGCGGCGAGAt                                                                                                      <  6:80868/149‑1 (MQ=255)
                                                  cAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGa                                                                                                   <  3:23013/149‑1 (MQ=255)
                                                     gTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCtt                                                                                                                                <  8:49334/117‑1 (MQ=255)
                                                     gTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCtt                                                                                                                                >  7:49334/1‑117 (MQ=255)
                                                        gCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCATGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAgg                                                                                             <  1:238702/149‑1 (MQ=255)
                                                        gCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAgg                                                                                             >  8:247165/1‑149 (MQ=255)
                                                        gCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAgg                                                                                             >  7:24219/1‑149 (MQ=255)
                                                        gCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGCTCGCTCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAgg                                                                                             >  5:363968/1‑149 (MQ=255)
                                                          aCTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCTGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtg                                                                                           >  1:303057/1‑149 (MQ=255)
                                                          aCTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtg                                                                                           >  1:4405/1‑149 (MQ=255)
                                                          aCTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtg                                                                                           <  1:221154/149‑1 (MQ=255)
                                                            tCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcag                                                                                                                        >  1:326826/1‑118 (MQ=255)
                                                            tCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcag                                                                                                                        <  2:326826/118‑1 (MQ=255)
                                                              cccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTg                                                                                     >  8:260425/1‑149 (MQ=255)
                                                                cGATGAGCAGGTCGTAGACCTATGGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCTTGGTGGCTTTGCCCTGCGTATCGCGGTCGAACAGCTGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTc                                                                                   <  2:363760/149‑1 (MQ=255)
                                                                  aTGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGAt                                                                                                                                                  >  3:348746/1‑86 (MQ=255)
                                                                  aTGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGAt                                                                                                                                                  <  4:348746/86‑1 (MQ=255)
                                                                      gCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGcc                                                                              >  5:315280/1‑148 (MQ=255)
                                                                      gCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGcc                                                                              >  3:37628/1‑148 (MQ=255)
                                                                      gCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGcc                                                                              >  2:244892/1‑148 (MQ=255)
                                                                         ggTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCcaca                                                                          >  7:119265/1‑149 (MQ=255)
                                                                             gTAGACCGAATGCTCGAGCTCGTGCTTGGTCCAAACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCAcacca                                                                       >  1:210748/1‑149 (MQ=255)
                                                                             gTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGCTTGATACAGGTGG‑‑TTGTCTTGCCCAcaccac                                                                      >  6:310237/1‑149 (MQ=255)
                                                                             gTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCAcaccac                                                                      >  4:216820/1‑149 (MQ=255)
                                                                                gACCGAATGCTCGAGCTCGTGCTTGTCCACAACGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTAGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCtt                                                                   >  4:45792/1‑149 (MQ=255)
                                                                                 aCCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCtt                                                                   >  4:366983/1‑148 (MQ=255)
                                                                                    gAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCtt                                                                                 >  6:371989/1‑131 (MQ=255)
                                                                                    gAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCtt                                                                                 <  5:371989/131‑1 (MQ=255)
                                                                                    gAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCAAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTg                                                               <  5:314986/149‑1 (MQ=255)
                                                                                       tGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGtt                                                            <  8:119265/149‑1 (MQ=255)
                                                                                         cTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTcg                                                          >  5:204002/1‑149 (MQ=255)
                                                                                          tCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCtt                                                                                 <  4:269892/125‑1 (MQ=255)
                                                                                          tCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCtt                                                                                 >  3:269892/1‑125 (MQ=255)
                                                                                            gAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGa                                                       <  7:260425/149‑1 (MQ=255)
                                                                                                    gCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGaa                                                >  4:290091/1‑148 (MQ=255)
                                                                                                         gTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATAcc                                            <  6:204002/148‑1 (MQ=255)
                                                                                                          tCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCt                                           >  3:144924/1‑148 (MQ=255)
                                                                                                           ccACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCt                                           >  6:334122/1‑147 (MQ=255)
                                                                                                            cacaCCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCGTTCTGGTTCGCGATTGCGAATACCTTAg                                        >  5:223644/1‑149 (MQ=255)
                                                                                                             acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGc                                       >  6:300634/1‑149 (MQ=255)
                                                                                                             acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGc                                       <  3:290091/149‑1 (MQ=255)
                                                                                                             acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGc                                       >  3:349872/1‑149 (MQ=255)
                                                                                                             acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGc                                       <  3:366983/149‑1 (MQ=255)
                                                                                                             acacCGCTGCCCATGGTGGCGTTGCCCTGCGCATCGAGGTCGATCAGCAGGCCACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATAAACGTGG‑‑ATGTCTTGCCAACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGc                                       >  5:188043/1‑149 (MQ=255)
                                                                                                               acCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCa                                     >  5:42534/1‑149 (MQ=255)
                                                                                                                     gCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAtt                                   >  4:62188/1‑145 (MQ=255)
                                                                                                                     gCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAtt                                   <  3:62188/145‑1 (MQ=255)
                                                                                                                      cccATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGc                              <  4:349872/149‑1 (MQ=255)
                                                                                                                        cATGTTGGCGTTTCCCTGCGGTTGTAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGATGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACGTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGa                            <  6:188043/149‑1 (MQ=255)
                                                                                                                        cATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGa                            <  1:244892/149‑1 (MQ=255)
                                                                                                                        cATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGa                            <  5:300634/149‑1 (MQ=255)
                                                                                                                        cATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGa                            <  2:303057/149‑1 (MQ=255)
                                                                                                                        cATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGa                            <  7:247165/149‑1 (MQ=255)
                                                                                                                        cATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGa                            <  2:4405/149‑1 (MQ=255)
                                                                                                                         aTGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAg                           >  2:172013/1‑149 (MQ=255)
                                                                                                                              ggCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTc                      >  4:68366/1‑149 (MQ=255)
                                                                                                                              ggCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTc                      >  8:254897/1‑149 (MQ=255)
                                                                                                                               gCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTc                      <  2:210748/148‑1 (MQ=255)
                                                                                                                                 gTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTTCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCa                   <  1:279731/149‑1 (MQ=255)
                                                                                                                                 gTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCt                                           <  8:144709/125‑1 (MQ=255)
                                                                                                                                 gTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAttgt                                 <  8:69468/135‑1 (MQ=255)
                                                                                                                                 gTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAttgt                                 >  7:69468/1‑135 (MQ=255)
                                                                                                                                  ttGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAg                  <  6:223644/149‑1 (MQ=255)
                                                                                                                                     tccTGCTGATCGACCTCGATCCGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCa               >  2:83653/2‑149 (MQ=255)
                                                                                                                                     cccTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCa               >  6:362945/1‑149 (MQ=255)
                                                                                                                                      ccTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTa                                         >  3:18571/1‑122 (MQ=255)
                                                                                                                                      ccTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTa                                         <  4:18571/122‑1 (MQ=255)
                                                                                                                                        tGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGAAAAGGCTTGCCCACATAACCTATATGGTTAGCGATTGAGAATACATTAGCCATTGTTGCGAGTGTTACCAgaaat              >  1:114941/1‑145 (MQ=255)
                                                                                                                                         gCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGcc           >  7:54356/1‑149 (MQ=255)
                                                                                                                                           ggATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCtt         <  5:334122/149‑1 (MQ=255)
                                                                                                                                              tCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTgcg      <  3:312960/149‑1 (MQ=255)
                                                                                                                                                  ggTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTgcg      <  2:3137/145‑1 (MQ=255)
                                                                                                                                                  ggTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTgcg      >  1:3137/1‑145 (MQ=255)
                                                                                                                                                  ggTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTGCGgcgc  <  3:68366/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGG‑‑TTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTGCGGCGC  >  NC_002947/1561‑1854

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: