Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 14,704 A→G V179V (GTA→GTG dnaJ → chaperone protein DnaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091314,7040AG98.5% 215.2 / NA 68V179V (GTA→GTGdnaJchaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (28/38);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCG  >  NC_000913/14461‑14940
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
cGGCGATATTTGTGGGGGCGGACGTGGTCGTCACCGTGCGGCGCGCGGTGCTGATTTCCGCTAGAACATGTAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGGCCAAAGAGAGCCGCATTCCGACTCTGGAAGAGTGTGACGGTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGGCATGGTTCGGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTTcagcag                                                                                                                                                                                                                                        <  2:445905/250‑1 (MQ=255)
   cGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCt                                                                                                                                                                                                                                    <  1:353713/251‑1 (MQ=255)
            tGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCacac                                                                                                                                                                                                                            <  1:89264/250‑1 (MQ=255)
                  gggACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGTATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAgg                                                                                                                                                                                                                     <  1:206518/250‑1 (MQ=255)
                     gcGTGGGCGGCAACGGGGGGCGCGTGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGTCCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACTCTGTCAGGGcc                                                                                                                                                                                                                  <  2:477125/250‑1 (MQ=255)
                      cGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCcg                                                                                                                                                                                                                 <  1:277735/251‑1 (MQ=255)
                       gTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCcgc                                                                                                                                                                                                                >  1:43044/1‑251 (MQ=255)
                            cgtcAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTa                                                                                                                                                                                                            <  1:379560/250‑1 (MQ=255)
                                aaCGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTg                                                                                                                                                                                                       >  1:137006/1‑251 (MQ=255)
                                        gcgcgcGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGa                                                                                                                                                                                               <  1:119651/251‑1 (MQ=255)
                                           cgcgGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGCAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACACTGATCAACGATcc                                                                                                                                                                                            >  2:229076/1‑251 (MQ=255)
                                             cgGTGCTGATTTACGCTATAAAATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCTTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCGGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGt                                                                                                                                                                                          <  2:279058/251‑1 (MQ=255)
                                                 gCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGca                                                                                                                                                                                       <  1:437208/250‑1 (MQ=255)
                                                      tttACACTATAACCTGGAGCTAACCCTCGAAATCGCTGTCCGTGGCGCGACCACAGAGATCCGCATTCCGACTCTGGGAGAGTGTGACGTTTGCCACGGCAGCGGTGCACAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGACGTGCAGCCAGCCTGTCAACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAAt                                                                                                                                                                                 >  2:668427/1‑251 (MQ=255)
                                                         aCGCTATAACATGGAGCTTACCCTCGAAGAAGCTGTACGTGGCGGTACCAAAGAGATCCGCATGCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgt                                                                                                                                                                               <  2:678971/250‑1 (MQ=255)
                                                          cGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtca                                                                                                                                                                             <  1:439230/251‑1 (MQ=255)
                                                           gCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCTGTGCAACAAATgtcat                                                                                                                                                                            <  1:514616/251‑1 (MQ=255)
                                                            cattaACATGGAGCTCATCCTCGCATAAGCGTTACGTGGCGTGAACAAATAGATCCGCATTCCGACCCTGTAAGTGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCTACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCAGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGGGCAACAAATgtcatg                                                                                                                                                                           <  2:150789/248‑1 (MQ=255)
                                                                   aTGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGtcatggtcatgg                                                                                                                                                                    >  1:629722/1‑251 (MQ=255)
                                                                     ggAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGAGCCGCATTCCGACTCTGGAAGAGTGTGACGGTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTGTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTcatggtcatggt                                                                                                                                                                   <  1:93515/250‑1 (MQ=255)
                                                                          tCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGTTCACAGAt                                                                                                                                                                                              >  2:188067/1‑218 (MQ=255)
                                                                          tCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGAt                                                                                                                                                                                              <  1:188067/218‑1 (MQ=255)
                                                                          ccACCCTCGAATAAGCTGTACGTGGCGTGACCAAAGAGAGCCGCATTCCGCATCTGGAAGAGTGTGACGTTTGCCACGGTAGCGTTGCAAAACCAGGTACACAGCCGCAGACCGGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCGGGGCCGCGGTACGCTGATCAAAGATCCGAGCAACAAATGGCATGGTCATGGTCGTGtt                                                                                                                                                             <  2:22600/250‑1 (MQ=255)
                                                                                 cGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAgcg                                                                                                                                                        <  1:594888/249‑1 (MQ=255)
                                                                                 cGAAGAAGCTGGACGGGGCGTGACCAAAGAGATCCGCATTCCGACTCTTGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGAACTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTGGAGCgca                                                                                                                                                      <  2:286495/251‑1 (MQ=255)
                                                                                    agaagCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCgcagc                                                                                                                                                    >  1:94133/1‑250 (MQ=255)
                                                                                    agaagCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGAACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGAgcagg                                                                                                                                                    >  2:474798/1‑249 (MQ=255)
                                                                                             aCGTGGCGTGACCAAAGAGATCCGTATTCCGACTCTTGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTTTCAGGGCCGCGGTACGATGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCATCAAAACGCTGt                                                                                                                                          <  2:393207/251‑1 (MQ=255)
                                                                                                 ggcggtaCCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTATCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCGGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGt                                                                                                                                      <  2:648939/245‑1 (MQ=255)
                                                                                                      gACCAAAGAGTTCCGCATTCCGACTCTGGAAGGGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTTGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTaaa                                                                                                                                  <  2:14059/250‑1 (MQ=255)
                                                                                                        ccAAAGAGATCCTCATTACGATTCTGGCAGAGGTGTACGTTTGCCTCGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGGCATGGGTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGGCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATc                                                                                                                               <  2:449685/251‑1 (MQ=255)
                                                                                                               gaTCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGtt                                                                                                                                                             >  2:594568/1‑214 (MQ=255)
                                                                                                               gaTCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGtt                                                                                                                                                             <  1:594568/214‑1 (MQ=255)
                                                                                                                   cGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAggggg                                                                                                                    >  1:36208/1‑250 (MQ=255)
                                                                                                                       ttCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACCAATGTCATGGTCATGGTCGTGTTGAGCGAAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGac                                                                                                                >  1:434790/1‑251 (MQ=255)
                                                                                                                       ttCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCTTGTTGAGCGCAGCAAAACGCGGTCCGTTAAAATCCCGGCAGGGGTGGac                                                                                                                >  1:434745/1‑251 (MQ=255)
                                                                                                                         ccGACTCTGGAAGAGTGTGACGTTTACCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTCAAATCCCGGCAGGGGTGGacac                                                                                                              >  1:107125/1‑251 (MQ=255)
                                                                                                                            aCTCTGGAAGAGGGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTAAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGACGCGGTACGCTGATCAAAGATCCGTGCAACAACTGTCATGGGCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGGGGGCACTgg                                                                                                           >  2:494426/1‑251 (MQ=255)
                                                                                                                                    agagTGTGACGTTTGCCACGGAAGCGGAGCAAGACCAGGTACACAGCCGCAGACTTGGCCGACCTGTAATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGAACAGCAAACCTGTCCACACTGTCAGAGCCGCGGTACGCTGATCAAAGACCCGTGCAACAAATGTCATGCTCATGGTCGTGTTGCGCGCACAAAAACGCTATCAGGTAAAATCCCGGCGGGGGTGGACCCAGGCGACCGcc                                                                                                   >  2:351149/1‑250 (MQ=255)
                                                                                                                                       gtgtGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGGCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGTACACTGGAGACTGCATcc                                                                                                <  2:94133/251‑1 (MQ=255)
                                                                                                                                         gtgACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGt                                                                                              >  1:581874/1‑251 (MQ=255)
                                                                                                                                           gACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCt                                                                                            >  1:367476/1‑251 (MQ=255)
                                                                                                                                             aCGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACCAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACCATGGGGACCGCATCCGTCtt                                                                                           >  1:532170/1‑251 (MQ=255)
                                                                                                                                              cGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGACGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTg                                                                                          >  1:549574/1‑251 (MQ=255)
                                                                                                                                                  tGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGAGGGAAACTGGAGACCGCATCCGTCTTTCggg                                                                                      >  2:642351/1‑251 (MQ=255)
                                                                                                                                                     cACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCg                                                                                    >  1:512015/1‑250 (MQ=255)
                                                                                                                                                        ggTAGCGGTGCAAAACAAGGCACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCAGCGGTACGCTGATCAACGATCCGTGCAACACATGTCATGGTCATGGTCGTGTTGAGCGCAGCTAAACGCTGTCCGTTAAAATCCCGGCAGGGGGGGACCCTGGAGAcgacatgcggattgaggtcgacgg                                                                                >  2:310865/1‑228 (MQ=255)
                                                                                                                                                            gCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAAGCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGaa                                                                            <  1:406188/251‑1 (MQ=255)
                                                                                                                                                                 gCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCggg                                                                       <  1:132665/251‑1 (MQ=255)
                                                                                                                                                                 gCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCgg                                                                        <  1:561286/250‑1 (MQ=255)
                                                                                                                                                                     aaCCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTACTGGGCAGGTGCAGATGCGTCAGGGATTCGTTGCCGTGCAGCAGACCGGTCCACACTGTCAGGGCCGCGGTACGCTGAACAAAGATCCCTGCAACAAATGTCATGGTCATGGGCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTCGACACTGGAGACCGCATCCGTCTTGCGGGc                                                                                     <  1:257726/233‑1 (MQ=255)
                                                                                                                                                                            gacacAGCCGCAGACCTGGCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGgcg                                                            <  1:429653/250‑1 (MQ=255)
                                                                                                                                                                             acacAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGAGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAACCCTGTCCACACTGTCAGAGCCGCGGTACGCTGATAAAAGATCCGTGCAACAAATGTCATGGACATGGTCATGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATACCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGGAGGAGAAGCGGGCGTGCATGgcgc                                                           >  1:630821/1‑251 (MQ=255)
                                                                                                                                                                                          cctGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGAt                                              >  1:44961/3‑251 (MQ=255)
                                                                                                                                                                                          cctGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGAt                                              <  1:370282/249‑1 (MQ=255)
                                                                                                                                                                                            tGTCCGACCTGTCCTGGCTCTGGACCGGTGCAGCTGCGTCAGGGCTTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGCTCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGCCCGCCTCCGTCTTGCGGGCGAAGGTACAGCGGGCGCGCATGGCGAACCGGCAgggcatcg                                            >  2:44897/1‑245 (MQ=255)
                                                                                                                                                                                               ccGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTa                                         <  1:306945/251‑1 (MQ=255)
                                                                                                                                                                                                cGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTAc                                        <  1:552910/251‑1 (MQ=255)
                                                                                                                                                                                                        cATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgttcaggt                                <  1:376501/251‑1 (MQ=255)
                                                                                                                                                                                                        cATGGTTCTGGTCAGGTGCAGAGGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgttcaggt                                <  1:484155/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGCCACTGGAAACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTcaggatccggttt                        >  2:323391/1‑249 (MQ=255)
                                                                                                                                                                                                                 ggTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTaa                       <  1:197023/251‑1 (MQ=255)
                                                                                                                                                                                                                       gTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCa                 <  1:637219/251‑1 (MQ=255)
                                                                                                                                                                                                                              tGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAACGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCCTCCGTCTTGCGGGCGAAGGGGAAGAGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAACACCCGAtt          >  2:163833/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gtcAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATTCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTc       <  1:288963/249‑1 (MQ=255)
                                                                                                                                                                                                                                 gtcAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTc       >  1:26773/3‑251 (MQ=255)
                                                                                                                                                                                                                                     gggATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGc   >  1:573966/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggATTCTTTGCCGTGCAGCATGCCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCTATCTGTACGTTCAGGTTTAGGTTAAACAGCACCCGATTTTCGAGCg  <  2:538526/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
CGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCG  >  NC_000913/14461‑14940

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: