Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,004,570 G→A 8.2% intergenic (+60/+90) PP_1787 → / ← PP_1788 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,004,5700GA 8.2% 241.7 / 8.3 85intergenic (+60/+90)PP_1787/PP_1788hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/37);  new base A (2/5);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-01

GCCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGTGCCATCTT  >  NC_002947/2004426‑2004715
                                                                                                                                                |                                                                                                                                                 
gcCCGTATTTTCAAGGTTTCAAAGTTCGCGTTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAAtct                                                                                                                                               <  2:75960/149‑3 (MQ=255)
gcCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAATATGAAATAGCTTATTGCAAAtct                                                                                                                                               <  2:70665/149‑3 (MQ=255)
gcCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGAttt                                                                                                                                               >  4:180435/1‑149 (MQ=255)
 cccGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTg                                                                                                                                              >  3:152624/1‑149 (MQ=255)
 cccGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTg                                                                                                                                              >  6:90260/1‑149 (MQ=255)
         ttCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAATAAg                                                                                                                                       <  3:38328/148‑1 (MQ=255)
           cAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTGCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAATAAGCt                                                                                                                                     <  2:52985/148‑1 (MQ=255)
            aaGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTAt                                                                                                                                   >  6:81039/1‑149 (MQ=255)
              ggTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTAtt                                                                                                                                  >  4:3257/1‑148 (MQ=255)
                  tCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTa                                                                                                                              >  8:360190/1‑148 (MQ=255)
                    aaaGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTaaa                                                                                                                            <  4:94481/148‑1 (MQ=255)
                      aGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAAttt                                                                                                                         <  1:80351/149‑1 (MQ=255)
                        ttCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGt                                                                                                                       <  4:152624/149‑1 (MQ=255)
                            cgGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAg                                                                                                                    >  8:360949/1‑148 (MQ=255)
                            cgGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAg                                                                                                                    <  7:360949/148‑1 (MQ=255)
                                 gAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAgg                                                                                                              >  7:217150/1‑149 (MQ=255)
                                   aaTTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGc                                                                                                             >  8:356188/1‑148 (MQ=255)
                                       tAAATGTGAGTTTTTTGACGTGACTGAATTGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATc                                                                                                         <  3:70221/148‑1 (MQ=255)
                                            gTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGtttt                                                                                                   <  7:356188/149‑1 (MQ=255)
                                            gTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGtttt                                                                                                   >  1:355131/1‑149 (MQ=255)
                                             tGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGttttt                                                                                                  >  5:236504/1‑149 (MQ=255)
                                               aGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGa                                                                                                                >  1:13967/1‑133 (MQ=255)
                                               aGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGa                                                                                                                <  2:13967/133‑1 (MQ=255)
                                               aGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGtttttt                                                                                                 <  4:348537/148‑1 (MQ=255)
                                                gTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTgg                                                                                               <  5:90260/149‑1 (MQ=255)
                                                 ttttttGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGa                                                                                              >  7:267467/1‑149 (MQ=255)
                                                 ttttttGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGa                                                                                              >  8:249583/1‑149 (MQ=255)
                                                     ttGACGTGAATGAAATGACGCTTTATGCTTAATATGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGTGGCGATCGGTTTTTTGGAtgt                                                                                           <  6:139225/148‑1 (MQ=255)
                                                       gACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAAttt                                                                                                                         >  5:223772/1‑116 (MQ=255)
                                                       gACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAAttt                                                                                                                         <  6:223772/116‑1 (MQ=255)
                                                       gACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAATAAg                                                                                                                                       <  1:215795/102‑1 (MQ=255)
                                                       gACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAATAAg                                                                                                                                       >  2:215795/1‑102 (MQ=255)
                                                               tGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCg                                                                                >  5:41456/1‑149 (MQ=255)
                                                                 aaaTGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGcc                                                                              >  3:272151/1‑149 (MQ=255)
                                                                      aCGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAg                                                                          <  5:118207/148‑1 (MQ=255)
                                                                            aatgcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTtgt                                                                   >  1:161697/1‑149 (MQ=255)
                                                                               gcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAt                                                                >  5:77933/1‑149 (MQ=255)
                                                                                cttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAt                                                                <  8:217150/148‑1 (MQ=255)
                                                                                cttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAt                                                                <  6:236504/148‑1 (MQ=255)
                                                                                cttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGATGAGGCGATCTGTTTTTTGTATGTGAACACTATCGCCTCAGGCTATGTGAtg                                                               >  6:233161/1‑149 (MQ=255)
                                                                                 ttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAtgt                                                              >  5:296354/1‑149 (MQ=255)
                                                                                  taataaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAtgtg                                                             >  2:86797/1‑149 (MQ=255)
                                                                                    ataaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGa                                                                 >  7:247453/1‑143 (MQ=255)
                                                                                    ataaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGa                                                                 <  8:247453/143‑1 (MQ=255)
                                                                                    ataaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATgtgt                                                            <  7:249583/148‑1 (MQ=255)
                                                                                    ataaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTa                                                           <  7:360190/149‑1 (MQ=255)
                                                                                      aaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTAtt                                                         >  1:345329/1‑149 (MQ=255)
                                                                                       aGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTAttt                                                        <  6:296354/149‑1 (MQ=255)
                                                                                            tttATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGa                                                    <  3:180435/148‑1 (MQ=255)
                                                                                                  gTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGtt                                             >  8:298155/1‑149 (MQ=255)
                                                                                                  gTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTAGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGtt                                             >  3:98230/1‑149 (MQ=255)
                                                                                                      aCCTCCCCTTGACATTCCCAGAACATGAAATAACTTATTGCAGAATTGCATTAAGCTATTTCTTAATTTCTTAGGTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGGCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGAt                                          <  4:134380/148‑1 (MQ=255)
                                                                                                       ccTCACCTTGACATGCCCCGAAAATGAAATTGCTTATTGCAGTTTTGAAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATGGTTTTTTTGTATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTTAt                                          <  5:233161/147‑1 (MQ=255)
                                                                                                         tCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCt                                        >  6:213661/1‑147 (MQ=255)
                                                                                                           cccTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGa                                                    <  2:364038/133‑1 (MQ=255)
                                                                                                           cccTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGa                                                    >  1:364038/1‑133 (MQ=255)
                                                                                                           cccTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTaa                                      <  6:41456/147‑1 (MQ=255)
                                                                                                             cTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCtt                                  <  2:355131/149‑1 (MQ=255)
                                                                                                                gACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTg                                 <  3:3257/147‑1 (MQ=255)
                                                                                                                  cATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCaa                              >  6:14739/1‑148 (MQ=255)
                                                                                                                    tGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATAc                           <  7:298155/149‑1 (MQ=255)
                                                                                                                       cccGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGc                                                                                   <  3:286870/90‑1 (MQ=255)
                                                                                                                       cccGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGc                                                                                   >  4:286870/1‑90 (MQ=255)
                                                                                                                       cccGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGa                        <  2:345329/149‑1 (MQ=255)
                                                                                                                        ccGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGAt                       >  1:44902/1‑149 (MQ=255)
                                                                                                                         cGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGAtt                      <  2:161697/149‑1 (MQ=255)
                                                                                                                             aaTGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGtt                  >  3:47068/1‑149 (MQ=255)
                                                                                                                             aaTGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGtt                  >  1:232237/1‑149 (MQ=255)
                                                                                                                              aTGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGttt                 <  6:77933/149‑1 (MQ=255)
                                                                                                                                gAAATAGCTTATTACAGATTTGCAATAATCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGAt                                                   <  1:10031/113‑1 (MQ=255)
                                                                                                                                gAAATAGCTTATTACAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGAt                                                   >  2:10031/1‑113 (MQ=255)
                                                                                                                                    tAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATgg           >  6:176600/1‑149 (MQ=255)
                                                                                                                                      gCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTaa              <  2:136365/144‑1 (MQ=255)
                                                                                                                                      gCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAAtaga          >  1:136365/1‑145 (MQ=255)
                                                                                                                                      gCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGt          <  5:14739/148‑1 (MQ=255)
                                                                                                                                      gCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGTg         >  5:260585/1‑149 (MQ=255)
                                                                                                                                        ttATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGAt                                                                                             >  6:291495/1‑63 (MQ=255)
                                                                                                                                        ttATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGAt                                                                                             <  5:291495/63‑1 (MQ=255)
                                                                                                                                        ttATTGCAAATCTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGTGc        >  6:330001/1‑148 (MQ=255)
                                                                                                                                           ttGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGttt                 <  3:307720/136‑1 (MQ=255)
                                                                                                                                           ttGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGttt                 >  4:307720/1‑136 (MQ=255)
                                                                                                                                           ttGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGTGCCATc    <  4:272151/149‑1 (MQ=255)
                                                                                                                                             gCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCTCATCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATATAAGCTTGCAATACGGATTGGTTTTAATTGTGCCATCtt  >  5:271784/1‑149 (MQ=255)
                                                                                                                                               aGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTCTGATGTGTATTTATGATAGTGTTGATGTAAGCTTGCAATAc                           >  5:326088/1‑122 (MQ=255)
                                                                                                                                               aGATTTGCAATAAGATATTTCTAAATTTGTTTGCTGAGGCGATCGTTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATAc                           <  6:326088/122‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GCCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGGTGCCATCTT  >  NC_002947/2004426‑2004715

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