Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 632,017 (G)7→6 intergenic (+18/‑165) cstA → / → ybdD carbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913632,0110G.97.9% 427.2 / 2.1 94intergenic (+12/‑171)cstA/ybdDcarbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (37/55);  total (39/55)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 106 total aligned reads displayed.

CCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTG  >  NC_000913/631913‑632105
                                                                                                  |                                                                                              
ccGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAgg                                                                                               <  2:2584017/100‑1 (MQ=255)
 cGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggg                                                                                              >  2:38022/1‑100 (MQ=255)
 cGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggg                                                                                              >  1:2633286/1‑100 (MQ=255)
 cGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggg                                                                                              <  1:2412247/100‑1 (MQ=255)
  gAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAgggg                                                                                             >  2:1266781/1‑100 (MQ=255)
  gAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAgggg                                                                                             <  2:38368/100‑1 (MQ=255)
    aGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAgggggg                                                                                           <  2:246960/100‑1 (MQ=255)
     gCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggggggn                                                                                          <  1:1406343/100‑2 (MQ=255)
       cAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggggggatg                                                                                        >  2:797813/1‑97 (MQ=255)
       cAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggggggatg                                                                                        >  1:2754150/1‑97 (MQ=255)
       cAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggggggatg                                                                                        >  2:2656370/1‑97 (MQ=255)
        aaCGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggcgggatg                                                                                        >  2:527336/1‑96 (MQ=255)
            gCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGc                                                                                  <  2:689518/100‑1 (MQ=255)
            gCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGc                                                                                  >  2:2661512/1‑100 (MQ=255)
            gCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGc                                                                                  >  2:2139542/1‑100 (MQ=255)
             cGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCg                                                                                 >  1:1612213/1‑100 (MQ=255)
               aaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGaa                                                                               <  2:1396514/100‑1 (MQ=255)
                aagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAAt                                                                              >  2:1029121/1‑100 (MQ=255)
                aagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAAt                                                                              >  1:1135298/1‑100 (MQ=255)
                aagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAAt                                                                              >  2:2897574/1‑100 (MQ=255)
                 agaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATg                                                                             <  2:268829/100‑1 (MQ=255)
                  gaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGt                                                                            >  2:2427594/1‑100 (MQ=255)
                  gaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGt                                                                            <  1:2352254/100‑1 (MQ=255)
                  gaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGt                                                                            <  1:51010/100‑1 (MQ=255)
                   aaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGCGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGtt                                                                           <  2:1651082/100‑1 (MQ=255)
                   aaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGtt                                                                           >  1:2446508/1‑100 (MQ=255)
                   aaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGGGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGtt                                                                           <  1:2129072/100‑1 (MQ=255)
                    aaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTg                                                                          <  2:1942595/100‑1 (MQ=255)
                    aaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTg                                                                          >  1:937375/1‑100 (MQ=255)
                     aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTgg                                                                         <  1:811588/100‑1 (MQ=255)
                      cGCCGTATGAACCAATGCCGGAAAATGTCTAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGc                                                                        <  2:1030187/100‑1 (MQ=255)
                      cGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAGGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGc                                                                        <  1:2365341/100‑1 (MQ=255)
                      cGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGc                                                                        >  1:515315/1‑100 (MQ=255)
                      cGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGc                                                                        <  2:1654201/100‑1 (MQ=255)
                        ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCtt                                                                      <  2:2455327/100‑1 (MQ=255)
                        ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCtt                                                                      <  1:1230436/100‑1 (MQ=255)
                          gTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAg                                                                    <  1:1904662/100‑1 (MQ=255)
                               aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGccc                                                               <  1:1643895/100‑1 (MQ=255)
                                     tGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGtt                                                         >  1:2258509/1‑100 (MQ=255)
                                     tGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGtt                                                         >  2:258785/1‑100 (MQ=255)
                                      gCCGGAAAATGTCGAGGAGATCGTGGCGGAGGCAAAAGGCGCACACTAAAGTCAGAGGGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAg                                                             <  2:2755041/95‑1 (MQ=255)
                                      gCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAACGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTc                                                        <  2:879054/100‑1 (MQ=255)
                                       ccGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTcc                                                       >  1:29895/1‑100 (MQ=255)
                                        cGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTccc                                                      <  1:2271900/100‑1 (MQ=255)
                                        cGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTccc                                                      <  2:379574/100‑1 (MQ=255)
                                           aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctct                                                   <  2:1319782/100‑1 (MQ=255)
                                            aaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCtctc                                                  <  2:29895/100‑1 (MQ=255)
                                              aTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAc                                                <  2:1087503/100‑1 (MQ=255)
                                               tGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAcc                                               >  2:1711749/1‑100 (MQ=255)
                                                gTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAATTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                              <  2:2917894/100‑1 (MQ=255)
                                                 tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGTGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCt                                             <  2:1600487/100‑1 (MQ=255)
                                                 tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCt                                             <  1:2116702/100‑1 (MQ=255)
                                                  cGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTa                                            >  1:1200800/1‑100 (MQ=255)
                                                  cGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTa                                            <  1:2661527/100‑1 (MQ=255)
                                                     ggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTACCCTCTCACCCTAAcc                                         >  2:1743140/1‑100 (MQ=255)
                                                        gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCctc                                      >  1:1968225/1‑100 (MQ=255)
                                                          tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCctct                                     >  2:1525840/1‑100 (MQ=255)
                                                          tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCtctc                                    >  2:853378/1‑100 (MQ=255)
                                                           cGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTcc                                   >  2:73294/1‑100 (MQ=255)
                                                                cgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGt                              >  1:1266530/1‑100 (MQ=255)
                                                                 gcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTg                             <  2:898591/100‑1 (MQ=255)
                                                                 gcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTg                             <  1:1490817/100‑1 (MQ=255)
                                                                 gcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTg                             >  2:1062010/1‑100 (MQ=255)
                                                                  cAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                            <  2:2697456/100‑1 (MQ=255)
                                                                  cAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                            <  2:747768/100‑1 (MQ=255)
                                                                   aGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTggg                           <  1:1613388/100‑1 (MQ=255)
                                                                   aGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTggg                           <  1:2929597/100‑1 (MQ=255)
                                                                   aGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTagg                           <  1:1125242/100‑1 (MQ=255)
                                                                    ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgggg                          >  2:1700957/1‑100 (MQ=255)
                                                                    ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgggg                          <  1:305551/100‑1 (MQ=255)
                                                                      cAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                        >  2:107528/1‑100 (MQ=255)
                                                                      cAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGCTGGGGcn                        <  1:1577231/100‑2 (MQ=255)
                                                                       aaaaGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGa                       <  2:1047757/100‑1 (MQ=255)
                                                                        aaaGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAg                      >  2:2305771/1‑100 (MQ=255)
                                                                         aaGGCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                     <  1:1743151/100‑1 (MQ=255)
                                                                           ggCGCACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgggg                   <  1:2897590/100‑1 (MQ=255)
                                                                             cgcACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                 >  2:1539840/1‑100 (MQ=255)
                                                                             cgcACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                 <  2:1612204/100‑1 (MQ=255)
                                                                              gcACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTACCCCTCTCCCCGGTGGGGCGAGGGGACt                <  2:134897/100‑1 (MQ=255)
                                                                              gcACACTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACt                <  2:515312/100‑1 (MQ=255)
                                                                               cacaCTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTg               >  2:2852313/1‑100 (MQ=255)
                                                                               cacaCTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTg               <  1:2656385/100‑1 (MQ=255)
                                                                                 cacTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa              >  2:399232/1‑100 (MQ=255)
                                                                                 cacTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc             >  2:2965498/1‑100 (MQ=255)
                                                                                 cacTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTAAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc             <  2:2097961/100‑1 (MQ=255)
                                                                                  acTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAcc            <  1:1029126/100‑1 (MQ=255)
                                                                                   ttccgatctgagTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCg           <  2:2916575/91‑1 (MQ=37)
                                                                                   cTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCg           >  1:214335/1‑100 (MQ=255)
                                                                                   cTAAAGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCg           >  1:2698102/1‑100 (MQ=255)
                                                                                       aGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgcn       <  2:2600380/100‑2 (MQ=255)
                                                                                       aGTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgcg       <  1:1539849/100‑1 (MQ=255)
                                                                                        gTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgcgc      >  2:444042/1‑100 (MQ=255)
                                                                                        gTCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgcgc      >  1:2661562/1‑100 (MQ=255)
                                                                                         tCAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGcgcg     <  1:2305785/100‑1 (MQ=255)
                                                                                          cAGAGTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGt    >  2:1720991/1‑100 (MQ=255)
                                                                                           agagTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt   >  1:1307194/1‑100 (MQ=255)
                                                                                            gagTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGTCGAGGGGACTGACCGAGCGCGTTg  <  1:2139554/100‑1 (MQ=255)
                                                                                            gagTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgatcgaag  >  1:2916591/1‑93 (MQ=38)
                                                                                            gagTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTg  <  1:2170032/100‑1 (MQ=255)
                                                                                            gagTGAGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTg  <  2:501773/100‑1 (MQ=255)
                                                                                                  |                                                                                              
CCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTG  >  NC_000913/631913‑632105

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