Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,216,435 T→A 7.2% E19D (GAA→GAT yqjH ← putative siderophore interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,216,4350TA 7.2% 276.3 / 10.6 83E19D (GAA→GATyqjHputative siderophore interacting protein
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base A (3/3);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

ACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAA  >  NC_000913/3216295‑3216573
                                                                                                                                            |                                                                                                                                          
acaAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGaa                                                                                                                                    >  3:367557/1‑149 (MQ=255)
acaAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGaa                                                                                                                                    <  1:533198/149‑1 (MQ=255)
            tCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCAtt                                                                                                                        <  1:138259/149‑1 (MQ=255)
            tCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCAtt                                                                                                                        <  3:250159/149‑1 (MQ=255)
            tCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCAtt                                                                                                                        <  2:118335/149‑1 (MQ=255)
            tCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCAtt                                                                                                                        <  1:91872/149‑1 (MQ=255)
            tCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCAtt                                                                                                                        <  4:180375/149‑1 (MQ=255)
                 ttGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGaa                                                                                                                   >  3:435317/1‑149 (MQ=255)
                  tGAGGAAAGAAGAGTTTGTTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAAc                                                                                                                  <  1:101209/149‑1 (MQ=255)
                  tGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAAc                                                                                                                  <  3:493752/149‑1 (MQ=255)
                      gAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACgcgc                                                                                                              >  4:508370/1‑149 (MQ=255)
                      gAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACgcgc                                                                                                              >  2:264345/1‑149 (MQ=255)
                         agaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGc                                                                                                           <  3:267101/149‑1 (MQ=255)
                          gacgagTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCg                                                                                                          >  4:33100/4‑149 (MQ=255)
                           aagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCgg                                                                                                         >  1:413082/1‑149 (MQ=255)
                             gagTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCTCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGt                                                                                                       <  4:243898/149‑1 (MQ=255)
                                tttGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGc                                                                                                    >  1:225008/1‑149 (MQ=255)
                                   gCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCggg                                                                                                 <  2:314899/149‑1 (MQ=255)
                                       tgATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGgt                                                                                             >  4:195341/1‑149 (MQ=255)
                                        gATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGgtg                                                                                            >  3:176091/1‑149 (MQ=255)
                                         aTCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGgtgt                                                                                           >  1:452691/1‑149 (MQ=255)
                                         aTCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGgtgt                                                                                           >  2:511850/1‑149 (MQ=255)
                                            gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTAt                                                                                        >  2:150385/1‑149 (MQ=255)
                                            gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTAt                                                                                        <  3:33100/149‑1 (MQ=255)
                                            gTCAAAGCCACGCGAGGTAAAACCGTCCAGCGCCTCGCCGCCGACGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGGAGCGGGGGGTGTTAt                                                                                        >  2:247425/1‑149 (MQ=255)
                                             tCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTAtt                                                                                       <  3:316121/149‑1 (MQ=255)
                                               aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCa                                                                                     <  2:339490/149‑1 (MQ=255)
                                                aaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCAt                                                                                    <  1:547748/149‑1 (MQ=255)
                                                 aGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCAtt                                                                                   <  1:404908/149‑1 (MQ=255)
                                                     aCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTa                                                                               >  3:220593/1‑149 (MQ=255)
                                                             tAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAATACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGATCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                       >  3:333305/1‑149 (MQ=255)
                                                             tAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                       >  2:348297/1‑149 (MQ=255)
                                                             tAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                       >  4:208024/1‑149 (MQ=255)
                                                             tAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGATCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                       >  3:244838/1‑149 (MQ=255)
                                                             tAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGATCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                       >  1:446884/1‑149 (MQ=255)
                                                               aaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGt                                                                     >  2:220068/1‑149 (MQ=255)
                                                                aaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTg                                                                    >  2:265057/1‑149 (MQ=255)
                                                                aaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTg                                                                    >  2:527703/1‑149 (MQ=255)
                                                                aaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTg                                                                    >  1:144803/1‑149 (MQ=255)
                                                                aaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTACGTg                                                                    >  1:94576/1‑149 (MQ=255)
                                                                aaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGGCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTg                                                                    >  3:276318/1‑149 (MQ=255)
                                                                 aCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGa                                                                   <  2:452691/149‑1 (MQ=255)
                                                                  ccGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGAt                                                                  <  3:508370/149‑1 (MQ=255)
                                                                   cGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATg                                                                 <  2:50242/149‑1 (MQ=255)
                                                                   cGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATg                                                                 <  1:264345/149‑1 (MQ=255)
                                                                    gTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATgg                                                                <  3:417309/149‑1 (MQ=255)
                                                                     tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGt                                                               <  3:406370/149‑1 (MQ=255)
                                                                        aGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAAt                                                            <  1:400840/149‑1 (MQ=255)
                                                                           gcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAg                                                         >  1:361218/1‑149 (MQ=255)
                                                                             cTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGat                                                       <  4:70498/149‑1 (MQ=255)
                                                                                      gAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaaa                                              >  3:55795/1‑149 (MQ=255)
                                                                                      gAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaaa                                              <  1:58750/149‑1 (MQ=255)
                                                                                       aGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATc                                                           <  2:529982/135‑1 (MQ=255)
                                                                                       aGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATc                                                           >  1:529982/1‑135 (MQ=255)
                                                                                            actGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                       >  2:118810/3‑118 (MQ=255)
                                                                                            aaTGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                       <  1:118810/118‑1 (MQ=255)
                                                                                                cgcTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGc                                    <  2:361218/149‑1 (MQ=255)
                                                                                                cgcTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGc                                    <  2:144803/149‑1 (MQ=255)
                                                                                                 gcTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGgtg                                                                                            >  1:191533/1‑92 (MQ=255)
                                                                                                 gcTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGgtg                                                                                            <  2:191533/92‑1 (MQ=255)
                                                                                                 gcTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCa                                   >  2:540079/1‑149 (MQ=255)
                                                                                                 gcTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCa                                   >  4:458911/1‑149 (MQ=255)
                                                                                                 gcTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCa                                   >  4:396322/1‑149 (MQ=255)
                                                                                                   tGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCaaa                                 >  2:9604/1‑149 (MQ=255)
                                                                                                       aaCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgat                             >  4:222040/1‑149 (MQ=255)
                                                                                                          cGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataat                          >  1:215455/1‑149 (MQ=255)
                                                                                                                     cgcTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc               <  1:511850/149‑1 (MQ=255)
                                                                                                                      gcTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTAATCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCa              <  2:94576/149‑1 (MQ=255)
                                                                                                                        tCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATg            >  3:401625/1‑149 (MQ=255)
                                                                                                                        tCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATg            >  2:517445/1‑149 (MQ=255)
                                                                                                                         cAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGa           <  1:9604/149‑1 (MQ=255)
                                                                                                                         cAACGCGTAACACAGGCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGa           <  1:220068/149‑1 (MQ=255)
                                                                                                                            cgcgTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATaa        <  1:439841/149‑1 (MQ=255)
                                                                                                                             gcgTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATaaa       >  1:227423/1‑149 (MQ=255)
                                                                                                                              cgTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAAt      <  3:55726/149‑1 (MQ=255)
                                                                                                                              cgTAACACAGTCAGATCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAAt      <  4:244838/149‑1 (MQ=255)
                                                                                                                              cgTAACACAGTCAGATCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAAt      <  4:333305/149‑1 (MQ=255)
                                                                                                                              cgTAACACAGTCAGATCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAAt      <  2:446884/149‑1 (MQ=255)
                                                                                                                               gTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATg     >  3:451675/1‑149 (MQ=255)
                                                                                                                               gTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATg     >  4:424981/1‑149 (MQ=255)
                                                                                                                               gTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATg     >  1:205998/1‑149 (MQ=255)
                                                                                                                                  acacAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCaa  <  2:102216/149‑1 (MQ=255)
                                                                                                                                  acacAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCaa  <  2:293532/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
ACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAA  >  NC_000913/3216295‑3216573

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