Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,704,334 G→A 5.7% A143V (GCG→GTG)  dppC ← dipeptide/heme ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,704,3340GA 5.7% 219.1 / 2.7 70A143V (GCG→GTG) dppCdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base A (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAACCAGCAGCGACAGGCGCGCACCGTACATC  >  NC_000913/3704186‑3704480
                                                                                                                                                    |                                                                                                                                                  
gTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCCGCGCCAGCAGCAGACTTGGCAgcg                                                                                                                                                    >  5:53247/1‑149 (MQ=255)
      gcgGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCa                                                                                                                                              <  6:309462/149‑1 (MQ=255)
       cgGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCAt                                                                                                                                             >  7:409116/1‑149 (MQ=255)
         gTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGa                                                                                                                                           >  7:438125/1‑149 (MQ=255)
            aaCTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGata                                                                                                                                         <  2:109894/148‑1 (MQ=255)
             aCTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATc                                                                                                                                       <  6:171812/149‑1 (MQ=255)
               ttCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCCGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCACCAGCATGATATCGa                                                                                                                                     >  2:191805/1‑149 (MQ=255)
               ttCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGa                                                                                                                                     >  5:340749/1‑149 (MQ=255)
               ttCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGa                                                                                                                                     >  4:213552/1‑149 (MQ=255)
               ttCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCACCAGCATGATATCGa                                                                                                                                     >  2:185899/1‑149 (MQ=255)
                 ccaccaGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGAc                                                                                                                                    <  7:106042/148‑1 (MQ=255)
                  caccaGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCa                                                                                                                                  >  6:196039/1‑149 (MQ=255)
                     caGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCAcgc                                                                                                                               >  4:306606/1‑149 (MQ=255)
                       gCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCg                                                                                                                                      >  4:319847/1‑140 (MQ=255)
                       gCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCg                                                                                                                                      <  3:319847/140‑1 (MQ=255)
                          cggcggCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATa                                                                                                                          <  5:246175/149‑1 (MQ=255)
                               gcgcGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAatgat                                                                                                                     <  7:358426/149‑1 (MQ=255)
                                        aaGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGa                                                                                                             <  6:166502/148‑1 (MQ=255)
                                         aGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGAc                                                                                                            <  3:170762/148‑1 (MQ=255)
                                          gcgcACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGAc                                                                                                            >  7:398752/1‑147 (MQ=255)
                                          gcgcACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGAAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGAc                                                                                                            >  8:132438/1‑147 (MQ=255)
                                             cacaTAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGgc                                                                                                       >  5:254731/1‑149 (MQ=255)
                                                  aGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAatgatg                                                                                                                    <  3:205487/131‑1 (MQ=255)
                                                  aGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAatgatg                                                                                                                    >  4:205487/1‑131 (MQ=255)
                                                    tGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCaaa                                                                                                <  8:13225/149‑1 (MQ=255)
                                                        gcaaggcaaCGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTaa                                                                                            <  5:214463/149‑1 (MQ=255)
                                                         caaggcaaCGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAAc                                                                                           <  5:196039/149‑1 (MQ=255)
                                                           aggcaaCGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCg                                                                                         >  2:165481/1‑149 (MQ=255)
                                                               aaCGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATc                                                                                                               <  1:362254/123‑1 (MQ=255)
                                                               aaCGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATc                                                                                                               >  2:362254/1‑123 (MQ=255)
                                                                 cGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATc                                                                                   <  6:340749/149‑1 (MQ=255)
                                                                   aaGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAg                                                                                 <  6:254731/149‑1 (MQ=255)
                                                                   aaGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAg                                                                                 <  1:165481/149‑1 (MQ=255)
                                                                      gTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGcc                                                                              >  4:202730/1‑149 (MQ=255)
                                                                       tCAGTGCCAGCGCGGCGTTCCCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGcc                                                                              >  1:97263/1‑148 (MQ=255)
                                                                       tCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCGCCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAatgat                                                                                                                     >  6:301703/1‑109 (MQ=255)
                                                                       tCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCGCCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAatgat                                                                                                                     <  5:301703/109‑1 (MQ=255)
                                                                         aGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCga                                                                            <  6:417969/148‑1 (MQ=255)
                                                                           tGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGaa                                                                         >  1:64951/1‑149 (MQ=255)
                                                                            gCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGACCGCCAAAGTAACCGGCGATCAGGCCGAGaa                                                                         >  7:385525/1‑148 (MQ=255)
                                                                              cAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGACCGCCAAAGTAACCGGCGATCAGGCCGAGAATaa                                                                      <  8:385525/149‑1 (MQ=255)
                                                                               aGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAAc                                                                     <  6:53247/149‑1 (MQ=255)
                                                                                gcgcGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCACCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACg                                                                    <  1:191805/149‑1 (MQ=255)
                                                                                gcgcGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCACCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACg                                                                    <  1:185899/149‑1 (MQ=255)
                                                                                 cgcgGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGc                                                                   <  3:213552/149‑1 (MQ=255)
                                                                                 cgcgGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGc                                                                   <  4:12383/149‑1 (MQ=255)
                                                                                  gcagcgTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGcc                                                                  <  8:409116/146‑1 (MQ=255)
                                                                                   cggcgTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGcc                                                                  >  1:323373/1‑148 (MQ=255)
                                                                                      cgTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCAt                                                               >  6:200011/1‑148 (MQ=255)
                                                                                       gTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATaa                                                             <  8:438125/149‑1 (MQ=255)
                                                                                          aCCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAatgat                                                                                                                     <  3:21415/90‑1 (MQ=255)
                                                                                          aCCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAatgat                                                                                                                     >  4:21415/1‑90 (MQ=255)
                                                                                          aCCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAAt                                                            >  5:281260/1‑147 (MQ=255)
                                                                                           ccAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCa                                                          <  2:64951/148‑1 (MQ=255)
                                                                                           ccAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACATGGCAGCGCCAGCATGCTATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAg                                                         <  6:259931/149‑1 (MQ=255)
                                                                                            cAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGc                                                                   <  6:277209/138‑1 (MQ=255)
                                                                                                  aCGGGCCGAAAATTGCCACCAGCACCAGCGCCAGAAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAac                                                  <  3:306606/149‑1 (MQ=255)
                                                                                                      gCCGAAAATTGCCACCAGCACCAGAGCCAGCAGCAGACTTGGCAGCGCCAGAATGATATCGACCACGCGCATAATGATGTTATCGAGCAGGCCGCCAAAGTACCCGGCGATCAGGCCGAGAATAACGCCCATAATAAGCGATAccacta                                              >  4:1496/1‑147 (MQ=255)
                                                                                                          aaaaTTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTAc                                          <  2:323373/149‑1 (MQ=255)
                                                                                                                   accagcaccagcGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGc                                  >  8:57529/1‑148 (MQ=255)
                                                                                                                      agcaccagcGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAAc                              <  6:118817/149‑1 (MQ=255)
                                                                                                                          ccagcGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTAcc                                         <  6:217939/134‑1 (MQ=255)
                                                                                                                          ccagcGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTAcc                                         >  5:217939/1‑134 (MQ=255)
                                                                                                                           cagcGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAACcag                           >  7:326259/1‑147 (MQ=255)
                                                                                                                           cagcGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAACcag                           >  4:381721/1‑147 (MQ=255)
                                                                                                                                    cagcagACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAACCAGCAGCGACAGgc                >  7:39302/1‑149 (MQ=255)
                                                                                                                                      gcagACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGTCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGACAACCAGCAGCGACAGgcgc              <  3:202730/149‑1 (MQ=255)
                                                                                                                                            ttGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAACCAGCAGCGACAGGCGCGCACCg        >  6:159636/1‑149 (MQ=255)
                                                                                                                                               gCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAACCAGCAGCGACAGGCGCGCACCGTAc     <  6:281260/149‑1 (MQ=255)
                                                                                                                                                cAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAACCAGCAGCGACAGGCGCGCACCGTACa    <  5:159636/149‑1 (MQ=255)
                                                                                                                                                 aGCGCCAGCATGATATCGACCACGCGCAAAATGAAGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATAAGCGATAACCCAACTACCAGACAGCCAACCAGCAGCGACAGGCGCGCACCGTACAt   >  5:185966/1‑149 (MQ=255)
                                                                                                                                                  gcgcCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAACCAGCAGCGACAGGCGCGCACCGTACAt   >  3:167626/1‑148 (MQ=255)
                                                                                                                                                  gcgcCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAACCAGCAGCGACAGGCGCGCACCGTACATc  <  2:97263/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACGGGCCGAAAATTGCCACCAGCACCAGCGCCAGCAGCAGACTTGGCAGCGCCAGCATGATATCGACCACGCGCATAATGATGTTATCGACCAGGCCGCCAAAGTAACCGGCGATCAGGCCGAGAATAACGCCCATAATCAGCGATAACACAACTACCAGACAGCCAACCAGCAGCGACAGGCGCGCACCGTACATC  >  NC_000913/3704186‑3704480

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