Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,164,284 A→T 5.4% K200* (AAG→TAG) ‡ PP_3665 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,164,2840AT 5.4% 342.1 / 8.6 111K200* (AAG→TAG) ‡PP_3665AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (51/54);  new base T (2/4);  total (53/58)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

Only 100 of 111 total aligned reads displayed.

TTGCAGGCGCACTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTGGAGGAGCCGCTGGGGCT  >  NC_002947/4164137‑4164432
                                                                                                                                                   |                                                                                                                                                    
ttGCAGGCGCACTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAg                                                                                                                                                    >  1:729709/1‑150 (MQ=255)
        gcACTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGCACTGGACCGCNACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTTGCCCGCGACACGGGCGGCCAGGACATTGCGGCCAAGGCCGCCGa                                                                                                                                            >  1:157970/1‑150 (MQ=255)
        gcACTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCTTGgccgcgat                                                                                                                                            <  1:499590/150‑4 (MQ=255)
        gcACTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCTTGgccgcgat                                                                                                                                            <  2:1685019/150‑4 (MQ=255)
        gcACTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGa                                                                                                                                            <  1:197081/150‑1 (MQ=255)
        gcACTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGa                                                                                                                                            >  2:1632842/1‑150 (MQ=255)
          aCTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTAGACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCTTGgccgcgat                                                                                                                                            >  2:1022238/1‑145 (MQ=255)
          aCTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCTTGgccgcgat                                                                                                                                            <  1:1022238/148‑4 (MQ=255)
           cTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCt                                                                                                                                         <  2:384538/150‑1 (MQ=255)
            ttCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCTATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTg                                                                                                                                        >  1:702751/1‑150 (MQ=255)
            ttCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTg                                                                                                                                        >  2:746290/1‑150 (MQ=255)
             tCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGt                                                                                                                                       <  1:1855978/150‑1 (MQ=255)
             tCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATAGCGGCCAAGGCCGCCGAGCTGt                                                                                                                                       <  1:67143/150‑1 (MQ=255)
               cccGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTg                                                                                                                                     >  2:805693/1‑150 (MQ=255)
                   gTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACTGAAACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCtgt                                                                                                                                 <  2:448298/150‑1 (MQ=255)
                   gTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCTATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCtgt                                                                                                                                 <  2:702751/150‑1 (MQ=255)
                     aTCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGc                                                                                                                               <  2:1432954/150‑1 (MQ=255)
                     aTCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGc                                                                                                                               <  1:366028/150‑1 (MQ=255)
                      aggTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCg                                                                                                                              >  1:1221495/3‑150 (MQ=255)
                         tGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGAcc                                                                                                                           >  1:1264430/1‑150 (MQ=255)
                         tGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATCGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGAcc                                                                                                                           <  2:1038993/150‑1 (MQ=255)
                          gATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACcg                                                                                                                          >  2:674482/1‑150 (MQ=255)
                              tcCAGCCAGCTGTTCGCACTGGACCGCGACCGTTACACCTGCTCNGGCGNCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACGAGGGCGCCAAGGCCGCCGAGCTGTTGCTGTGGGACCGCGTa                                                                                                                      >  1:18151/1‑150 (MQ=255)
                                  gccagcTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTAcgcg                                                                                                                  <  2:815284/150‑1 (MQ=255)
                                    cagcTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGgt                                                                                                                <  1:657864/150‑1 (MQ=255)
                                     agcTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGgtg                                                                                                               >  1:827853/1‑150 (MQ=255)
                                          ttCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGgtgagcgt                                                                                                          <  1:1632842/150‑1 (MQ=255)
                                            cGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGgtgagcgtga                                                                                                        <  1:674482/150‑1 (MQ=255)
                                            cGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGgtgagcgtga                                                                                                        <  2:1550227/150‑1 (MQ=255)
                                            cGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGgtgagcgtga                                                                                                        <  1:1612977/150‑1 (MQ=255)
                                                cTGGACCGCGACCGTTACACCTGTTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAgcgc                                                                                                    <  2:729709/150‑1 (MQ=255)
                                                 ggtaCCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGcc                                                                                                   >  1:1221463/4‑150 (MQ=255)
                                                  ggACCGCGACCTTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTg                                                                                                                                     <  2:1503664/115‑1 (MQ=255)
                                                  ggACCGCGACCTTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTg                                                                                                                                     >  1:1503664/1‑115 (MQ=255)
                                                  ggACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCa                                                                                                  <  2:1109011/150‑1 (MQ=255)
                                                  ggACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCa                                                                                                  >  2:294459/1‑150 (MQ=255)
                                                   gACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAg                                                                                                 >  1:1703179/1‑150 (MQ=255)
                                                     ccGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTAAGCGGTGAGCTTGAGCTCCAGCg                                                                                               >  2:585974/1‑150 (MQ=255)
                                                      cgagACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCgt                                                                                              >  2:1818194/4‑150 (MQ=255)
                                                            cGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCgtgt                                                                                            >  1:1611740/1‑146 (MQ=255)
                                                              ttACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCTTGgccgc                                                                                                                                               <  1:1227008/93‑1 (MQ=255)
                                                              ttACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCTTGgccgc                                                                                                                                               >  2:1227008/1‑93 (MQ=255)
                                                              taacacCTGCGCAGGACGCGTCGCGGACACGGGCCTGATTGTGCAACTGTTGGCGCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCt                                                                                      <  2:18151/148‑1 (MQ=255)
                                                                acacCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTgc                                                                                    <  2:1071714/150‑1 (MQ=255)
                                                                acacCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGCTCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTgc                                                                                    <  1:881777/150‑1 (MQ=255)
                                                                  acCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGCGCGCCGGCGTGTGCCGCTgcgc                                                                                  >  2:125465/1‑150 (MQ=255)
                                                                          aGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGcc                                                                          >  1:784936/1‑149 (MQ=255)
                                                                               gcgTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCa                                                                     <  2:314074/150‑1 (MQ=255)
                                                                                 gTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAgc                                                                   >  1:1460704/1‑150 (MQ=255)
                                                                                   ggcggcCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAgcgc                                                                 >  2:1078742/1‑150 (MQ=255)
                                                                                    gcggcCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATTGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCt                                                                <  2:1697196/150‑1 (MQ=255)
                                                                                      ggcCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCa                                                                                 <  1:547710/131‑1 (MQ=255)
                                                                                      ggcCATGGACATGATGTTGCAACCGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCa                                                              <  1:1120942/150‑1 (MQ=255)
                                                                                         tttGGACATGAGGTTGCAACCGGTGGCCCGCGACAGGGGCGGCCAGGACATCGCGGCCACGGCCGCTGCGCTGTGGCTGTCCGACCGCGTGCGCGGTCAGCGTGACCGCCACACTGTGCCGCGGCGCACCGTGCTGGGCAGCCCTCAGcc                                                           <  2:827853/148‑1 (MQ=255)
                                                                                         cATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCgtgt                                                                                            >  2:1599168/1‑117 (MQ=255)
                                                                                         cATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCgtgt                                                                                            <  1:1599168/117‑1 (MQ=255)
                                                                                         cATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCATcc                                                           <  1:648574/150‑1 (MQ=255)
                                                                                          aTGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGccc                                                          >  2:1923239/1‑150 (MQ=255)
                                                                                          aTGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGccc                                                          >  1:615216/1‑150 (MQ=255)
                                                                                             gACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCaaa                                                       <  2:512184/150‑1 (MQ=255)
                                                                                              aCATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAAc                                                      <  1:1132218/150‑1 (MQ=255)
                                                                                               cATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGTGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACt                                                     >  1:1645085/1‑150 (MQ=255)
                                                                                                atgatgTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCa                                                                     >  1:1337344/1‑133 (MQ=255)
                                                                                                atgatgTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCa                                                                     <  2:1337344/133‑1 (MQ=255)
                                                                                                 tgatgTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGa                                                   <  2:1221463/150‑1 (MQ=255)
                                                                                                 tgatgTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGa                                                   <  2:1629440/150‑1 (MQ=255)
                                                                                                    tgTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGc                                                            <  1:1532212/138‑1 (MQ=255)
                                                                                                    tgTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGc                                                            >  2:1532212/1‑138 (MQ=255)
                                                                                                      ttGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCgt                                                                                              >  2:1472513/1‑102 (MQ=255)
                                                                                                      ttGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCgt                                                                                              <  1:1472513/102‑1 (MQ=255)
                                                                                                            cTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGgtgagcgtgag                                                                                                       <  1:1085826/87‑1 (MQ=255)
                                                                                                            cTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGgtgagcgtgag                                                                                                       >  2:1085826/1‑87 (MQ=255)
                                                                                                            cTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGgtggtg                                        >  1:1328264/1‑150 (MQ=255)
                                                                                                             tggtggCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGtggtgg                                       >  2:1071988/1‑150 (MQ=255)
                                                                                                                tggCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCCCCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGa                                    <  1:805693/150‑1 (MQ=255)
                                                                                                                  gCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATc                                  >  1:1880787/1‑150 (MQ=255)
                                                                                                                    ccGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCAt                                >  2:1485108/1‑150 (MQ=255)
                                                                                                                    ccGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCAt                                >  1:1175447/1‑150 (MQ=255)
                                                                                                                    ccGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGCGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCAt                                >  2:1683560/1‑150 (MQ=255)
                                                                                                                     cgcgACAGAGGCTGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATg                               <  1:1021035/150‑1 (MQ=255)
                                                                                                                       cgACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGa                             <  2:1175447/150‑1 (MQ=255)
                                                                                                                        gACACGGGCCGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGTCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGTTCATGGAg                            <  1:585974/150‑1 (MQ=255)
                                                                                                                         acacGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAgg                           >  2:200932/1‑150 (MQ=255)
                                                                                                                         acacGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAgg                           <  2:1460704/150‑1 (MQ=255)
                                                                                                                         acacGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAgg                           <  2:52170/150‑1 (MQ=255)
                                                                                                                         acacGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAgg                           >  1:1424844/1‑150 (MQ=255)
                                                                                                                             gggCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCaa                       >  2:1583036/1‑150 (MQ=255)
                                                                                                                             gggCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGAGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCaa                       >  1:417055/1‑150 (MQ=255)
                                                                                                                                   ccAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGc                                                            >  2:1805234/1‑107 (MQ=255)
                                                                                                                                   ccAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGc                                                            <  1:1805234/107‑1 (MQ=255)
                                                                                                                                    cAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTggag                <  1:1485108/150‑1 (MQ=255)
                                                                                                                                      ggACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTggagga              <  1:860062/150‑1 (MQ=255)
                                                                                                                                       gACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGCGCCAGGTGGTGGCGATCATGGATGCCAATCTGgaggag             <  2:1221495/150‑1 (MQ=255)
                                                                                                                                       gACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTGgaggag             >  1:1917332/1‑150 (MQ=255)
                                                                                                                                        aCATCGCGGCCAAGGCCGCCGAGCCGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTGGAGGAGc            <  2:1707350/150‑1 (MQ=255)
                                                                                                                                           tCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTGGAGGAGCCGc         >  1:1493384/1‑150 (MQ=255)
                                                                                                                                           tCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTGGAGGAGCCGc         >  2:846095/1‑150 (MQ=255)
                                                                                                                                                 ccAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTGGAGGAGCCGCTgg      >  2:259999/1‑147 (MQ=255)
                                                                                                                                                 ccAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTGGAGGAGCCGCTgg      <  1:259999/147‑1 (MQ=255)
                                                                                                                                                  cAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTGGAGGAGCCGCTGGGGCt  <  1:720140/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGCAGGCGCACTTCCCCGGTATCGTGATCTCCAGCCAGCTGTTCGAACTGGACCGCGACCGTTACACCTGCTCAGGCGGCGTGGCGGCCATGGACATGATGTTGCAACTGGTGGCCCGCGACACGGGCGGCCAGGACATCGCGGCCAAGGCCGCCGAGCTGTTGCTGTGCGACCGCGTACGCGGTGAGCGTGAGCGCCAGCGTGTGCCGCTGCGCACCCTGCTCGGCAGCGCTCAGCCCAAACTGAGCCAGGTGGTGGCGATCATGGAGGCCAATCTGGAGGAGCCGCTGGGGCT  >  NC_002947/4164137‑4164432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: