Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,354,220 C→A 5.2% G100V (GGT→GTT)  sapC ← antimicrobial peptide transport ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,354,2200CA 5.2% 472.1 / 11.4 134G100V (GGT→GTT) sapCantimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (61/66);  new base A (2/5);  total (63/71)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 142 total aligned reads displayed.

GATAATTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAG  >  NC_000913/1354072‑1354367
                                                                                                                                                    |                                                                                                                                                   
gATAATTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGccac                                                                                                                                                     <  1:349209/149‑1 (MQ=255)
   aaTTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGccaccca                                                                                                                                                  <  1:943144/149‑1 (MQ=255)
     ttGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACg                                                                                                                                                <  3:764695/149‑1 (MQ=255)
        ccAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGtt                                                                                                                                             <  4:1143602/149‑1 (MQ=255)
         cAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTg                                                                                                                                            <  2:485905/149‑1 (MQ=255)
          aaCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTg                                                                                                                                            >  1:1052659/1‑148 (MQ=255)
               cagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCat                                                                                                                                      <  2:539246/149‑2 (MQ=255)
               cagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCat                                                                                                                                      <  3:728396/149‑2 (MQ=255)
                agcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGc                                                                                                                                     <  3:1276800/149‑1 (MQ=255)
                agcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGc                                                                                                                                     >  1:81331/1‑149 (MQ=255)
                agcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCatt                                                                                                                                     <  1:703163/149‑3 (MQ=255)
                  cGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTc                                                                                                                                   <  3:693077/149‑1 (MQ=255)
                        gATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCa                                                                                                                             <  4:359347/149‑1 (MQ=255)
                        gATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCa                                                                                                                             >  1:645945/1‑149 (MQ=255)
                        gATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCa                                                                                                                             >  2:1055591/1‑149 (MQ=255)
                         aTCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAg                                                                                                                            >  1:1555632/1‑149 (MQ=255)
                                gTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACgg                                                                                                                     >  2:1304455/1‑149 (MQ=255)
                                   aGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCtt                                                                                                                  <  2:362038/149‑1 (MQ=255)
                                        tCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCattt                                                                                                                                    <  2:1171223/126‑4 (MQ=255)
                                          cAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACAt                                                                                                           <  1:185873/149‑1 (MQ=255)
                                           aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATc                                                                                                          <  1:1278487/149‑1 (MQ=255)
                                           aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATc                                                                                                          >  3:54961/1‑149 (MQ=255)
                                            gAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATcg                                                                                                         >  4:143285/1‑149 (MQ=255)
                                             aaTATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATcgc                                                                                                        <  3:717975/149‑1 (MQ=255)
                                                 tGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTcc                                                                                                    >  2:248369/1‑149 (MQ=255)
                                                  ggTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCg                                                                                                   <  2:787982/149‑1 (MQ=255)
                                                   gTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGa                                                                                                  <  2:1403702/149‑1 (MQ=255)
                                                    ttAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAg                                                                                                 >  3:568013/1‑149 (MQ=255)
                                                          aCCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgtc                                                                                           <  2:1441374/149‑1 (MQ=255)
                                                             gCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGt                                                                                        >  4:348677/1‑149 (MQ=255)
                                                                   cGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAg                                                                                  <  4:880542/149‑1 (MQ=255)
                                                                   cGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAg                                                                                  <  4:78940/149‑1 (MQ=255)
                                                                   cGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAg                                                                                  <  4:1203358/149‑1 (MQ=255)
                                                                      aGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGgaa                                                                               >  3:45328/1‑149 (MQ=255)
                                                                      aGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGgaa                                                                               <  1:124803/149‑1 (MQ=255)
                                                                       ggCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGgaag                                                                              >  4:1419263/1‑149 (MQ=255)
                                                                        gCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACa                                                                                                            <  3:765590/118‑1 (MQ=255)
                                                                         ccGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGgaagaa                                                                            <  2:615180/149‑1 (MQ=255)
                                                                          cGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGaagaaa                                                                           <  2:1560023/149‑1 (MQ=255)
                                                                          cGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGaagaaa                                                                           <  3:143285/149‑1 (MQ=255)
                                                                          cGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGaagaaa                                                                           <  4:1088284/149‑1 (MQ=255)
                                                                          cGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGaagaaa                                                                           <  3:498967/149‑1 (MQ=255)
                                                                          cGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGaagaaa                                                                           >  4:648925/1‑149 (MQ=255)
                                                                           gtgtgtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGaagaaag                                                                          >  2:62561/1‑149 (MQ=255)
                                                                           gtgtgtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTAAGTCCCCAGGaagaaa                                                                           <  1:684328/148‑1 (MQ=255)
                                                                             gtgtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAagaaagaa                                                                        <  2:1555632/149‑1 (MQ=255)
                                                                              tgtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAgaaagaaa                                                                       <  4:721110/149‑1 (MQ=255)
                                                                                tCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACt                                                                     <  3:1419263/149‑1 (MQ=255)
                                                                                  gCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTc                                                                   >  2:1074975/1‑149 (MQ=255)
                                                                                   ccccGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAgaaagaaa                                                                       >  3:953936/1‑144 (MQ=255)
                                                                                   ccccGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAgaaagaaa                                                                       <  4:953936/144‑1 (MQ=255)
                                                                                   ccccGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCg                                                                  >  3:458622/1‑149 (MQ=255)
                                                                                   ccccGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCCAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCg                                                                  >  4:94933/1‑149 (MQ=255)
                                                                                     ccGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGc                                                                 <  1:502659/148‑1 (MQ=255)
                                                                                     ccGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGc                                                                 >  2:502659/1‑148 (MQ=255)
                                                                                        gCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATa                                                             <  1:248369/149‑1 (MQ=255)
                                                                                                  ccAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCt                                                                                                                   <  1:1519010/85‑1 (MQ=255)
                                                                                                  ccAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCa                                                                                   <  3:11701/117‑1 (MQ=255)
                                                                                                  ccAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCa                                                                                   >  4:11701/1‑117 (MQ=255)
                                                                                                   cAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATg                                                  >  1:901211/1‑149 (MQ=255)
                                                                                                        ccAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGAcggc                                             <  3:980720/149‑1 (MQ=255)
                                                                                                        ccAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGAcggc                                             >  4:186798/1‑149 (MQ=255)
                                                                                                          aGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATc                                                                                                          >  4:582306/1‑86 (MQ=255)
                                                                                                          aGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATc                                                                                                          <  3:582306/86‑1 (MQ=255)
                                                                                                          aGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGAcggcgg                                           >  2:796379/1‑149 (MQ=255)
                                                                                                             ccACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaa                                        >  1:1511492/1‑149 (MQ=255)
                                                                                                             ccACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaa                                        >  2:66358/1‑149 (MQ=255)
                                                                                                              cacaGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCa                                                                                          <  2:804104/98‑1 (MQ=255)
                                                                                                                 aGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaataat                                    <  3:306739/149‑1 (MQ=255)
                                                                                                                 aGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaataat                                    >  1:1010320/1‑149 (MQ=255)
                                                                                                                   aTCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGaagaaaga                                                                         <  3:871301/110‑1 (MQ=255)
                                                                                                                      gtcgCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGTTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAgaaagaaa                                                                       <  3:896718/109‑1 (MQ=255)
                                                                                                                      gtcgCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAgaaagaaa                                                                       >  4:896718/1‑109 (MQ=255)
                                                                                                                      gtcgCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATa                               >  2:1486723/1‑149 (MQ=255)
                                                                                                                       tcgCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATaa                              <  1:298949/149‑1 (MQ=255)
                                                                                                                        cgcgGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAAc                             <  1:1469236/149‑1 (MQ=255)
                                                                                                                        cgcgGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAAc                             <  2:1032105/149‑1 (MQ=255)
                                                                                                                         gcgGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAAcc                            <  4:831924/149‑1 (MQ=255)
                                                                                                                             cAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCGGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGa                                                 <  4:545060/124‑1 (MQ=255)
                                                                                                                             cAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCGGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGa                                                 >  3:545060/1‑124 (MQ=255)
                                                                                                                              aaGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaataat                                    <  2:930285/136‑1 (MQ=255)
                                                                                                                              aaGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGaa                       >  3:1543859/1‑149 (MQ=255)
                                                                                                                               aGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGaaa                      <  4:45328/149‑1 (MQ=255)
                                                                                                                               aGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGaaa                      <  4:119444/149‑1 (MQ=255)
                                                                                                                                 cGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAAtt                    <  4:867899/149‑1 (MQ=255)
                                                                                                                                  gTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAg                                                            >  2:1161938/1‑108 (MQ=255)
                                                                                                                                  gTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAg                                                            <  1:1161938/108‑1 (MQ=255)
                                                                                                                                   tCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGc                  >  2:582883/1‑149 (MQ=255)
                                                                                                                                   tCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGc                  <  2:645945/149‑1 (MQ=255)
                                                                                                                                     accacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTg                >  2:486918/1‑149 (MQ=255)
                                                                                                                                     accacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTg                >  1:1144462/1‑149 (MQ=255)
                                                                                                                                      ccacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATa                               >  4:221053/1‑133 (MQ=255)
                                                                                                                                      acacaAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAAGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATa                               <  3:221053/132‑1 (MQ=255)
                                                                                                                                          aaaTGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGa           >  2:341083/1‑149 (MQ=255)
                                                                                                                                          aaaTGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGa           <  2:81331/149‑1 (MQ=255)
                                                                                                                                              gcgcCAACCGTGGGTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGGCGGCGGCAATAATTGATAACCGAGAAAt                     >  2:1348455/1‑135 (MQ=255)
                                                                                                                                                gcCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGt     >  3:871297/1‑149 (MQ=255)
                                                                                                                                                   acccacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAgaaagaaa                                                                       <  1:384535/80‑1 (MQ=255)
                                                                                                                                                   acccacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAgaaagaaa                                                                       <  1:1102874/80‑1 (MQ=255)
                                                                                                                                                   acccacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAgaaagaaa                                                                       >  2:1102874/1‑80 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATAATTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAG  >  NC_000913/1354072‑1354367

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