Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,280,219 A→C 6.2% E80D (GAA→GAC yjcE → putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,280,2190AC 6.2% 325.9 / 13.9 97E80D (GAA→GACyjcEputative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (40/50);  new base C (3/3);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

GCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGG  >  NC_000913/4280072‑4280367
                                                                                                                                                   |                                                                                                                                                    
gcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAAc                                                                                                                                                     >  3:666829/1‑149 (MQ=255)
     tGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGt                                                                                                                                                >  3:270088/1‑149 (MQ=255)
     tGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGt                                                                                                                                                >  1:268796/1‑149 (MQ=255)
          aTCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCgaga                                                                                                                                           <  3:111425/149‑1 (MQ=255)
            cGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGACCATGTTCgagaaa                                                                                                                                         <  1:718443/149‑3 (MQ=255)
                                   ggCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTc                                                                                                                                    >  4:904097/1‑131 (MQ=255)
                                   ggCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTc                                                                                                                                    <  3:904097/131‑1 (MQ=255)
                                    gCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGcgctggcgctgg                                                                                                                 <  2:123873/149‑1 (MQ=255)
                                      cGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtg                                                                                                               >  4:450979/1‑149 (MQ=255)
                                      cGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtg                                                                                                               <  4:1142102/149‑1 (MQ=255)
                                        aCGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCATGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggt                                                                                                             <  1:103498/149‑1 (MQ=255)
                                         cGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtg                                                                                                            <  3:1106038/149‑1 (MQ=255)
                                               gtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCAcc                                                                                                      >  3:150659/1‑149 (MQ=255)
                                               gtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCAcc                                                                                                      >  2:684191/1‑149 (MQ=255)
                                               gtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCAcc                                                                                                      >  4:493759/1‑149 (MQ=255)
                                                tttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCg                                                                                                     <  4:371954/149‑1 (MQ=255)
                                                tttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCg                                                                                                     >  4:471259/1‑149 (MQ=255)
                                                tttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCg                                                                                                     >  4:64953/1‑149 (MQ=255)
                                                   gCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTgg                                                                                                  <  2:216374/149‑1 (MQ=255)
                                                     aTGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGACCATGTTCGAGAAAttcac                                                                                                                                    >  2:902591/1‑110 (MQ=255)
                                                     aTGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGACCATGTTCGAGAAAttcac                                                                                                                                    <  1:902591/113‑4 (MQ=255)
                                                     aTGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTc                                                                                                <  2:1008998/149‑1 (MQ=255)
                                                       gtgGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCgg                                                                                              >  2:45903/1‑149 (MQ=255)
                                                          gAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCAt                                                                                           <  1:71495/149‑1 (MQ=255)
                                                            gTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTg                                                                                         <  1:45903/149‑1 (MQ=255)
                                                             tttGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTgg                                                                                        >  1:103150/1‑149 (MQ=255)
                                                             tttGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTgg                                                                                        >  2:1128418/1‑149 (MQ=255)
                                                              ttGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGc                                                                                       <  1:1128418/149‑1 (MQ=255)
                                                              ttGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGc                                                                                       <  1:739129/149‑1 (MQ=255)
                                                                 aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTc                                                                                    >  1:338685/1‑149 (MQ=255)
                                                                  tCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                   <  3:493759/149‑1 (MQ=255)
                                                                  tCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                   >  1:169949/1‑149 (MQ=255)
                                                                   ccTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCt                                                                                  <  1:684191/149‑1 (MQ=255)
                                                                    cTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCtt                                                                                 >  3:604339/1‑149 (MQ=255)
                                                                     tGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTa                                                                                <  1:222359/149‑1 (MQ=255)
                                                                     tGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTa                                                                                >  2:177321/1‑149 (MQ=255)
                                                                         ctctTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTa                                                                            <  2:277367/149‑1 (MQ=255)
                                                                         ctctTTTTAGTCTTGTTTATCCCGCCATTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAGCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTTCTTATTTa                                                                            >  3:961086/1‑149 (MQ=255)
                                                                            tttttAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTg                                                                         <  2:1037906/149‑1 (MQ=255)
                                                                             ttttAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTgg                                                                        <  2:169949/149‑1 (MQ=255)
                                                                                 aGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGggtgg                                                                    >  3:1038499/1‑149 (MQ=255)
                                                                                   tCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTgg                                                                                                  >  2:443966/1‑117 (MQ=255)
                                                                                   tCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTAt                                                                               <  1:324463/136‑1 (MQ=255)
                                                                                   tCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTgg                                                                        <  2:1016317/143‑1 (MQ=255)
                                                                                   tCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTgg                                                                        >  1:1016317/1‑143 (MQ=255)
                                                                                   tCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGgtggtg                                                                  >  4:379218/1‑149 (MQ=255)
                                                                                   tCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCACTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTAt                                                                               >  2:324463/1‑136 (MQ=255)
                                                                                   tCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTAGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTgg                                                                                                  <  1:443966/117‑1 (MQ=255)
                                                                                     ttGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtg                                                                                                            <  2:62910/105‑1 (MQ=255)
                                                                                     ttGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtg                                                                                                            >  1:62910/1‑105 (MQ=255)
                                                                                     ttGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGggt                                                                      >  2:924109/1‑143 (MQ=255)
                                                                                     ttGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGggt                                                                      <  1:924109/143‑1 (MQ=255)
                                                                                        tttATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCggg                                                             >  4:107888/1‑149 (MQ=255)
                                                                                            tCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCGTTg                                                                                         <  1:98595/117‑1 (MQ=255)
                                                                                            tCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTg                                                                                         >  2:98595/1‑117 (MQ=255)
                                                                                              ccgccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGggt                                                                      <  2:410553/134‑1 (MQ=255)
                                                                                              ccgccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGAATTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGggt                                                                      >  1:410553/1‑134 (MQ=255)
                                                                                               cgccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCtt                                                                                     >  3:855254/1‑118 (MQ=255)
                                                                                               cgccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGAATTGGCtt                                                                                     <  4:855254/118‑1 (MQ=255)
                                                                                                gccgTTGCTGTTCGCTGATGGCTGGTAAACGCCGACGCGTGAATTTCTCGAACATGGTCTAGATATTTTGGGCCTCGCGCAGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCAGTATTTACAGGGTGATGCCGGGCATTCCGc                                                     <  4:961086/149‑1 (MQ=255)
                                                                                                 ccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCt                                                    <  2:488857/149‑1 (MQ=255)
                                                                                                 ccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCt                                                    <  4:957071/149‑1 (MQ=255)
                                                                                                 ccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCt                                                    <  3:471259/149‑1 (MQ=255)
                                                                                                    ttGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGAt                                                 <  3:817963/149‑1 (MQ=255)
                                                                                                    ttGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGAt                                                 <  3:300520/149‑1 (MQ=255)
                                                                                                      gCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATcc                                               >  1:260770/1‑149 (MQ=255)
                                                                                                      gCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATcc                                               <  1:10372/149‑1 (MQ=255)
                                                                                                        tGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCt                                             >  1:166732/1‑149 (MQ=255)
                                                                                                        tGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCt                                             <  4:604339/149‑1 (MQ=255)
                                                                                                          ttCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                           <  2:260770/149‑1 (MQ=255)
                                                                                                          ttCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                           >  2:846147/1‑149 (MQ=255)
                                                                                                                gATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTgc                                     <  3:379218/149‑1 (MQ=255)
                                                                                                                      tGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggc                               >  2:933955/1‑149 (MQ=255)
                                                                                                                       ggAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggcg                              >  2:731939/1‑149 (MQ=255)
                                                                                                                           aaCGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtgg                                                                                                           <  2:146586/68‑1 (MQ=255)
                                                                                                                           aaCGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtgg                                                                                                           >  1:146586/1‑68 (MQ=255)
                                                                                                                             cGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTg                        <  1:177321/149‑1 (MQ=255)
                                                                                                                              gCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGc                       <  1:846147/149‑1 (MQ=255)
                                                                                                                              gCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGc                       <  2:685223/149‑1 (MQ=255)
                                                                                                                              gCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGc                       >  3:253697/1‑149 (MQ=255)
                                                                                                                               ccGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGGGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCGCTGCCTTTGCGCTGGCGGCGGTGCt                      <  1:731939/149‑1 (MQ=255)
                                                                                                                                 gACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCttt                    <  1:963175/149‑1 (MQ=255)
                                                                                                                                 gACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCttt                    <  2:906312/149‑1 (MQ=255)
                                                                                                                                  aCCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTtc                   <  4:115046/149‑1 (MQ=255)
                                                                                                                                   cccGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTtct                  <  4:666829/149‑1 (MQ=255)
                                                                                                                                    ccGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTtctc                 >  3:670075/1‑149 (MQ=255)
                                                                                                                                      gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGc                                                     <  4:289685/111‑1 (MQ=255)
                                                                                                                                      gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGc                                                     >  3:289685/1‑111 (MQ=255)
                                                                                                                                      gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                           >  3:859500/1‑121 (MQ=255)
                                                                                                                                      gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                           <  2:750220/121‑1 (MQ=255)
                                                                                                                                      gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                           >  1:750220/1‑121 (MQ=255)
                                                                                                                                      gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                           <  4:859500/121‑1 (MQ=255)
                                                                                                                                            ttCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGAt         <  3:450979/149‑1 (MQ=255)
                                                                                                                                             tctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATg        >  3:2178/1‑149 (MQ=255)
                                                                                                                                              ctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATcc                                               <  2:165912/109‑1 (MQ=255)
                                                                                                                                              ctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATcc                                               >  1:165912/1‑109 (MQ=255)
                                                                                                                                                 gAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCtgt    >  2:599362/1‑149 (MQ=255)
                                                                                                                                                   aCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTgg  <  4:1038499/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGG  >  NC_000913/4280072‑4280367

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