Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,111,823 C→T 5.4% pseudogene (268/963 nt) yghE ← pseudogene, secretion pathway protein, L‑type protein homology; putative general secretion pathway for protein export (GSP)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,111,8230CT 5.4% 437.8 / 8.8 129pseudogene (268/963 nt)yghEpseudogene, secretion pathway protein, L‑type protein homology; putative general secretion pathway for protein export (GSP)
Reads supporting (aligned to +/- strand):  ref base C (61/61);  new base T (4/3);  total (65/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

GCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCAT  >  NC_000913/3111679‑3111970
                                                                                                                                                |                                                                                                                                                   
gCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGgcg                                                                                                                                                 >  2:1081821/1‑149 (MQ=255)
gCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGgcg                                                                                                                                                 >  4:1405367/1‑149 (MQ=255)
    agCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTcc                                                                                                                                             <  2:1405658/149‑1 (MQ=255)
    agCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTcc                                                                                                                                             <  1:206847/149‑1 (MQ=255)
         tAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATccc                                                                                                                                        >  1:1276378/1‑149 (MQ=255)
             ggTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGc                                                                                                                                    <  2:677383/149‑1 (MQ=255)
               tCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCAt                                                                                                                                  >  3:1296487/1‑149 (MQ=255)
                 aGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTc                                                                                                                                >  2:398491/1‑149 (MQ=255)
                   gCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCAt                                                                                                                              >  1:437546/1‑149 (MQ=255)
                    cGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATc                                                                                                                             <  4:201227/149‑1 (MQ=255)
                      ttCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGa                                                                                                                           >  4:677886/1‑149 (MQ=255)
                            gTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGc                                                                                                                     <  1:1081821/149‑1 (MQ=255)
                             tCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTc                                                                                                                                              <  4:82376/123‑1 (MQ=255)
                             tCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTc                                                                                                                                              >  3:82376/1‑123 (MQ=255)
                             tCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCa                                                                                                                    >  4:814862/1‑149 (MQ=255)
                                 cgcCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATa                                                                                                                <  4:1343525/149‑1 (MQ=255)
                                  gcCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGa                                                                                                                  >  3:1046696/1‑146 (MQ=255)
                                  gcCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGa                                                                                                                  <  4:1046696/146‑1 (MQ=255)
                                  gcCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATaa                                                                                                               <  2:1334709/149‑1 (MQ=255)
                                      gCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGCACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATg                                                                                                                         <  4:1429033/135‑1 (MQ=255)
                                      gCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATg                                                                                                                         >  3:1429033/1‑135 (MQ=255)
                                      gCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCg                                                                                                           >  4:1405624/1‑149 (MQ=255)
                                       cccATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGt                                                                                                          <  3:814862/149‑1 (MQ=255)
                                          aTGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTgc                                                                                                       <  3:307788/149‑1 (MQ=255)
                                          aTGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTgc                                                                                                       <  3:272230/149‑1 (MQ=255)
                                            gCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCGGGCTGGCGCtgg                                                                                                                                             <  1:185890/109‑3 (MQ=255)
                                            gCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCGGGCTGGCGCtgg                                                                                                                                             >  2:185890/1‑107 (MQ=255)
                                               gATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAAc                                                                                                              <  3:1381599/137‑1 (MQ=255)
                                                aTAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCATATAACCCGTTGCGCATCg                                                                                                 >  1:207201/1‑149 (MQ=255)
                                                aTAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCa                                                                                                    <  1:462145/146‑1 (MQ=255)
                                                    cTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAg                                                                                             <  4:1296487/149‑1 (MQ=255)
                                                     ttttCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATc                                                                                                  <  4:1200127/143‑1 (MQ=255)
                                                     ttttCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGc                                                                                            >  2:884427/1‑149 (MQ=255)
                                                     ttttCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGc                                                                                            >  3:1173397/1‑149 (MQ=255)
                                                     ttttCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCAGCGATGCGGCAGATAACCCGTTGCGCATc                                                                                                  >  3:1200127/1‑143 (MQ=255)
                                                            gTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCgg                                                                                                                                     >  2:326647/1‑101 (MQ=255)
                                                            gTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCgg                                                                                                                                     <  1:326647/101‑1 (MQ=255)
                                                            gTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCATATAACCCGTTGCGCATCGCAAGCCATAAAc                                                                                     <  2:207201/149‑1 (MQ=255)
                                                                 aTGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGTCAtt                                                                                                                                 <  4:700604/100‑1 (MQ=255)
                                                                 aTGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCAtt                                                                                                                                 >  3:700604/1‑100 (MQ=255)
                                                                  tGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACaa                                                                               <  3:1405367/149‑1 (MQ=255)
                                                                    ccGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAgg                                                                             >  4:396793/1‑149 (MQ=255)
                                                                              ttAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGc                                                                                            >  3:177383/1‑124 (MQ=255)
                                                                              ttAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATc                                                                   >  1:36201/1‑149 (MQ=255)
                                                                              ttAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATc                                                                   <  4:1173397/149‑1 (MQ=255)
                                                                                aGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCa                                                                 <  2:437546/149‑1 (MQ=255)
                                                                                aGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCa                                                                 >  2:498802/1‑149 (MQ=255)
                                                                                   gAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTc                                                              >  3:651716/1‑149 (MQ=255)
                                                                                    aaaGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGtt                                                                                                         >  1:381701/1‑105 (MQ=255)
                                                                                       gCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACaa                                                                               <  3:531069/128‑1 (MQ=255)
                                                                                        cGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGgcg                                                         >  2:1271316/1‑149 (MQ=255)
                                                                                         ggTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATcc                                                                  >  2:507011/1‑139 (MQ=255)
                                                                                         ggTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATcc                                                                  <  1:507011/139‑1 (MQ=255)
                                                                                          gTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGcgcg                                                       >  2:115737/1‑149 (MQ=255)
                                                                                           tAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATg                                                                         <  1:1000693/130‑1 (MQ=255)
                                                                                           tAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATg                                                                         >  2:1000693/1‑130 (MQ=255)
                                                                                             aaTCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAAt                                                                          <  3:961987/127‑1 (MQ=255)
                                                                                              aTCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCa                                                                                                                                   >  1:954352/1‑69 (MQ=255)
                                                                                              aTCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCa                                                                                                                                   <  2:954352/69‑1 (MQ=255)
                                                                                                ccTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACaa                                                                               >  2:1244393/1‑119 (MQ=255)
                                                                                                ccTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACaa                                                                               >  4:679381/1‑119 (MQ=255)
                                                                                                ccTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACaa                                                                               <  3:679381/119‑1 (MQ=255)
                                                                                                ccTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACaa                                                                               <  1:1244393/119‑1 (MQ=255)
                                                                                                 cTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAg                                                >  2:1174514/1‑149 (MQ=255)
                                                                                                 cTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAg                                                <  4:1148152/149‑1 (MQ=255)
                                                                                                    cTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAgcgc                                             >  1:1328261/1‑149 (MQ=255)
                                                                                                      cATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCgg                                           <  3:719346/149‑1 (MQ=255)
                                                                                                        ttGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGt                                         <  3:1142415/149‑1 (MQ=255)
                                                                                                         tGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTg                                        <  2:1328261/149‑1 (MQ=255)
                                                                                                           aGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGAt                                      >  4:260711/1‑149 (MQ=255)
                                                                                                           aGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGAt                                      <  1:115737/149‑1 (MQ=255)
                                                                                                            gCTGGGCGCGGATCATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAACGAATGGACATCCAGTCCGGCGCGACGGTAGCGCggg                                          <  4:1236771/144‑1 (MQ=255)
                                                                                                            gCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCggg                                          >  3:1236771/1‑144 (MQ=255)
                                                                                                             cTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGACTGGATGTCCATTcct                                                            <  1:1391207/125‑2 (MQ=255)
                                                                                                                 gcgcGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGtt                                                                                                         <  1:998168/76‑1 (MQ=255)
                                                                                                                 gcgcGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGtt                                                                                                         >  2:998168/1‑76 (MQ=255)
                                                                                                                    cgGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTc                             <  4:651716/149‑1 (MQ=255)
                                                                                                                      gAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAAt                                                                          >  1:507227/1‑102 (MQ=255)
                                                                                                                      gAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAAt                                                                          <  2:507227/102‑1 (MQ=255)
                                                                                                                      gAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCt                           <  3:260711/149‑1 (MQ=255)
                                                                                                                       aaCATGACCGGAGCCAGCGCCAGCCCTGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTg                          >  2:923868/1‑149 (MQ=255)
                                                                                                                         cATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTcc                            <  4:194470/145‑1 (MQ=255)
                                                                                                                            gACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCca                                                                                  <  3:1436253/88‑1 (MQ=255)
                                                                                                                            gACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCca                                                                                  >  4:1436253/1‑88 (MQ=255)
                                                                                                                            gACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAAcc                     >  3:538273/1‑149 (MQ=255)
                                                                                                                             aCCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCa                    >  3:470786/1‑149 (MQ=255)
                                                                                                                                    ggagcgccgggCTGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAAc                                                                                     <  4:773443/75‑1 (MQ=255)
                                                                                                                                    ggagcgccgggCTGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAAc                                                                                     >  3:773443/3‑77 (MQ=255)
                                                                                                                                    ccagcgccagcCCGGCGCTCCATCCCCGGCACTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCGAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCa             >  4:899300/1‑149 (MQ=255)
                                                                                                                                      agcgccgggCTGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGttt                               <  2:385737/129‑1 (MQ=255)
                                                                                                                                      agcgccgggCTGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGttt                               >  1:385737/1‑129 (MQ=255)
                                                                                                                                      agcgccagcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCaac           >  3:87703/1‑148 (MQ=255)
                                                                                                                                      agcgccagcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAAt           >  1:469109/1‑149 (MQ=255)
                                                                                                                                       gcgccagcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATg          >  4:37254/1‑149 (MQ=255)
                                                                                                                                         gccagcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGtt        <  4:470786/149‑1 (MQ=255)
                                                                                                                                         gccagcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGtt        <  2:36201/149‑1 (MQ=255)
                                                                                                                                           cagcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCa                                                                                          >  3:1218718/1‑65 (MQ=255)
                                                                                                                                             gcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGAc                                                     >  3:943453/1‑100 (MQ=255)
                                                                                                                                             gcCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGAc                                                     <  4:943453/100‑1 (MQ=255)
                                                                                                                                               ccGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCag  <  4:87703/149‑2 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCAGAGCTGTAACGGTCAGGCGTTCAGCGTCTGCGCCAGCCCATGCTGATAGCTTTTCCGGTGCAATGCCGACTACGGTTAGCGAAAGCGGTAAATCCTGCTCATTGAGCTGGGCGCGGAACATGACCGGAGCCAGCGCCAGCCCGGCGCTCCATCCCCGGCATTCATCGATGCGGCAGATAACCCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCAT  >  NC_000913/3111679‑3111970

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: