Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,530,671 A→C 10.6% intergenic (+19/+44) hslO → / ← yhgE molecular chaperone Hsp33/putative transport protein YhgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,530,6710AC10.6% 287.3 / 26.4 113intergenic (+19/+44)hslO/yhgEmolecular chaperone Hsp33/putative transport protein YhgE
Reads supporting (aligned to +/- strand):  ref base A (57/44);  new base C (6/6);  total (63/50)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

Only 100 of 114 total aligned reads displayed.

TAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTA  >  NC_000913/3530523‑3530815
                                                                                                                                                    |                                                                                                                                                
tAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCagagag                                                                                                                                                 <  2:623289/150‑1 (MQ=255)
 aGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAgg                                                                                                                                                <  1:342819/150‑1 (MQ=255)
 aGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAgg                                                                                                                                                >  1:81147/1‑150 (MQ=255)
 aGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAgg                                                                                                                                                <  2:309604/150‑1 (MQ=255)
   cATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGc                                                                                                                                              >  2:881539/1‑150 (MQ=255)
     tCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCTGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGctc                                                                                                                                            <  2:81147/150‑1 (MQ=255)
      ccTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGctct                                                                                                                                           <  1:606029/150‑1 (MQ=255)
       cTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTg                                                                                                                                          <  1:1600533/150‑1 (MQ=255)
        tGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGc                                                                                                                                         >  2:229105/1‑150 (MQ=255)
         ggCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGcc                                                                                                                                        >  2:332094/1‑150 (MQ=255)
          gCGGAATATGTCTAAATTGACATGCATTGTGATTACTGCGGTAAGCTCTATATGTTCAATGCGATGGATATTGCTGTTATCCGCAACTACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCg                                                                                                                                       <  1:670665/150‑1 (MQ=255)
          gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCg                                                                                                                                       <  2:224117/150‑1 (MQ=255)
          gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCg                                                                                                                                       <  2:1389866/150‑1 (MQ=255)
          gCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGACATGGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCg                                                                                                                                       <  1:806184/150‑1 (MQ=255)
           cGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCgg                                                                                                                                      <  1:1541135/150‑1 (MQ=255)
              aaGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGCCCTCTCTGCCGGaca                                                                                                                                   <  1:370433/150‑3 (MQ=255)
                  tGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTAc                                                                                                                               <  1:1156006/150‑1 (MQ=255)
                  tGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTAc                                                                                                                               <  1:365840/150‑1 (MQ=255)
                      caaaTTGACTCGCATTGTGATTACGCCGGAACCCCCAAGCTGTTCAATCCCATGGATATTGCTCAAACCCCCAACAAGCCGGGTCCGGATTACCGGCACGTCCATTATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTa                                                                                                                          <  2:1538137/149‑1 (MQ=255)
                        aaTTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTaa                                                                                                                         <  2:437061/150‑1 (MQ=255)
                          ttGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAAGAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATa                                                                                                                       <  1:1217943/150‑1 (MQ=255)
                            gACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATaca                                                                                                                     >  2:1601031/1‑150 (MQ=255)
                                    ttGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACAcc                                                                                                                   <  2:416868/144‑1 (MQ=255)
                                    ttGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACAcc                                                                                                                   >  1:416868/1‑144 (MQ=255)
                                      gtgATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATgcgc                                                                                                           >  2:781577/1‑150 (MQ=255)
                                       tgATTACTGCGGTAACCACTATCTGTTCAATGCGCTGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCt                                                                                                          <  1:720524/150‑1 (MQ=255)
                                       tgATTACTGCGGTAAACACTATCTGTTCAATGCGAAGGATATTGCTCAAATACGCAACAACGCGTATCCGGCAGATCCGCAAGTTCATTAATGACTAAGTCCGGCAGAGAGGGCTCTGCCTGATTTTACGACTAATACACCTGATgcgca                                                                                                          >  2:37593/1‑149 (MQ=255)
                                        gATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCtt                                                                                                         >  2:884547/1‑150 (MQ=255)
                                        gATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCtt                                                                                                         <  2:930653/150‑1 (MQ=255)
                                            aCTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTc                                                                                                     <  2:33101/150‑1 (MQ=255)
                                             cTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTcc                                                                                                    >  2:738973/1‑149 (MQ=255)
                                             cTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCa                                                                                                    >  2:1389517/1‑150 (MQ=255)
                                             cTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTCTACGACTAATACACCTGATGCGCTTTGTcc                                                                                                    >  1:273327/1‑149 (MQ=255)
                                              tGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCAc                                                                                                   <  2:730485/150‑1 (MQ=255)
                                                 ggTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCa                                                                                                >  2:174001/1‑150 (MQ=255)
                                                     aCCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCa                                                                                            <  2:1094258/150‑1 (MQ=255)
                                                     aCCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCa                                                                                            <  1:942979/150‑1 (MQ=255)
                                                        aCTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTTTGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGAt                                                                                         <  1:998404/150‑1 (MQ=255)
                                                           aTCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTca                                                                                      >  1:1036501/1‑150 (MQ=255)
                                                           aTCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTca                                                                                      >  1:70624/1‑150 (MQ=255)
                                                            tCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTcac                                                                                     >  1:863786/1‑150 (MQ=255)
                                                             cTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATa                                                                                                                       <  1:1413498/115‑1 (MQ=255)
                                                             cTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATa                                                                                                                       >  2:1413498/1‑115 (MQ=255)
                                                              tGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTATTACACCTGATGCGCTTTGTCACTCATTCAGATTCACAt                                                                                   >  1:431133/1‑150 (MQ=255)
                                                               gTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACAtt                                                                                  >  1:267032/1‑150 (MQ=255)
                                                                  cAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAg                                                                               <  1:313655/150‑1 (MQ=255)
                                                                  cAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAg                                                                               >  1:107029/1‑150 (MQ=255)
                                                                   aaTGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCGGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGc                                                                              <  2:1297035/150‑1 (MQ=255)
                                                                   aaTGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGc                                                                              >  1:609767/1‑150 (MQ=255)
                                                                   aaTGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGc                                                                              <  2:595883/150‑1 (MQ=255)
                                                                    aTGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGcc                                                                             >  2:794884/1‑150 (MQ=255)
                                                                    aTGCGATGGATATTGATGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGcc                                                                             >  2:1173666/1‑150 (MQ=255)
                                                                       cGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGCCTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCttt                                                                          >  2:687155/1‑150 (MQ=255)
                                                                               atTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATcc                                                                  <  1:881539/150‑1 (MQ=255)
                                                                               atTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATcc                                                                  <  1:884547/150‑1 (MQ=255)
                                                                                 tGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCAc                                                                <  1:1496931/150‑1 (MQ=255)
                                                                                  gCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACt                                                               >  2:793026/1‑150 (MQ=255)
                                                                                   cTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCCCCATCCACTa                                                              >  2:1429711/1‑150 (MQ=255)
                                                                                     gAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACa                                                            >  2:1226570/1‑150 (MQ=255)
                                                                                        aTCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGa                                                         >  2:1361143/1‑150 (MQ=255)
                                                                                         tCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGAt                                                        <  2:333580/150‑1 (MQ=255)
                                                                                         tCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGAt                                                        >  2:1071106/1‑150 (MQ=255)
                                                                                          ccGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATc                                                       >  1:490859/1‑150 (MQ=255)
                                                                                            gCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATcgc                                                     >  2:747630/1‑150 (MQ=255)
                                                                                             caacaaCGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATcgcg                                                    >  2:864134/1‑150 (MQ=255)
                                                                                             caacaaCGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATcgcg                                                    >  2:894774/1‑150 (MQ=255)
                                                                                               acaacGCGTCTCCGGCAGATCCGCAATTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGcc                                                  <  1:942499/150‑1 (MQ=255)
                                                                                                        ctcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAg                                         <  1:229105/150‑1 (MQ=255)
                                                                                                         tcCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAgg                                        >  1:714571/1‑150 (MQ=255)
                                                                                                              cAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCGCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTGTAGGCTTCg                                   <  2:1204354/150‑1 (MQ=255)
                                                                                                               aGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGt                                  >  2:620795/1‑150 (MQ=255)
                                                                                                                gATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTc                                 <  2:1064158/150‑1 (MQ=255)
                                                                                                                  tCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTccc                               >  1:1402918/1‑150 (MQ=255)
                                                                                                                   ccGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTcccc                              <  1:231999/150‑1 (MQ=255)
                                                                                                                          ttCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATAAACCTGATGCACTTTCTCACTAATTCAGATTCACATTAAGCCTTTCACCATACACTACAAGATAGAGCCACGCTTTAGGCTTCGTCCCCAGtttta                       >  2:949494/1‑149 (MQ=255)
                                                                                                                           tCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTcc                      <  2:1609807/150‑1 (MQ=255)
                                                                                                                            cATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCa                     <  1:793026/150‑1 (MQ=255)
                                                                                                                            cATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCa                     >  1:273589/1‑150 (MQ=255)
                                                                                                                             aTTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCAt                    <  2:863786/150‑1 (MQ=255)
                                                                                                                              ttAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGCTTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATc                   >  2:471208/1‑150 (MQ=255)
                                                                                                                                aaTGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCt                 >  1:1401598/1‑150 (MQ=255)
                                                                                                                                 aTGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCtt                >  1:671230/1‑150 (MQ=255)
                                                                                                                                   gAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCg              <  2:671230/150‑1 (MQ=255)
                                                                                                                                    aGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCATCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGc             <  1:781577/150‑1 (MQ=255)
                                                                                                                                    aGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGc             <  1:1071106/150‑1 (MQ=255)
                                                                                                                                     gTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCa            >  2:1335154/1‑150 (MQ=255)
                                                                                                                                      tATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTGTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAg           >  1:429543/1‑150 (MQ=255)
                                                                                                                                      aaaatCCGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAg           >  2:494376/5‑150 (MQ=255)
                                                                                                                                        tGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCt         >  1:533854/1‑150 (MQ=255)
                                                                                                                                         gTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtg        >  1:98573/1‑150 (MQ=255)
                                                                                                                                         ctCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTAACACATTCACATTAACATTCAGCCTATCACCATCAACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtg        >  2:1237783/2‑150 (MQ=255)
                                                                                                                                         atCCGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtg        >  1:529108/2‑150 (MQ=255)
                                                                                                                                           ccGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg      >  2:1506356/1‑150 (MQ=255)
                                                                                                                                           ccGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg      >  1:38070/1‑150 (MQ=255)
                                                                                                                                           ccGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg      >  1:1127369/1‑150 (MQ=255)
                                                                                                                                           ccGGCAGAGCCCTCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCtgtg      >  1:1121823/1‑150 (MQ=255)
                                                                                                                                            cGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGa     >  1:266507/1‑150 (MQ=255)
                                                                                                                                            cGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGa     >  1:1570029/1‑150 (MQ=255)
                                                                                                                                             ggCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAg    <  1:1226570/150‑1 (MQ=255)
                                                                                                                                               cAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTa  >  2:192273/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TAGCATCCTGGCGGAAGATGGCGAAATTGACATGCATTGTGATTACTGCGGTAACCACTATCTGTTCAATGCGATGGATATTGCTGAAATCCGCAACAACGCGTCTCCGGCAGATCCGCAAGTTCATTAATGAGTATGTCCGGCAGAGAGGGCTCTGCCGGATTTTACGACTAATACACCTGATGCGCTTTGTCACTCATTCAGATTCACATTCAGCCTTTCACCATCCACTACAAGATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTA  >  NC_000913/3530523‑3530815

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