Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,561,205 A→T 5.7% intergenic (+307/‑163) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,561,2050AT 5.7% 277.5 / 4.7 100intergenic (+307/‑163)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (47/47);  new base T (4/2);  total (51/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

TTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAA  >  NC_000913/2561059‑2561352
                                                                                                                                                  |                                                                                                                                                   
ttAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAg                                                                                                                                                   >  2:1232132/1‑149 (MQ=255)
  aaGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                 >  2:626234/1‑149 (MQ=255)
   aGAATTCTAGGCGTCGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGGCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTg                                                                                                                                                <  2:1059659/149‑1 (MQ=255)
    gAATTCTAGGCGATGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTgg                                                                                                                                               >  1:1443470/1‑149 (MQ=255)
       ttCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGaaa                                                                                                                                            <  2:263507/149‑1 (MQ=255)
        tatAGGCGACGGGAGTTTGATAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGACGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAAc                                                                                                                                           <  2:1494644/147‑1 (MQ=255)
          tAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGc                                                                                                                                         <  2:686826/149‑1 (MQ=255)
          tAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGc                                                                                                                                         <  1:1389747/149‑1 (MQ=255)
           aGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGcc                                                                                                                                        >  1:834155/1‑149 (MQ=255)
            ggCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCGGGaaatgatt                                                                                                                                       <  1:818169/149‑6 (MQ=255)
               gACGGGAGGTTTTTAGCCTCCCGGTCGGGCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTAGTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCAtt                                                                                                                                    <  2:878422/149‑1 (MQ=255)
                 cGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTTACCATGCTGATTAGCGCAATCATTTCCAGCAAGCGGGAACCGCCCATtcc                                                                                                                                  >  1:1043896/1‑148 (MQ=255)
                  gggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                 <  1:487641/149‑1 (MQ=255)
                    gAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactc                                                                                                                               >  1:1400320/1‑149 (MQ=255)
                    gAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAAACATTTCCAGCAAGCTGGAAACGCCCAAtcactc                                                                                                                               >  2:705987/1‑149 (MQ=255)
                    caGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCAAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAACCGCCCATtcactc                                                                                                                               >  2:838556/2‑149 (MQ=255)
                     aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactca                                                                                                                              >  2:990364/1‑149 (MQ=255)
                     aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGGTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactca                                                                                                                              >  1:868587/1‑149 (MQ=255)
                      gTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCACAATCATTTCCAGCAAGCTGGAAACGCCCATtcactcac                                                                                                                             <  1:142719/149‑1 (MQ=255)
                       ttttttAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACg                                                                                                                            <  2:444282/149‑1 (MQ=255)
                          tttAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCGATTCACTCACGcac                                                                                                                         >  2:981957/1‑149 (MQ=255)
                          tttAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGcac                                                                                                                         <  1:1438710/149‑1 (MQ=255)
                           ttAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGcaca                                                                                                                        >  1:1557682/1‑149 (MQ=255)
                             aGCCTCCCGGTCGGTGTATGAGAGTCGCATATAGTTAAACAGAAGAATTATTTCGCCGACTGATAGGTCGGCTGGTTGAGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAc                                                                                                                      <  1:1310213/149‑1 (MQ=255)
                              gCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAcc                                                                                                                     <  1:30795/149‑1 (MQ=255)
                               ccTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCt                                                                                                                    <  2:1531750/149‑1 (MQ=255)
                                  cccGGTCGGTCATAGAGAGGCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGaa                                                                                                                 <  2:637434/149‑1 (MQ=255)
                                      gtcggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagc                                                                                                             <  2:1277122/149‑1 (MQ=255)
                                         ggtgATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGa                                                                                                          <  1:190523/149‑1 (MQ=255)
                                         ggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTTACCATGCTGATTAGAGCAATCATTTCCAGCAAGGTGGAAACGCCCATTCACTCACGCACCACCGGAActgctgc                                                                                                          >  2:497496/1‑146 (MQ=255)
                                         ggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGa                                                                                                          <  2:437771/149‑1 (MQ=255)
                                          gtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGAc                                                                                                         <  2:1163756/149‑1 (MQ=255)
                                             aTAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGCCcaa                                                                                                      >  2:586984/1‑149 (MQ=255)
                                             aTAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCAGTCACTCACGCACCACCTGAACAGCAGAGcaa                                                                                                      >  2:477770/1‑149 (MQ=255)
                                               agagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTATCGCAATCATTTCAAGCAAGCTGTAAACGCCCATTCACTCAAGCACAACATGAACCGCAGACcaaca                                                                                                    >  2:1186840/1‑149 (MQ=255)
                                                 atagTCGCCTAGAGTTAAACAGAAGAATTATTTCGACGCCAGATAGGTAGGCTGGTTGCGAGCATTAACCATGCTGATTAGCTCAATCAGTTCCAGCTAGCTGTACACGCCCAGTCACAAACGCAAAACCTGCACAGCAGACcaacacc                                                                                                  >  2:1059651/3‑147 (MQ=255)
                                                     gcGCCTAGAGTTAAACAGATGAATTATTTGGCCGCCTGCTAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACCAACAACCTGAACAGCTGACCATCTATCcat                                                                                              <  2:1068427/148‑1 (MQ=255)
                                                        ccTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGCCGGCTGGTTGCGAGCATTAACCATGCTGATTATCGAAATCATTTACAGCAATCTGGAAACGCCCATTCACTCACGCACACCCTGAACAGCAGACCAACAAACcatcat                                                                                           >  1:1100279/1‑149 (MQ=255)
                                                         cTAGAGTTAATCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTCGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACACCCTGAACAGCAGACCCACAAACCCTCATg                                                                                          >  1:1203700/1‑149 (MQ=255)
                                                         cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATg                                                                                          >  2:546912/1‑149 (MQ=255)
                                                          tAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTGGGCTGGTTGCGAGCATTAACCATGCGGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATgg                                                                                         <  2:1140043/149‑1 (MQ=255)
                                                               ttAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTc                                                                                    >  1:200345/1‑149 (MQ=255)
                                                                  aaCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGtt                                                                                 <  2:1443470/149‑1 (MQ=255)
                                                                     agaagaATTATTTGTCCGCCTGATAGGTCGGCTGGTTGCGAGAATTAACCATGCTGATTATCGCAATCATTTCCATCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                              <  1:1062567/149‑1 (MQ=255)
                                                                     agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTTGCGCAATCATTTCCAGCAAGCTGGAACCGCCCATTCACTCACTCACCACCTGCACAGCAGACCACCACACCCTCATGGCATTCGATCAt                                                                              >  2:1019777/1‑149 (MQ=255)
                                                                     agaagaATTATTTCGCAGCCTGAGAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCTCAATCATATCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                              <  1:290066/149‑1 (MQ=255)
                                                                      gatgaaGTATTTCGCCGCCTGATAGGTCGGCTGGATGCGAGCATTAACCATGCTGATTAGCGCACTCATTTCCAGCAAGATGGAAACGCGCATGCCCTCCCGCACACCCTGCACAGCAGCCCAACAAACCAACAGGGCAATCGTGCATg                                                                             >  2:1469150/4‑149 (MQ=255)
                                                                      gaagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAACAGACCAACAAACCATCATGGCATTCGTTCATg                                                                             <  2:308950/149‑1 (MQ=255)
                                                                        agaaATATTTCGCCGCCTGATAGGGCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATggg                                                                           <  2:1353527/149‑1 (MQ=255)
                                                                          aaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGTGCAGCAGACCAACAACCCATCATGGCATTCGTTCTTggggt                                                                         >  2:991171/1‑148 (MQ=255)
                                                                             tattTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgcg                                                                      <  2:369035/149‑1 (MQ=255)
                                                                               tttCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGcgcg                                                                    <  1:8665/149‑1 (MQ=255)
                                                                               tttCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCAGGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAACCGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCAGTCGTTCATGGGGAGctcg                                                                    >  1:1356143/1‑149 (MQ=255)
                                                                                   gccgccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGtat                                                                <  1:626234/149‑1 (MQ=255)
                                                                                     cgccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTatat                                                              <  2:451713/149‑1 (MQ=255)
                                                                                       ccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATAtct                                                            >  2:214224/1‑149 (MQ=255)
                                                                                       ccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATAtct                                                            <  1:1428399/149‑1 (MQ=255)
                                                                                           aTAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACCAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCaaa                                                        >  1:1297953/1‑149 (MQ=255)
                                                                                               gTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCTTGGGGAGCGCGGTATATCTCAAAGGGt                                                    >  1:926226/1‑149 (MQ=255)
                                                                                               gTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGt                                                    <  2:1202876/149‑1 (MQ=255)
                                                                                                tCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAACCGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTg                                                   >  2:901751/1‑149 (MQ=255)
                                                                                                tCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACACACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTg                                                   >  1:458902/1‑149 (MQ=255)
                                                                                                  ggCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCTCGCACAACCTGACCAGCAGACCAACTAACCATCATGGCATTCGTTCAGGGGGAGCGGGGGATTTGTCTAAGtgtggt                                                 >  2:1015179/1‑147 (MQ=255)
                                                                                                  ggCTGGTTGAGAGCATTAACCATGCTGATTATCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACTCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGTAGCGCGGTATATCTCAAAGGGTGAt                                                 <  1:362908/149‑1 (MQ=255)
                                                                                                           cGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                        <  1:1342362/149‑1 (MQ=255)
                                                                                                           cGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                        <  1:838556/149‑1 (MQ=255)
                                                                                                            gAGCATTAACCATGCTGCTTAGCGCAATCCTTTCCCGCAAGCTGGAAACGCCCATTCACTACCGCACAACCTGAACAGCAGACCAACCAACCCTCAGGGCATTCGTTGATGGGGAGCGCGGTATCTCTCAAAGGGTGATGTCAGATtgg                                       >  2:1137382/1‑149 (MQ=255)
                                                                                                            gAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                       >  1:127366/1‑149 (MQ=255)
                                                                                                              gCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtg                                     <  2:1557682/149‑1 (MQ=255)
                                                                                                               cATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtgg                                    <  2:834155/149‑1 (MQ=255)
                                                                                                                aTTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAACCGCCCATTCACTCACGCACAACCTGAACAGCAGCCCAACAAAGCATCAGGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCATTATGGTggg                                   >  2:1154579/1‑149 (MQ=255)
                                                                                                                aTTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                       <  1:892536/145‑1 (MQ=255)
                                                                                                                aTTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                       >  2:892536/1‑145 (MQ=255)
                                                                                                                 ttAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGc                                  <  1:598436/149‑1 (MQ=255)
                                                                                                                     ccattgtGATTAGAGCAAGCATTTCCAGCAAGCTGGATAGGCCCATTCACTCACGCACAACCTGAACAGCAGACAATCAATACTTCATTGCACTCGTTCATGGTGAGCGCTGTATATCTCAAACGGTGATTTCTTATTGGTTGGCATTc                              <  2:720733/143‑1 (MQ=255)
                                                                                                                     ccATGCTGATTAGCGCAATCATTTCCAGCAAGCTCGAAACGCCCATTACCTCACGCACAACCTGAACTGCAGACCAACAAACCATCATGGCCTTCGTTCATCGGGAGCGCGCTATAACTCAAAGGGTGATATCTTATAGGTGGGCCTTc                              >  1:353947/1‑149 (MQ=255)
                                                                                                                        tGCTCATTAGCGCAATCATATCCAGCAAGCTGGAAACGCCCAATCACTCACGCACAAGCTGAACAGCAGACCAACAAACCATCATGGAATTCGTACAAGGGGAGCCGGGTTTATCTCAAAGCGTCATTGCTTATTGGTGGGCCTTCGGa                           >  2:1302309/1‑149 (MQ=255)
                                                                                                                         gCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGaa                          <  1:457586/149‑1 (MQ=255)
                                                                                                                         gCTGAATAGCGCAATCATTTCAAGCAAGCTGGAAACGCCCATTCACTCACGCACAACATGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGaa                          <  1:1186840/149‑1 (MQ=255)
                                                                                                                            gATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGc                       <  1:990364/149‑1 (MQ=255)
                                                                                                                              ttAGCGCAATCATTTCCAGCAAGCTGGAAACGCCGATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCtt                     <  1:981957/149‑1 (MQ=255)
                                                                                                                                  cgcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGATCAGCAGACCAACAAACCATCATCGCATGCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAg                 <  1:256518/149‑1 (MQ=255)
                                                                                                                                    cAATCATTTCCAGCAAGCTTGAACCGCCCATTCACTCCCGCTCAACCTGAACAGCCGACCAACAAACGCGCATGGCATTAGTTCATGGGGAGCGCGGTATATCTCAAAGTGTGATATCTTAATGGTGGGCCTTCGGAAGGCTTCGAGCa               >  2:1454635/1‑149 (MQ=255)
                                                                                                                                    cAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACACCCTGAACAGCAGACCAACACACCATCATGGGATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGATGGCTTCGAGCa               >  1:629208/1‑149 (MQ=255)
                                                                                                                                     gggcgtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCaa              >  2:216311/6‑149 (MQ=255)
                                                                                                                                     aaTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCaa              >  1:255461/1‑149 (MQ=255)
                                                                                                                                       tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAAtt            >  2:1280700/1‑149 (MQ=255)
                                                                                                                                       tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAAtt            <  2:1297953/149‑1 (MQ=255)
                                                                                                                                       tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTAATGGGGCGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGGGGGCCTGCGGAAGGCTTCGAGCAAtt            >  2:871862/1‑149 (MQ=255)
                                                                                                                                       gcgtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                        <  2:1208695/118‑1 (MQ=255)
                                                                                                                                       gcgtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                        >  1:1208695/4‑121 (MQ=255)
                                                                                                                                        cgtttCCAGCTTGCTGGAAATGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTGTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTg           >  2:771987/3‑149 (MQ=255)
                                                                                                                                         aTTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTGATAGGTGGGCCTTCGGAAGGCTTCGAGCAATTGc          >  1:719747/1‑149 (MQ=255)
                                                                                                                                          tttCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCAGGATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTGGGAAGGCTTCGAGCAATTGcc         >  1:252921/1‑149 (MQ=255)
                                                                                                                                          tttCCAGCAAGCTGGAAACGCCCATTCACTCACGCAAAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGca         >  2:1318964/1‑148 (MQ=255)
                                                                                                                                            tCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTACATGGGGAGCGCGGTATATCTCAAAGGGTCATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCac       >  1:807694/1‑148 (MQ=255)
                                                                                                                                             ccAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGc      <  1:29175/149‑1 (MQ=255)
                                                                                                                                               cgCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCtt    <  2:775417/148‑1 (MQ=255)
                                                                                                                                               aGCAAGCTGGAACCGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGACCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCtt    >  1:480661/1‑149 (MQ=255)
                                                                                                                                                 cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAGGGGGAGCGCGGTATATCTCAATGGGTGATTTCGTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTaa  >  1:1379531/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAA  >  NC_000913/2561059‑2561352

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