Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,037,575 A→T 8.7% intergenic (+169/+30) PP_1810 → / ← PP_5478 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,037,5750AT 8.7% 218.7 / 8.5 77intergenic (+169/+30)PP_1810/PP_5478hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/32);  new base T (2/5);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

TTGTAGGCATGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACC  >  NC_002947/2037429‑2037712
                                                                                                                                                  |                                                                                                                                          
ttGTAGGCATGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCa                                                                                                                                            >  3:185636/1‑147 (MQ=255)
 tGTAGTCATGGGATGGATTGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCGCGCTGCCACTGTGTGGGGTGCTTGGCGAATGCAGTGGTTTGCCCTCGCCACGTCTGGCGTGATGAAGCCGATTTCTTCTGACTGCAATCGAGGCATa                                                                                                                                          <  8:130369/148‑1 (MQ=255)
   tAGGCATGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGc                                                                                                                                        <  2:102264/148‑1 (MQ=255)
     ggCATGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATCTCTTCTGACTGCAAGCGAGGCATAGCCTc                                                                                                                                     >  3:295483/1‑149 (MQ=255)
        aTGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCACCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGcc                                                                                                                                  >  1:334203/1‑149 (MQ=255)
         tGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAACGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCt                                                                                                                                 >  7:119774/1‑149 (MQ=255)
            gatggatgGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCTATGCCTCGCTTgca                                                                                                                              <  6:118738/149‑3 (MQ=255)
             atggatggCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCTATGCCTCGCTTgcag                                                                                                                             <  6:307621/149‑4 (MQ=39)
              tggatggCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCTATGCCTCGCTTgcagt                                                                                                                            <  7:73678/149‑5 (MQ=39)
                gatggCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATCTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGtctc                                                                                                                           <  4:295483/148‑1 (MQ=255)
                     cACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCa                                                                                                                    >  4:63882/1‑149 (MQ=255)
                     cACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCa                                                                                                                    >  5:154856/1‑149 (MQ=255)
                     cACCGGCGGTGCCGGTGTTCGCGGGCGCCCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCa                                                                                                                    >  1:205894/1‑149 (MQ=255)
                          gcggTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTccc                                                                                                               >  1:316810/1‑149 (MQ=255)
                           cggTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCt                                                                                                              >  3:272887/1‑149 (MQ=255)
                            ggTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCtt                                                                                                             <  1:23112/149‑1 (MQ=255)
                               gCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAg                                                                                                          <  4:272887/149‑1 (MQ=255)
                                         gcgGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTcgc                                                                                                                         >  4:228887/1‑125 (MQ=255)
                                              cgcgcCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTc                                                                                           <  7:150722/149‑1 (MQ=255)
                                              cgcgcCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTc                                                                                           <  4:185636/149‑1 (MQ=255)
                                              cgcgcCCGCTCCCACAGTGTGGGGTGCTTGGCCAACGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTc                                                                                           <  8:119774/149‑1 (MQ=255)
                                               gcgcCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCATCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCTATGCATAGCCTCTCATGGTCTCGCTTCATTCCCTTCATTGGTGCCTGTTTTTc                                                                                           >  6:332810/1‑148 (MQ=255)
                                                         ccACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGc                                                                                                                                   <  3:356014/99‑1 (MQ=255)
                                                         ccACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGc                                                                                                                                   >  4:356014/1‑99 (MQ=255)
                                                         ccACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTc                                                                                           <  8:166795/138‑1 (MQ=255)
                                                         ccACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTc                                                                                           >  7:166795/1‑138 (MQ=255)
                                                         ccACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGAt                                                                                >  5:287319/1‑149 (MQ=255)
                                                           acaGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTc                                                                                           <  8:153469/136‑1 (MQ=255)
                                                           acaGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTc                                                                                           >  7:153469/1‑136 (MQ=255)
                                                               tgtgGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAAtc                                                                                     <  7:279586/138‑1 (MQ=255)
                                                               tgtgGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAAtc                                                                                     >  8:279586/1‑138 (MQ=255)
                                                                  ggtgTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGTCCTGATGAAGGCCATTTCTTCTGTCTGCTAGCGAGGCATAGACTCTCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAAtctc                                                                                   <  5:45483/134‑1 (MQ=255)
                                                                  ggggTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGTTCTCGCGTCATTCCCTACAGTGGTGCCTGTTTTTAGCAAtctc                                                                                   >  6:45483/1‑137 (MQ=255)
                                                                    ggTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTtcatca                                                                     <  3:63882/149‑1 (MQ=255)
                                                                     gTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAg                                                                    <  1:329164/149‑1 (MQ=255)
                                                                        cTTGGCCAATGCAGTGGTTTGCCCTCGCCCCGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGc                                                                 >  4:41385/1‑149 (MQ=255)
                                                                            gCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCtttt                                                             >  8:17075/1‑149 (MQ=255)
                                                                            gCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCtttt                                                             <  6:287319/149‑1 (MQ=255)
                                                                             ccAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTg                                                            <  2:308264/149‑1 (MQ=255)
                                                                              cAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGc                                                           >  5:58481/1‑149 (MQ=255)
                                                                                    aGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGc                                                     <  2:316810/149‑1 (MQ=255)
                                                                                    aGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGc                                                     <  7:17075/149‑1 (MQ=255)
                                                                                      tGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGtt                                                                          >  6:52045/1‑126 (MQ=255)
                                                                                      tGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGtt                                                                          <  5:52045/126‑1 (MQ=255)
                                                                                        gTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCcag                                                  >  3:81712/1‑148 (MQ=255)
                                                                                         tttGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAAtc                                                                                     <  7:90482/112‑1 (MQ=255)
                                                                                         tttGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAAtc                                                                                     >  8:90482/1‑112 (MQ=255)
                                                                                              ccTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTaa                                           >  6:23889/1‑149 (MQ=255)
                                                                                                tCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATc                                         >  5:227870/1‑149 (MQ=255)
                                                                                                     aCGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTg                                    <  6:154856/149‑1 (MQ=255)
                                                                                                            gCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCa                              <  4:81712/148‑1 (MQ=255)
                                                                                                             catgatgaATCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTTTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAg                             <  5:332810/146‑1 (MQ=255)
                                                                                                                  tgaAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCtt                       >  2:50560/1‑149 (MQ=255)
                                                                                                                      gCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTgg                   >  7:303687/1‑149 (MQ=255)
                                                                                                                         cATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATc                                         >  5:320014/1‑124 (MQ=255)
                                                                                                                         cATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATc                                         <  6:320014/124‑1 (MQ=255)
                                                                                                                         cATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtcg                >  2:204950/1‑149 (MQ=255)
                                                                                                                         cATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtcg                >  4:357878/1‑149 (MQ=255)
                                                                                                                         cATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTACCTTCAGTGGTGCCAGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATACGCTCCGCCCAGCATCTTGTGgtcg                >  6:202475/1‑149 (MQ=255)
                                                                                                                           tttCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtcgt               >  1:73653/1‑149 (MQ=255)
                                                                                                                            ttcttcTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtcgtc              <  4:139075/148‑1 (MQ=255)
                                                                                                                            ttcttcTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtcgtc              >  3:139075/1‑148 (MQ=255)
                                                                                                                               ttctGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTc           >  5:100822/1‑148 (MQ=255)
                                                                                                                               ttctGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTc           >  7:18971/1‑148 (MQ=255)
                                                                                                                                 ctGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCCGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTcc          <  2:205894/147‑1 (MQ=255)
                                                                                                                                   gACTGCAAGCGAGGCATAGCCTCGCCTGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGt      >  5:243365/1‑148 (MQ=255)
                                                                                                                                    aCTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCAt        <  2:73653/146‑1 (MQ=255)
                                                                                                                                    aCTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGGCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCAt        <  2:334203/146‑1 (MQ=255)
                                                                                                                                     cTGCAAGCGATGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTc     <  6:58481/148‑1 (MQ=255)
                                                                                                                                      tGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGGTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCAc   >  7:191667/1‑149 (MQ=255)
                                                                                                                                      tGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtc                 >  6:297048/1‑135 (MQ=255)
                                                                                                                                      tGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtc                 <  5:297048/135‑1 (MQ=255)
                                                                                                                                      tGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCAc   >  3:192161/1‑149 (MQ=255)
                                                                                                                                        cAATCGAGGCATAGCCTCTCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGTTAATCCGCTGCTCCCAGGATCTTGTGGTCGTCTTCCATGTCAcc  <  5:202475/148‑1 (MQ=255)
                                                                                                                                         aggCGAGGCTATGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTtgtg                    >  5:282547/3‑129 (MQ=255)
                                                                                                                                         aggCGAGGCTATGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTtgtg                    <  6:282547/127‑1 (MQ=255)
                                                                                                                                         aggCGAGGCTATGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtcgt               >  3:186802/3‑134 (MQ=255)
                                                                                                                                         aggCGAGGCTATGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGgtcgt               <  4:186802/132‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
TTGTAGGCATGGGATGGATGGCACCGGCGGTGCCGGTGTTCGCGGGCGCGCCCGCTCCCACAGTGTGGGGTGCTTGGCCAATGCAGTGGTTTGCCCTCGCCACGTCTGGCCTGATGAAGCCCATTTCTTCTGACTGCAAGCGAGGCATAGCCTCGCCTGGTCTCGCGTCATTCCCTTCAGTGGTGCCTGTTTTTCGCAATCTCGATGTCGTTCATCAGGGCTTTTGCGGTGGCCAGCAGGTAATCCGCTGCGCCCAGCATCTTGTGGTCGTCTTCCATGTCACC  >  NC_002947/2037429‑2037712

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: