Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,513,778 T→G 14.3% V321G (GTG→GGG)  ydcT → putative ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,513,7780TG14.3% 201.9 / 3.0 78V321G (GTG→GGG) ydcTputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (22/41);  new base G (11/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGA  >  NC_000913/1513542‑1514012
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            
tGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGa                                                                                                                                                                                                                                            <  1:681292/239‑1 (MQ=255)
        tGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTcc                                                                                                                                                                                                                                    <  1:61927/239‑1 (MQ=255)
           cAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTg                                                                                                                                                                                                                                  <  1:789453/238‑1 (MQ=255)
            agagAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACCGGTTGTGAGTCAGGGCAATAAAACAGGCGGAGAACTAGCCGTCAAGCCCACGCCCGGCCAAACGGCGAGGGTTTCCTg                                                                                                                                                                                                                                  >  2:515059/1‑237 (MQ=255)
              agaAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGGGAAAAACTGGTTGTGAGGCAGGCCAATATGAAAGGCGAAGAACTGCCCGCCACGGTCACGCCCGGAACACAGGTGGTGGGTTTCTGGTcg                                                                                                                                                                                                                              >  1:286126/1‑239 (MQ=255)
                 aaCTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTcgcg                                                                                                                                                                                                                            <  1:730019/238‑1 (MQ=255)
                      tgtgGCATGACGGGAAGCTTCGCCCTGCGACCGGACCATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCgtgatgt                                                                                                                                                                                                                      <  1:780810/239‑1 (MQ=255)
                        tgGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGGGGCAGGGCAATATGACAGGCGGAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGGGGTGGTTTTCTGGTTGCggggtgtga                                                                                                                                                                                                                    >  1:659623/1‑239 (MQ=255)
                                cGGGAAGCTTCGACAGAGAACCGGTACATATCGGCCGCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCCGTAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGTCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACACCAGGTGATGGTTTCCTGCTCGCGTGATGTGATGGTGCCg                                                                                                                                                                                                           <  2:1203665/239‑1 (MQ=255)
                                agggAAGCTGCGCCCTGCGACCGGAACATATCCGCGTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCTGTGCAATATCAGGGCGAGGAAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCg                                                                                                                                                                                                           <  2:717505/238‑1 (MQ=255)
                                          cGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAAATCGTTTTGAACTGAAATTGGAAGGCGGTGAAAAAATGCCTGTGAGTGAAGGCAATATTAAAGGGGGAGAAATGACCGCCCCGCCCAAGCCCCGACCACAGGGGAGGGGTTCCCGGTAGCGGGCGGGGATGGGGCCggcgggggggg                                                                                                                                                                                                 >  2:1049208/1‑229 (MQ=255)
                                           gCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTgagga                                                                                                                                                                                                <  1:118508/239‑1 (MQ=255)
                                           gCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTgagga                                                                                                                                                                                                <  1:564168/239‑1 (MQ=255)
                                            cccTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTgaggag                                                                                                                                                                                               <  2:1014104/239‑1 (MQ=255)
                                                 cGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAAATGCTTGTGAGTCAGGCCAATATGACAGGCGGAGAACGGCCTGCCCCGCTCAAGCCCGGACAAAAGGTGAGGGGTTCCTGGGCGCGTGCGGTGATGGTGCCCGCTGGtgtggggggggg                                                                                                                                                                                           >  2:885081/1‑228 (MQ=255)
                                                                 gCCTCAACACCCCTGGTGAACTGCCGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAGCGGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTa                                                                                                                                                                          <  1:584408/239‑1 (MQ=255)
                                                                 gCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTa                                                                                                                                                                          <  1:107661/239‑1 (MQ=255)
                                                                 gCCCCAACACCCCGGGTGAGCTGCAGGCCAATGGCACGGGCCAGGCGGTGCAAGATCAGGGCGCGGGAACTCGGTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGGGATGGTTTCCTGGGCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTa                                                                                                                                                                          <  2:910080/239‑1 (MQ=255)
                                                                           cccTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGt                                                                                                                                                                                                                <  1:390034/192‑1 (MQ=255)
                                                                             cTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACCGCTTGTGAGTCAGGCCAATATGAAAGGCGAAGAAATGCCCGCCCCGCTCACGCCCGGACAAAAGGGGATGGGTTTCTGGTCGCGGGATGTGGTGGGGCCGGCGGCTGAGGAGAGGGGAATGGCGGTGACTGTATTGCAATCaacg                                                                                                                                                             >  2:24530/1‑236 (MQ=255)
                                                                                    aCTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTc                                                                                                                                                      <  1:1147986/239‑1 (MQ=255)
                                                                                         aGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAgg                                                                                                                                                  <  2:397207/238‑1 (MQ=255)
                                                                                            ccAATGGCACGATGCAGTGGGTGGAATATCAGTGCGCAGCATGTCGTTTTGAACTGAAATTGAAGGGCGGTGAAAAACTGCTTGTGAGTCAGGGCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCGTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGt                                                                                                                                                 <  2:538988/236‑1 (MQ=255)
                                                                                            ccAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGt                                                                                                                                                 <  1:885081/236‑1 (MQ=255)
                                                                                                tGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTa                                                                                                                                          >  1:238397/1‑239 (MQ=255)
                                                                                                     cGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCCTGTGAGTCAGGCCAATATGACAAGGGAAGGACTGCCTGCCACGCTCACGCCCGGACAACAGGGGAGGGTTTTCTGGTGGCGTGGTGTGGTGGCGCCGCCGGGTGGGGAGAGGGGAATGGCGGAGGATGGATTGCGATCTGCTTCACGGGCAGGGGTGGGGAAgggt                                                                                                                                     >  2:609905/1‑237 (MQ=255)
                                                                                                       aTCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGGCAATATGACAGGGGAAGAAATGGCTGCCACGCTCACGCCCGGACAAAAGGGGAAGGTTTCCTGGGCGCGGGAGGGGGTGGTGCCGGTGGGTGAGGGGAGGGGGAAGGCGGTGGAAGTATTTCAATCACCTTTCCGGCCAGGGCGGGGGAAgggtg                                                                                                                                    >  2:870819/1‑235 (MQ=255)
                                                                                                       aTCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAgg                                                                                                                                                  <  1:796295/224‑1 (MQ=255)
                                                                                                        tCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACCGCGTGTGAGTCAGGCCAATATGACAGGGGGAGGACTGCCTGCCACCCTCACCCCCGGGCAAAAAGGGACGGTTTTCTGGGCCCGTGATGTGATGGGGCCGCCGGGTGAGGAGGGGGGGAAGGCGAGGGATGGATTGGAAGTCACCTCCCGGTCCGGCTCGGGGGAAgtggcc                                                                                                                                   >  2:189177/1‑234 (MQ=255)
                                                                                                              cGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGt                                                                                                                               <  1:567809/236‑1 (MQ=255)
                                                                                                              cGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGACCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGt                                                                                                                               <  1:652412/236‑1 (MQ=255)
                                                                                                                 tGCCATATCAGGGCGCGGCAACTCGTTTTGAAATGAAATTGAACGGCGGTGAAAAAATGGCTGTGATGCAGGCCAATATGACAGGCGAAGAAATGCCTGCCACGCTCACGCCCGGCCAACAGGTGATGGGTTTCTGGGCGCGGGAAGTGCTGGGGCCGCCGGGTGAGGGGAGGTGGATGGCGGTGGATGTTTTGTCATCACCTTCACGTCCAGGGCTGGGGAAGGGGTCCGGGttct                                                                                                                           >  2:738475/1‑237 (MQ=255)
                                                                                                                 tGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCCGGCCAATAAGACAGGCGGAGAACTGCCTGCCCCGCCCACGCCCCGCCCACAGGGGAGGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCCTGTTGTGGAGAGGTGTATGGCGCTGGATGGATTGTCATCACCTTCACGGACCGGGCTGGGGAAGGTGGCCGGtttac                                                                                                                           >  2:1000741/1‑235 (MQ=255)
                                                                                                                 tGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTc                                                                                                                                                           >  1:842616/1‑205 (MQ=255)
                                                                                                                 tGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTc                                                                                                                                                           <  2:842616/205‑1 (MQ=255)
                                                                                                                 tGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGGGATGGTTTCCTGGTCGCGTGGTGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCCCGTCCAGGTCTGGGTAAGGGGTCCGGGttct                                                                                                                           >  1:288520/1‑237 (MQ=255)
                                                                                                                 tGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGGGAAAAACTGCCTGGGAGGCAGGCCAATATAACAAGCGGAGAACTGCCCGCCACGCTCACGCCCGGGCCAAAGGGGATGGGTTCCTGGGCGCCTGATGTGATGGTGCCGCTGGTTGAGGAGAGGGGGATGGCGATGGATGGATTGCAATCCACTTCACGTCCCGGGCTGGGGGAGGGGGCCGGTttcttt                                                                                                                         >  2:189434/1‑238 (MQ=255)
                                                                                                                  gCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTttcttc                                                                                                                         <  1:738475/238‑1 (MQ=255)
                                                                                                                     atagCAGGGCGCGGCAACTCGGTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAGGTGTCCGGTTTCTTCTGGCa                                                                                                                    <  1:5144/239‑1 (MQ=255)
                                                                                                                      tatCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCa                                                                                                                    <  1:401502/239‑1 (MQ=255)
                                                                                                                           gggCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGACGAATTGCCTGCCACGCTCACGCACGGCCAACTGGCGATGGTTTCCTCGTCGCGTGATGTGCTGGTGCCGCTGGTTGAGGAGAGGGGAATGGCGCTGCATGTCTTGCAATCCCCTTCACGTCCCGGTGTGGGGAAAGTGTCCGGTTTCTTCTGGcgttat                                                                                                                >  2:321323/1‑233 (MQ=255)
                                                                                                                           gggCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCCGGTCTGGGTAAGGTGTCCGGTATCTTCTGGCATAAt                                                                                                                >  1:280730/1‑238 (MQ=255)
                                                                                                                               gcgGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCgg                                                                                                            <  1:870819/238‑1 (MQ=255)
                                                                                                                                cgGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCggg                                                                                                           >  1:930509/1‑238 (MQ=255)
                                                                                                                                     aCTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCggggctggg                                                                                                     <  1:136730/239‑1 (MQ=255)
                                                                                                                                        cGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCggggctggggct                                                                                                  <  1:321323/239‑1 (MQ=255)
                                                                                                                                           tttGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGGGATGGTTTCCTGGTCGCGTGATGTGATGGGGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGGGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGtt                                                                                               >  1:1063789/1‑239 (MQ=255)
                                                                                                                                             tGAACTGAAATTGAGCGGCGGTGAAAAAGTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGtttt                                                                                             <  2:827583/239‑1 (MQ=255)
                                                                                                                                             tGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGtttt                                                                                             <  1:24530/239‑1 (MQ=255)
                                                                                                                                                           aCGGCGGTGCAAAACTGCTTGTGAGTCCGGCCCATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGCCAACCGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCCGGTTGCGGAGAGGTGACTGGCGCTGAATGTATTGCCATCCCCTTCACGTCCAGGTCTGGGGAAGGCGGCCGGTTTCTTCTGGCCTAATCCCGGGCTGGCGCCGTTTTTACTTCCGCTTgg                                                                                >  2:1208264/1‑238 (MQ=255)
                                                                                                                                                           aCGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTgg                                                                                <  1:288810/238‑1 (MQ=255)
                                                                                                                                                               cggTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGGGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGccc                                                                             >  1:1217442/1‑237 (MQ=255)
                                                                                                                                                                 gTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGACTGCCACGCTCAAGCCCGGACAACAGGGGATGGTTTCCTGGTCGCGGGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGGATTGGAATCACCTTCACGTGCAGGGCTGGGGAAGGGGGCCGGTTTCTTCTGGGATAAACCGGGGCGGGGGCGGTTTGTACTGCTGCTTGGGCCggca                                                                         >  2:1125024/1‑236 (MQ=255)
                                                                                                                                                                     aaaaCTGCTTGTGAGTCAGGCCAATATGACAGGCTAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAAtg                                                                      <  2:768895/238‑1 (MQ=255)
                                                                                                                                                                           gCTTGTGCGTCAGGCCAATATGACAGTAGTAGAACTGCCTCCAACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTTATGGTGCCGCTGACTGAGGAGAGGTGATTGGCGATGAATGTATTGCAATCTCCTTCATGTCCAGGTCTGGGTAAGGTGGCCGGTTTCTTCTGGCATAATACGGGGCTGGGGCTGTTTTTATTGCTGCTTGGCTCTCTACTTTGGTTTg                                                               >  2:449006/1‑239 (MQ=255)
                                                                                                                                                                              tgggAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGATCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCa                                                            <  1:234010/236‑1 (MQ=255)
                                                                                                                                                                                   gTCAGGCCAATATGACGGGCGAAGAACTGCCGGCCAAGCTCACGCCCGGACAGCAGGTGATGGTTTCCTGGGCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGt                                                        <  2:1041442/238‑1 (MQ=255)
                                                                                                                                                                                             tatGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGc                                              >  2:300727/1‑238 (MQ=255)
                                                                                                                                                                                                  cAGGCGAAGAACTGCCTTCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGGTGAGGAGAGGTGAATGGCGATGAGTGTATTGCAATCACCTTCTCGTCCAGGTCTGGGTAAGGTGTTCGGTTTCTTCTGGCATAATCCGGGGCTCGGGCTGTTTTTACTGCTGCTTGCCGCGCTCACGTGGTTTGGCATTGTCTTTTGTGGCTCgct                                         >  2:33211/1‑238 (MQ=255)
                                                                                                                                                                                                  cAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCggggctggggc                                                                                                   >  1:92928/1‑180 (MQ=255)
                                                                                                                                                                                                  cAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCggggctggggc                                                                                                   <  2:92928/180‑1 (MQ=255)
                                                                                                                                                                                                  cAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCgct                                         >  1:346608/1‑238 (MQ=255)
                                                                                                                                                                                                  cAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCgct                                         >  1:1203915/1‑238 (MQ=255)
                                                                                                                                                                                                    ggCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTc                                            >  1:665127/1‑233 (MQ=255)
                                                                                                                                                                                                    ggCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTc                                            <  2:665127/233‑1 (MQ=255)
                                                                                                                                                                                                       gaagaaCTGCCGGCCACTCTCACGCCCGGGCAACAGGTTATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGa                                    <  1:1059443/238‑1 (MQ=255)
                                                                                                                                                                                                        aagaaCTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCa                                                            >  1:845568/1‑213 (MQ=255)
                                                                                                                                                                                                        aagaaCTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGAGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCa                                                            <  2:845568/213‑1 (MQ=255)
                                                                                                                                                                                                                     cgcGCTCACGCCCGGACAACAGGGGATGGTTTCCTGGTCGCGTGATGTGATGTGGTCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCAGAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCGCGCTGCTGACGCTGTTATGGCAg                      <  2:1218354/236‑1 (MQ=255)
                                                                                                                                                                                                                     cACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAg                      <  2:930509/238‑1 (MQ=255)
                                                                                                                                                                                                                     cACGCTCACGCCCGGACAACAGGTGATGGTTGCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAg                      <  1:569069/238‑1 (MQ=255)
                                                                                                                                                                                                                       cGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAg                      <  2:363653/236‑1 (MQ=255)
                                                                                                                                                                                                                       cGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAg                      >  1:363653/1‑236 (MQ=255)
                                                                                                                                                                                                                         cTCACGCCCGGACAACATGTGATGGTTTCCTGGGCGCGTGATGTGATGGTGCCGCTGGTTgagga                                                                                                                                                                                                <  2:811259/65‑1 (MQ=255)
                                                                                                                                                                                                                                  ggACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTATACttt        >  1:267983/1‑239 (MQ=255)
                                                                                                                                                                                                                                   gACAACAGGTGATGGTTTCCTGGTCGCGTGGTGTGATGGGGCCGATGGTTGAGGAGAGGTGAATg                                                                                                                                                                                      >  1:811259/1‑65 (MQ=255)
                                                                                                                                                                                                                                       acaGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTTTATttttcatg   >  2:1075435/1‑235 (MQ=255)
                                                                                                                                                                                                                                        caGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGGTTTTATACTTTTTACGa  >  1:1188810/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            
TGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACATATCCGCCTCAACACCCCTGGTGAACTGCAGGCCAATGGCACGATCCAGGCGGTGCAATATCAGGGCGCGGCAACTCGTTTTGAACTGAAATTGAACGGCGGTGAAAAACTGCTTGTGAGTCAGGCCAATATGACAGGCGAAGAACTGCCTGCCACGCTCACGCCCGGACAACAGGTGATGGTTTCCTGGTCGCGTGATGTGATGGTGCCGCTGGTTGAGGAGAGGTGAATGGCGATGAATGTATTGCAATCACCTTCACGTCCAGGTCTGGGTAAGGTGTCCGGTTTCTTCTGGCATAATCCGGGGCTGGGGCTGTTTTTACTGCTGCTTGGCCCGCTAATGTGGTTTGGCATTGTCTATTTCGGCTCGCTGCTGACACTGTTATGGCAGGGATTTTATACTTTTGACGA  >  NC_000913/1513542‑1514012

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