Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,111,962 G→A 15.1% pseudogene (129/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,9620GA15.1% 359.4 / 44.4 144pseudogene (129/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 G (50/72);  new base A (7/15);  total (57/87)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 146 total aligned reads displayed.

CCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCT  >  NC_000913/3111863‑3112061
                                                                                                   |                                                                                                   
     gcgcATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTgcg                                                                                                >  2:2105002/1‑100 (MQ=255)
      cgcATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTgcgc                                                                                               >  2:2688782/1‑100 (MQ=255)
      cgcATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTgcgc                                                                                               <  2:545697/100‑1 (MQ=255)
      cgcATCGCAAGCCATAAACCAACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTgcgc                                                                                               >  1:1909558/1‑100 (MQ=255)
       gcATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTgcgc                                                                                               >  1:2877378/1‑99 (MQ=255)
       gcATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCa                                                                                              >  1:40368/1‑100 (MQ=255)
       gcATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCa                                                                                              >  2:1875833/1‑100 (MQ=255)
       gcATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCa                                                                                              <  1:2285767/100‑1 (MQ=255)
        cATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCag                                                                                             <  2:289077/100‑2 (MQ=255)
         atcacAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagt                                                                                            <  2:106863/96‑3 (MQ=255)
         aTCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCata                                                                                            <  2:872809/100‑1 (MQ=255)
         aTCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCata                                                                                            <  2:2477611/100‑1 (MQ=255)
         aTCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagt                                                                                            <  1:1172948/100‑3 (MQ=255)
          tCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCatat                                                                                           >  2:1672983/1‑100 (MQ=255)
          tCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCata                                                                                            >  2:950192/1‑99 (MQ=255)
          tCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagtg                                                                                           <  2:1490812/100‑4 (MQ=255)
           cGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagtgg                                                                                          <  2:1929259/100‑5 (MQ=255)
            gCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGc                                                                                         >  2:161704/1‑100 (MQ=255)
            gCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGc                                                                                         >  1:2716709/1‑100 (MQ=255)
            gCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagtggt                                                                                         <  1:743121/100‑6 (MQ=255)
            gCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagtggt                                                                                         <  1:2421094/100‑6 (MQ=255)
             cAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCt                                                                                        >  1:605300/1‑100 (MQ=255)
             cAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAACATTGCGCagtggtt                                                                                        <  1:1889227/100‑7 (MQ=255)
              aaGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTg                                                                                       >  2:1007048/1‑100 (MQ=255)
                gCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGtt                                                                                     <  2:809048/100‑1 (MQ=255)
                 ccATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTg                                                                                    <  2:2918697/100‑1 (MQ=255)
                 ccATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTg                                                                                    <  1:1138328/100‑1 (MQ=255)
                  cATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTTCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTgg                                                                                   >  2:550950/1‑100 (MQ=255)
                    tAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTg                                                                                 <  1:2030259/100‑1 (MQ=255)
                    tAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTg                                                                                 <  2:1824361/100‑1 (MQ=255)
                      aaCCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGaa                                                                               <  2:946475/100‑1 (MQ=255)
                      aaCCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGaa                                                                               <  2:2669320/100‑1 (MQ=255)
                       aCCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAAt                                                                              <  2:2183291/100‑1 (MQ=255)
                       aCCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAAt                                                                              <  1:2305513/100‑1 (MQ=255)
                        cccACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATc                                                                             <  2:644497/100‑1 (MQ=255)
                              aGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCt                                                                       <  1:1128282/100‑1 (MQ=255)
                              aGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCt                                                                       <  1:1929910/100‑1 (MQ=255)
                               ggAATGGATATCCATTTCGGCTCGACGGTAGTGCAGGTGATGCCGGTTTCCTGCAACCACTGCGCAATGTTGCCCATATGCTGTTGGTGAATCACt                                                                          <  1:1792362/96‑1 (MQ=255)
                               ggAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCt                                                                       >  1:2296793/1‑99 (MQ=255)
                                 aaTGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACg                                                                    <  1:1470550/100‑1 (MQ=255)
                                  aTGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACgg                                                                   <  2:213355/100‑1 (MQ=255)
                                    ggACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGtt                                                                 >  1:2160431/1‑100 (MQ=255)
                                     gACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTg                                                                <  2:2511025/100‑1 (MQ=255)
                                     gACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTg                                                                >  2:2137432/1‑100 (MQ=255)
                                       cATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGcc                                                              <  1:394010/100‑1 (MQ=255)
                                       cATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGcc                                                              <  2:751517/100‑1 (MQ=255)
                                       cATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGcc                                                              <  2:2907151/100‑1 (MQ=255)
                                         tCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAg                                                            <  1:1095925/100‑1 (MQ=255)
                                         tCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAg                                                            >  1:1279584/1‑100 (MQ=255)
                                         tCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAg                                                            >  1:2726536/1‑100 (MQ=255)
                                         tCCAGTCCGGCGCGACGGTAGCACGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAg                                                            >  1:1693556/1‑100 (MQ=255)
                                          ccAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGt                                                           >  1:2130939/1‑100 (MQ=255)
                                           cAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGtt                                                          <  2:529867/100‑1 (MQ=255)
                                           cAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGtt                                                          <  2:291434/100‑1 (MQ=255)
                                           cAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGtt                                                          >  1:2291189/1‑100 (MQ=255)
                                              tCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCt                                                       <  1:1139225/100‑1 (MQ=255)
                                              tCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCt                                                       >  2:921433/1‑100 (MQ=255)
                                                cGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGc                                                        <  1:627668/97‑1 (MQ=255)
                                                 ggCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGt                                                    >  2:548183/1‑100 (MQ=255)
                                                  gcgcGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTc                                                   <  2:605303/100‑1 (MQ=255)
                                                  gcgcGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTc                                                   <  2:2588567/100‑1 (MQ=255)
                                                   cgcgACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCg                                                  >  1:943657/1‑100 (MQ=255)
                                                    gcgACGGTAGCGCGGGTGATTCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGGTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCg                                                  >  1:805120/1‑99 (MQ=255)
                                                     cgACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGat                                                <  1:519859/100‑1 (MQ=255)
                                                        cGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATAtt                                             >  2:545803/1‑100 (MQ=255)
                                                         ggTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTc                                            <  2:2096186/100‑1 (MQ=255)
                                                          gTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCa                                           <  1:303730/100‑1 (MQ=255)
                                                          gTAGCGCGGGCGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCa                                           <  1:303060/100‑1 (MQ=255)
                                                              cgcgGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACgg                                       >  1:1048569/1‑100 (MQ=255)
                                                               gcgGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACggg                                      <  1:1181131/100‑1 (MQ=255)
                                                                 gggTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGc                                    <  2:1952418/100‑1 (MQ=255)
                                                                  ggTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCg                                   <  2:2712545/100‑1 (MQ=255)
                                                                     gATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGaaa                                >  2:2352680/1‑100 (MQ=255)
                                                                      aTGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGTGGCGAAAt                               >  1:2901237/1‑100 (MQ=255)
                                                                      aTGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAAt                               <  2:1033916/100‑1 (MQ=255)
                                                                          cGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAg                           <  2:1693563/100‑1 (MQ=255)
                                                                           gTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGt                          <  2:1433857/100‑1 (MQ=255)
                                                                           gTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGt                          <  2:2742154/100‑1 (MQ=255)
                                                                               tCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCAt                      >  1:2083850/1‑100 (MQ=255)
                                                                                ccTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATc                     <  2:2100931/100‑1 (MQ=255)
                                                                                ccTGCAACCACCGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATc                     <  2:394272/100‑1 (MQ=255)
                                                                                 cTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGTTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCg                    <  1:2565610/100‑1 (MQ=255)
                                                                                  tGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGa                   >  2:2303067/1‑100 (MQ=255)
                                                                                  tGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGa                   >  1:2421825/1‑100 (MQ=255)
                                                                                   acagcatatGCGCAACATTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGat                  >  2:1662205/9‑100 (MQ=255)
                                                                                    cAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGata                 <  2:2834555/100‑1 (MQ=255)
                                                                                      gcatatGCGCAACATTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATc               >  2:2154257/6‑100 (MQ=255)
                                                                                      gcatatGCGCAACATTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATc               >  1:379471/6‑100 (MQ=255)
                                                                                      gcatatGCGCAACATTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATc               >  1:760547/6‑100 (MQ=255)
                                                                                        cACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCt             <  2:2291204/100‑1 (MQ=255)
                                                                                        atatGCGCAACATTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCt             >  1:2558584/4‑100 (MQ=255)
                                                                                          cTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTgg           <  2:302739/100‑1 (MQ=255)
                                                                                            gcgcAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGtt         >  2:2169779/1‑100 (MQ=255)
                                                                                             cgcAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTc        <  1:1778029/100‑1 (MQ=255)
                                                                                              gcAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCa       >  2:1268393/1‑100 (MQ=255)
                                                                                               cAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACCGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAg      >  1:1554224/1‑100 (MQ=255)
                                                                                                aaTGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGc     <  1:2888488/100‑1 (MQ=255)
                                                                                                aaTGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGc     <  1:2794967/100‑1 (MQ=255)
                                                                                                aaTGTTGCGCATATACTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGc     <  1:2368024/100‑1 (MQ=255)
                                                                                                   gTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGctct  >  2:2748440/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CCGTTGCGCATCGCAAGCCATAAACCCACAAGGAATGGACATCCAGTCCGGCGCGACGGTAGCGCGGGTGATGCCGTTTTCCTGCAACCACTGCGCAATGTTGCGCATATGCTGTTGGTGAATCACTGCTACGGTTGCCAGTTGCTGGTCGATATTCAACGGGGCGAAATGCAGTTCATCGATATCCTGGTTCAGCTCT  >  NC_000913/3111863‑3112061

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