Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,800 C→A 7.3% A580S (GCG→TCG)  yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,8000CA 7.3% 1291.5 / 66.1 370A580S (GCG→TCG) yicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (165/178);  new base A (12/15);  total (177/193)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

Only 100 of 372 total aligned reads displayed.

TAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCA  >  NC_000913/3832652‑3832948
                                                                                                                                                    |                                                                                                                                                    
 aGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCa                                                                                                                                                     >  4:266176/1‑149 (MQ=255)
       gCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCa                                                                                                                                               >  1:452758/1‑149 (MQ=255)
        cacaTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcgg                                                                                                                                              <  2:417272/149‑3 (MQ=255)
        cacaTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcgg                                                                                                                                              <  1:305311/149‑3 (MQ=255)
         acaTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcggg                                                                                                                                             <  2:451246/149‑4 (MQ=255)
           aTCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCtt                                                                                                                                           >  3:533605/1‑149 (MQ=255)
            tCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTc                                                                                                                                          >  1:386913/1‑149 (MQ=255)
            tCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTc                                                                                                                                          >  1:383095/1‑149 (MQ=255)
               cccGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCa                                                                                                                                         <  1:100497/147‑1 (MQ=255)
                   cTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAATCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATAc                                                                                                                                   <  4:347584/149‑1 (MQ=255)
                    ttCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACa                                                                                                                                  >  3:381305/1‑149 (MQ=255)
                     tCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAg                                                                                                                                 <  1:596704/149‑1 (MQ=255)
                      cAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGa                                                                                                                                <  4:310951/149‑1 (MQ=255)
                       aGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGAt                                                                                                                               <  4:374413/149‑1 (MQ=255)
                       aGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGAt                                                                                                                               >  3:333332/1‑149 (MQ=255)
                        gTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGCAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATa                                                                                                                              >  3:148236/1‑149 (MQ=255)
                         tGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATAc                                                                                                                             >  3:557854/1‑149 (MQ=255)
                         tGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATAc                                                                                                                             <  4:508397/149‑1 (MQ=255)
                              acCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGcat                                                                                                                        >  4:486162/1‑149 (MQ=255)
                                cGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCa                                                                                                                                         >  3:33970/1‑130 (MQ=255)
                                     cAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGCGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCgg                                                                                                                 <  2:255983/149‑1 (MQ=255)
                                          aTCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCAttt                                                                                                            >  2:174443/1‑149 (MQ=255)
                                           tCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                                                                           >  3:414904/1‑149 (MQ=255)
                                           tCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                                                                           <  3:100426/149‑1 (MQ=255)
                                           tCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                                                                           <  4:239822/149‑1 (MQ=255)
                                            cACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCa                                                                                                          <  4:331036/149‑1 (MQ=255)
                                               gTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGtt                                                                                                       <  1:280307/149‑1 (MQ=255)
                                                 tGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGc                                                                                                     <  1:488964/149‑1 (MQ=255)
                                                 tGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGc                                                                                                     <  4:433083/149‑1 (MQ=255)
                                                     gCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcgg                                                                                                                                              >  4:558726/1‑102 (MQ=255)
                                                       cTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaag                                                                                               <  4:541274/149‑1 (MQ=255)
                                                         aaCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                                                                             <  3:548030/149‑1 (MQ=255)
                                                         aaCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                                                                             <  3:170680/149‑1 (MQ=255)
                                                          aCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGaagaag                                                                                            <  1:551737/149‑1 (MQ=255)
                                                            aTGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAgcg                                                                                          >  2:48089/1‑149 (MQ=255)
                                                              gTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGca                                                                                        <  1:583935/149‑1 (MQ=255)
                                                                 ttGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGcacca                                                                                     >  2:433301/1‑149 (MQ=255)
                                                                  tGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgc                                                                                                                                                  <  2:359593/87‑1 (MQ=255)
                                                                  tGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgc                                                                                                                                                  >  1:359593/1‑87 (MQ=255)
                                                                   gACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACcaca                                                                                   <  3:487121/149‑1 (MQ=255)
                                                                   gACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACcaca                                                                                   >  3:561816/1‑149 (MQ=255)
                                                                     cGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATc                                                                                 <  1:509200/149‑1 (MQ=255)
                                                                       gTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                                                                               >  1:196301/1‑149 (MQ=255)
                                                                       gTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                                                                               <  4:60128/149‑1 (MQ=255)
                                                                        tCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcaca                                                                              >  2:94555/1‑149 (MQ=255)
                                                                         cAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAg                                                                             <  2:278005/149‑1 (MQ=255)
                                                                           aGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGa                                                                           <  1:132070/149‑1 (MQ=255)
                                                                                aTCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcat                                                                      <  1:107522/149‑1 (MQ=255)
                                                                                  cacaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCA                                                                    <  4:474154‑M1/149‑1 (MQ=255)
                                                                                   acaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCAT                                                                   <  4:146919‑M1/149‑1 (MQ=255)
                                                                                   acaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCAT                                                                   <  1:259014‑M1/149‑1 (MQ=255)
                                                                                    caAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATC                                                                  <  3:213487‑M1/149‑1 (MQ=255)
                                                                                     aaGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCG                                                                 >  3:440046‑M1/1‑149 (MQ=255)
                                                                                       gCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa                                                               <  1:94555/149‑1 (MQ=255)
                                                                                        ccGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAg                                                              <  3:95876/149‑1 (MQ=255)
                                                                                        ccGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAg                                                              >  4:495668/1‑149 (MQ=255)
                                                                                           gATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGcc                                                           >  3:244679/1‑149 (MQ=255)
                                                                                            aTCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGccc                                                          >  2:185820/1‑149 (MQ=255)
                                                                                              cGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCAc                                                        <  2:312679/149‑1 (MQ=255)
                                                                                                tCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACgg                                                      >  1:99796/1‑149 (MQ=255)
                                                                                                   gggAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGcac                                                   >  4:446064/1‑149 (MQ=255)
                                                                                                         tCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATa                                             >  1:238163/1‑149 (MQ=255)
                                                                                                          tcatcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATaa                                            >  3:573056/2‑149 (MQ=255)
                                                                                                           catcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtca                                                                       <  3:308646/121‑1 (MQ=255)
                                                                                                           catcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAg                                           <  3:189989/149‑1 (MQ=255)
                                                                                                            atcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCCCGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGa                                          >  4:361383/1‑149 (MQ=255)
                                                                                                            atcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGcatca                                                                                                                      <  3:514818/73‑1 (MQ=255)
                                                                                                            atcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGt                                                                                                   >  4:291696/1‑92 (MQ=255)
                                                                                                            atcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtca                                                                       >  2:103994/1‑120 (MQ=255)
                                                                                                             tcatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGAt                                         >  3:461226/1‑149 (MQ=255)
                                                                                                              catcatGCCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGAtt                                        <  1:330544/149‑1 (MQ=255)
                                                                                                               atcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGAttt                                       <  4:158037/149‑1 (MQ=255)
                                                                                                                tcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTg                                      <  2:123119/149‑1 (MQ=255)
                                                                                                                tcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTg                                      >  4:117874/1‑149 (MQ=255)
                                                                                                                tcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTg                                      >  4:137507/1‑149 (MQ=255)
                                                                                                                tcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTg                                      >  2:558919/1‑149 (MQ=255)
                                                                                                                 catGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGc                                     <  4:282937/149‑1 (MQ=255)
                                                                                                                  atGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGca                                                                                        >  1:461321/1‑97 (MQ=255)
                                                                                                                  atGGCCCGCATCATCGGCGTACCCCGCGCGTTCACACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCACCGTAGGCCCACGGCACACGATAAGATTTGCt                                    >  1:169747/1‑149 (MQ=255)
                                                                                                                      cccGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCg                                >  3:285755/1‑149 (MQ=255)
                                                                                                                      accGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCg                                >  3:135497/2‑149 (MQ=255)
                                                                                                                       ccGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgt                               >  3:271083/1‑149 (MQ=255)
                                                                                                                        cGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgtg                              <  4:519693/149‑1 (MQ=255)
                                                                                                                                    gTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATg                  >  1:15679/1‑149 (MQ=255)
                                                                                                                                     tACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGc                                     <  4:412246/129‑1 (MQ=255)
                                                                                                                                        gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTGAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTg              >  2:511893/4‑149 (MQ=255)
                                                                                                                                        cccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGt                                                                                                   >  1:86721/1‑64 (MQ=255)
                                                                                                                                        cccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGt                                                                                                   <  2:86721/64‑1 (MQ=255)
                                                                                                                                        cccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTg              <  2:15679/149‑1 (MQ=255)
                                                                                                                                        cccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTg              <  1:420976/149‑1 (MQ=255)
                                                                                                                                        cccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTg              >  2:498056/1‑149 (MQ=255)
                                                                                                                                         ctgcgcgtgCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACa                                                  >  2:542867‑M2/3‑111 (MQ=255)
                                                                                                                                         ccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTgg             >  4:308415/1‑149 (MQ=255)
                                                                                                                                         ccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTgg             >  3:236976/1‑149 (MQ=255)
                                                                                                                                          tgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACt                                                                         >  1:215290/2‑88 (MQ=255)
                                                                                                                                           gcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCgt                               >  2:62309/1‑129 (MQ=255)
                                                                                                                                                 tCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaa     <  2:203893/149‑1 (MQ=255)
                                                                                                                                                 tCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaa     <  2:242657/149‑1 (MQ=255)
                                                                                                                                                 tCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaa     <  2:386816/149‑1 (MQ=255)
                                                                                                                                                 tCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGcaaa     <  4:533605/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAGAACTGCACATCGCCCGCTTCAGTGAACACCGGCGCAACCATCACGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCACACGATAAGATTTGCTACCGTGTAAACGGCTATGGCTGGAGAGCAAACCA  >  NC_000913/3832652‑3832948

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