Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,280,224 G→T 6.1% G82V (GGT→GTT)  yjcE → putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,280,2240GT 6.1% 346.3 / 10.2 98G82V (GGT→GTT) yjcEputative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (41/51);  new base T (3/3);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

TGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGGC  >  NC_000913/4280077‑4280368
                                                                                                                                                   |                                                                                                                                                
tGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGt                                                                                                                                                 >  1:268796/1‑149 (MQ=255)
tGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGt                                                                                                                                                 >  3:270088/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                                                                                                       >  4:493759/1‑149 (MQ=255)
                                          gtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCAcc                                                                                                       >  3:150659/1‑149 (MQ=255)
                                          gtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCAcc                                                                                                       >  2:684191/1‑149 (MQ=255)
                                           tttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCg                                                                                                      <  4:371954/149‑1 (MQ=255)
                                           tttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCg                                                                                                      >  4:64953/1‑149 (MQ=255)
                                           tttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCg                                                                                                      >  4:471259/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:739129/149‑1 (MQ=255)
                                                         ttGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGc                                                                                        <  1:1128418/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                                                                                 >  2:177321/1‑149 (MQ=255)
                                                                tGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTa                                                                                 <  1:222359/149‑1 (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                                                                                                             >  1:62910/1‑105 (MQ=255)
                                                                                ttGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtg                                                                                                             <  2:62910/105‑1 (MQ=255)
                                                                                ttGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGggt                                                                       <  1:924109/143‑1 (MQ=255)
                                                                                ttGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGggt                                                                       >  2:924109/1‑143 (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                                                     <  4:957071/149‑1 (MQ=255)
                                                                                            ccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCt                                                     <  2:488857/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                                              <  4:604339/149‑1 (MQ=255)
                                                                                                   tGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCt                                              >  1:166732/1‑149 (MQ=255)
                                                                                                     ttCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                            >  2:846147/1‑149 (MQ=255)
                                                                                                     ttCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                            <  2:260770/149‑1 (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                        >  3:253697/1‑149 (MQ=255)
                                                                                                                         gCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGc                        <  2:685223/149‑1 (MQ=255)
                                                                                                                          ccGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGGGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCGCTGCCTTTGCGCTGGCGGCGGTGCt                       <  1:731939/149‑1 (MQ=255)
                                                                                                                            gACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCttt                     <  2:906312/149‑1 (MQ=255)
                                                                                                                            gACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCttt                     <  1:963175/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                                            <  4:859500/121‑1 (MQ=255)
                                                                                                                                 gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                            >  3:859500/1‑121 (MQ=255)
                                                                                                                                 gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                            <  2:750220/121‑1 (MQ=255)
                                                                                                                                 gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                            >  1:750220/1‑121 (MQ=255)
                                                                                                                                       ttCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGAt          <  3:450979/149‑1 (MQ=255)
                                                                                                                                        tctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATg         >  3:2178/1‑149 (MQ=255)
                                                                                                                                         ctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATcc                                                >  1:165912/1‑109 (MQ=255)
                                                                                                                                         ctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATcc                                                <  2:165912/109‑1 (MQ=255)
                                                                                                                                            gAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCtgt     >  2:599362/1‑149 (MQ=255)
                                                                                                                                              aCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTgg   <  4:1038499/149‑1 (MQ=255)
                                                                                                                                               cATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGgc  >  3:648805/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGGC  >  NC_000913/4280077‑4280368

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: