Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,639,580 G→A 9.6% intergenic (‑170/‑141) czcR‑III ← / → PP_1439 response regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,639,5800GA 9.6% 198.3 / 9.2 73intergenic (‑170/‑141)czcR‑III/PP_1439response regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/29);  new base A (3/4);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

GGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAAC  >  NC_002947/1639434‑1639726
                                                                                                                                                   |                                                                                                                                                   
ggACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGtt                                                                                                                                                   >  7:204887/1‑148 (MQ=255)
ggACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTa                                                                                                                                                   <  1:150913/149‑1 (MQ=255)
 gACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAg                                                                                                                                                  >  7:161886/1‑149 (MQ=255)
  aCGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAg                                                                                                                                                  >  6:61738/1‑148 (MQ=255)
    gATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTAGCGAGTAAACCCGCGCCGGTAggta                                                                                                                                               >  5:250755/1‑147 (MQ=255)
       tCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACgg                                                                                                                                            >  8:193988/1‑149 (MQ=255)
        cTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACgg                                                                                                                                            <  1:157063/148‑1 (MQ=255)
        cTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACgg                                                                                                                                            <  5:61738/148‑1 (MQ=255)
        cTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACgg                                                                                                                                            <  6:72959/148‑1 (MQ=255)
         ttGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGt                                                                                                                                          >  4:304165/1‑149 (MQ=255)
           gACCCTGTAGTGACATGAGCGGCAAGCGCCCCAATGTTACGGCCTGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACggg                                                                                                                                           <  8:267311/146‑1 (MQ=255)
             cccTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTAc                                                                                                                                      >  3:85431/1‑149 (MQ=255)
                     tGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAag                                                                                                                              <  6:146204/149‑1 (MQ=25)
                       aCATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAaga                                                                                                                             <  3:19980/148‑1 (MQ=25)
                       aCATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAagag                                                                                                                            <  3:233260/149‑1 (MQ=255)
                       aCATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAACGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAagag                                                                                                                            <  3:49378/149‑1 (MQ=255)
                        cATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAgg                                                                                                                           <  4:85431/149‑1 (MQ=255)
                          tCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGcc                                                                                                                        <  5:93384/149‑1 (MQ=255)
                               gTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGaa                                                                                                                    >  1:332955/1‑148 (MQ=255)
                                tCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGaa                                                                                                                    >  4:149529/1‑147 (MQ=255)
                                  aaGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGt                                                                                                                <  2:244264/149‑1 (MQ=255)
                                                      aTGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTg                                                                                              >  1:358066/1‑147 (MQ=255)
                                                      aTGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCa                                                                                            <  8:161886/149‑1 (MQ=255)
                                                        taaaaatGGCCTGGGGTGAATTTAACGTTTTTTTAATTCTGCTGTTATGTGGCTGTGTGTTCGGGCCTCTTCGCTAGTAACCCCGCGCCTGTAGGAACGGTTTTACCCGCTAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAg                                                                                           <  2:332955/147‑1 (MQ=255)
                                                           aaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGaa                                                                                                                               >  4:243700/1‑110 (MQ=18)
                                                           aaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGaa                                                                                                                               <  3:243700/110‑1 (MQ=18)
                                                             aaGGCCTGCGGTGCAGATAACGGATTTGTAATTATGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAAACCGCGACTGTAGGAACGGGTTTACACGCGAAGAGGCCCGAACAGTCAACACAAGGACACATTGCa                                                                                            >  5:337920/1‑142 (MQ=38)
                                                                       gTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                               <  3:321509/113‑1 (MQ=255)
                                                                       gTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                               >  4:321509/1‑113 (MQ=255)
                                                                          cAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTTATCCCCGCCTCa                                                                        >  7:194240/1‑149 (MQ=255)
                                                                           aGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCAc                                                                       >  8:292107/1‑149 (MQ=255)
                                                                               aaCGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCGCTCTTCGCGATTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcaca                                                                                                           >  1:256158/1‑110 (MQ=255)
                                                                                     tttGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGt                                                               >  3:325403/1‑147 (MQ=255)
                                                                                       tGTAATTCTGCTTTTATCTGGCTGTCTGTTCGGCCCTCTTCGCGAGAAAACCCCCGCCTGTAGCAACGGGTTAACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTTACCTTGCCGTTag                                                            <  7:193988/148‑1 (MQ=255)
                                                                                        gTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCTCTCTTCGCGAGTAAACCAGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCACGAACAGTCCAGACAAGGCCACATTGCAg                                                                                           >  1:50222/1‑117 (MQ=39)
                                                                                        gTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTAGCGAGTAAACCCGCGCATGTAGGAAAGGGTTTACACGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAAACCACGCCTCACCTTGCCGTTagc                                                           >  5:228564/1‑148 (MQ=255)
                                                                                         tAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagc                                                           >  4:322792/1‑147 (MQ=255)
                                                                                          aaTTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagcca                                                         <  4:325403/148‑1 (MQ=255)
                                                                                            ttCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCGCACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagccagcc                                                      >  7:123296/1‑149 (MQ=255)
                                                                                              ctgctgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTACAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCt                                                    >  1:307108/1‑149 (MQ=255)
                                                                                                   gTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGc                                                 >  8:348894/1‑147 (MQ=255)
                                                                                                    ttATCTGGCTGTCTGTTCGGGCCAATTAGCGCGTAAACCCGCGCCGGTTGGGACGGGGTTTACCGCGCAGAGGCCCGCACACTGCCCACAAGGCCCCAATGCACGGCGGGAATACCCCCCGCCCCGTGCCGGTAGCCACCCCGCGCTAg                                              >  2:156708/1‑149 (MQ=255)
                                                                                                          ggCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGc                                        >  3:305717/1‑149 (MQ=255)
                                                                                                           gCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAg                                              >  8:352495/1‑142 (MQ=255)
                                                                                                           gCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAg                                              <  7:352495/142‑1 (MQ=255)
                                                                                                           gCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATgg                                         <  3:304165/147‑1 (MQ=255)
                                                                                                           gCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATACCCGCCGCACCTTGCCGTTAGCCAGCCCGCGCTAGGATGGcc                                       >  6:290752/1‑149 (MQ=255)
                                                                                                               tctgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAAGTGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGGTAGCCAGCCCTCGCTAGGATGGCCgt                                     <  2:44784/147‑1 (MQ=255)
                                                                                                                 tgtTCGGGCCTCTTCGCGATTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTa                                               <  3:79817/135‑1 (MQ=255)
                                                                                                                 tgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTa                                               >  4:79817/1‑135 (MQ=255)
                                                                                                                 tgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCt                                 >  8:292971/1‑149 (MQ=255)
                                                                                                                   ttCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCTAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTc                                  <  3:149529/146‑1 (MQ=255)
                                                                                                                     cGGCCCTCTTCGCCAGTAACCCCGCGCCTGTCGGACCGGTTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAt                              <  1:156708/148‑1 (MQ=255)
                                                                                                                       ggCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCgg                           >  5:80281/1‑149 (MQ=255)
                                                                                                                        gCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGa                                                                                         <  5:28675/86‑1 (MQ=255)
                                                                                                                        gCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGa                                                                                         >  6:28675/1‑86 (MQ=255)
                                                                                                                           ctctTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGtta                        <  5:290752/149‑1 (MQ=255)
                                                                                                                            tctTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTTCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttat                       >  4:188484/1‑149 (MQ=255)
                                                                                                                            tctTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttat                       >  1:145161/1‑149 (MQ=255)
                                                                                                                                cgcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttatt                      <  7:292107/146‑1 (MQ=255)
                                                                                                                                 gcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTATGATGGCCGTGTCTGATCGGTTATTACGa                  <  5:6632/149‑1 (MQ=255)
                                                                                                                                   gggTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGatg                <  7:240470/147‑3 (MQ=37)
                                                                                                                                   gggTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAgg                >  8:240470/3‑149 (MQ=37)
                                                                                                                                   gggTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATAGGTTATTACGAgg                >  1:124917/3‑149 (MQ=37)
                                                                                                                                    aGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTCGAGGCa              <  2:358066/149‑1 (MQ=255)
                                                                                                                                     gTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGc               >  8:98871/1‑148 (MQ=255)
                                                                                                                                      tAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCa              >  3:222392/1‑148 (MQ=255)
                                                                                                                                        aaCCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAac           <  4:222392/149‑1 (MQ=255)
                                                                                                                                          cccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCGCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAaca          >  1:65968/1‑149 (MQ=255)
                                                                                                                                          cccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAacac         >  2:364790/1‑149 (MQ=255)
                                                                                                                                           ccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAacac         >  3:131294/1‑148 (MQ=255)
                                                                                                                                               gcCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGa    <  3:188484/149‑1 (MQ=255)
                                                                                                                                                 cTGTCGGACCGGTTTTACCCGCTAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAAc  <  8:194240/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAAC  >  NC_002947/1639434‑1639726

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: