Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,019,572 A→C 14.2% T138P (ACA→CCA)  ssnA → putative chlorohydrolase/aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,019,5720AC14.2% 169.2 / 2.2 65T138P (ACA→CCA) ssnAputative chlorohydrolase/aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (3/50);  new base C (8/2);  total (12/53)
Fisher's exact test for biased strand distribution p-value = 1.73e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCG  >  NC_000913/3019340‑3019809
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
cTCGCACAACCATTTTTACTCGGGGCTTTTCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGACGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTCTATCGGCGGTTGGCTATCcccacttc                                                                                                                                                                                                                                              >  2:806522/1‑235 (MQ=255)
        aCCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCa                                                                                                                                                                                                                                    <  2:1358502/239‑1 (MQ=255)
        aCCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCa                                                                                                                                                                                                                                    <  1:615674/239‑1 (MQ=255)
        aCCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACAGCGCCCCTTGCCCGGATTTCATCCCAACGCTGAAAAATCTCTGGTGGCGGCTCGAGCGCGCCCTTGATGAAGAGGCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCa                                                                                                                                                                                                                                    <  2:268893/239‑1 (MQ=255)
                 aCTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCCCTCCCCACATCGCGACGCCTTTTTtaaa                                                                                                                                                                                                                          >  1:572312/1‑239 (MQ=255)
                  cTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAg                                                                                                                                                                                                                         <  1:608539/239‑1 (MQ=255)
                    cGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCAGTCCCCACATTCCGCGAGCCCTTTTTTTAAAGt                                                                                                                                                                                                                        >  1:147565/1‑239 (MQ=255)
                     ggggCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTGGCTCTCCCCCTTCCCCGCCCCCTTTTTTAAAATTg                                                                                                                                                                                                                      >  1:480220/1‑239 (MQ=255)
                            tCCCGCGGAATTATGGCAAACCTCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTgcg                                                                                                                                                                                                               <  1:645639/239‑1 (MQ=255)
                            tCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCCCCCTTCCCCGCCCCATTTTTAAAAATTGGCCtc                                                                                                                                                                                                                 >  1:336441/1‑237 (MQ=255)
                                gcgGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCTGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGcgcgcg                                                                                                                                                                                                           <  2:643324/239‑1 (MQ=255)
                                     atcatcccAAGCATCGCCCCTTGCCCGTACTTCATCTCAACGCCGAAAAATCTCTGGTGGTGGCTCGTGCCCGCCTCTGATTAAGAGTCGCTCTATTACAGCGGACTGATTTGTACCCTGGCAGCGATTAAGCGCCGATGTTCATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTATCCCCATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGAc                                                                                                                                                                                                      <  2:1020656/232‑1 (MQ=255)
                                     aTTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGAc                                                                                                                                                                                                      <  1:1266862/239‑1 (MQ=255)
                                     aTTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGAc                                                                                                                                                                                                      <  1:860337/239‑1 (MQ=255)
                                     aTTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGAc                                                                                                                                                                                                      <  1:115925/239‑1 (MQ=255)
                                     aTTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGAc                                                                                                                                                                                                      <  2:286083/239‑1 (MQ=255)
                                                cATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCCCTCTCCCCACTTCCCGGCGCACTTTTTTAAATTGGCCCGCGCGCGGTGACCTGTGTTTGaa                                                                                                                                                                                           >  1:145530/1‑239 (MQ=255)
                                                   cGCCCCGTGCGCGGAATCCAGCGCAACGCTGAAAAAGCGCTGGTGGCGGCTCGATCGCGCACTTGAAGAAGAGTCGCTGTAATACAGCGGACTGATTTGTACTCTGGAAGCGAGTAAGAGTGGATGTACATCGGTTATAGATCGCCATGCCTCTCCGGCGTATAACGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAtcta                                                                                                                                                                                       <  2:1340629/239‑1 (MQ=255)
                                                      cccTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTg                                                                                                                                                                                    <  1:22612/239‑1 (MQ=255)
                                                                 aTTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGGGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTaacaa                                                                                                                                                                         <  1:969440/238‑1 (MQ=255)
                                                                 aTTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTaacaa                                                                                                                                                                         <  1:333550/238‑1 (MQ=255)
                                                                        ctcAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCa                                                                                                                                                                 <  2:556875/239‑1 (MQ=255)
                                                                                 gAAAAATCTCTGGTGGCGGCACGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGGAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATGTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTg                                                                                                                                                         <  1:1293771/238‑1 (MQ=255)
                                                                                     aaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:1093766/239‑1 (MQ=255)
                                                                                     aaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  2:1119850/239‑1 (MQ=255)
                                                                                     aaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:589082/239‑1 (MQ=255)
                                                                                     aaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:430480/239‑1 (MQ=255)
                                                                                     aaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:1308940/239‑1 (MQ=255)
                                                                                     aaTCTCTGGGGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  2:1107534/239‑1 (MQ=255)
                                                                                     aaTCTCTGATGGCTGCTCGATCGCGCCCTTGATGAAGAGTAGCTCTATTACAGCGGACTGATTTGTCCCCTGGAAGCGATTAAGAGCGGATGTCCATCTGTTATCGATCACCATGCCTCTGCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  2:1068358/239‑1 (MQ=255)
                                                                                     aaTCGCTGGGGGCGGCTCGAGCGCGTCCTTGACGAAGAGTCGCGCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:859148/239‑1 (MQ=255)
                                                                                              tggCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTag                                                                                                                                            <  1:977407/238‑1 (MQ=255)
                                                                                                 cggcTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTGGCTCTCCCCCCTTCTCCCCTCCTTTTTTTAAATTTGGCTGCGCGCGGGGGCCCGTTTTTTTACTATTTACCCTCAAAACGGCATCACACAGATGTGGGGAGGTGGaggaga                                                                                                                                        >  2:977407/1‑239 (MQ=255)
                                                                                                  ggcTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAgaagaa                                                                                                                                       <  1:1285163/239‑1 (MQ=255)
                                                                                                  ggcTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAgaagaa                                                                                                                                       <  2:550246/239‑1 (MQ=255)
                                                                                                   gcTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTGGGTCTCCCCCCTTGCCGCCCCCGTTTTTAAAATTTGGCTGCGCGCGCAGGCCCTTTTTTTAACCACTTACCGCAACACACGCCAAAAAGAGCTTCAGGAGGGGGTATGAGaaa                                                                                                                                      >  2:1271450/1‑239 (MQ=255)
                                                                                                                             cGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTtgattgat                                                                                                            <  1:1231936/239‑1 (MQ=255)
                                                                                                                                tctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTaacaa                                                                                                                                                                         <  1:381511/175‑1 (MQ=255)
                                                                                                                                tctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAg                                                                                                         <  2:997020/239‑1 (MQ=255)
                                                                                                                                  tATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:111187/194‑1 (MQ=255)
                                                                                                                                          cGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGGTCTCCACAATTCGCGGCCCCATTTTTTAAATTGTGCTTCGGGCGCGGTCCCGTTTTTTAAATCTTGTCCGTAACAAAAGCCTCCAAGAGGTCCAGGAAGGTGTAGAAGAAAAAACCCGTTTTGCCCCTTTGTTtgttgacggggagaca                                                                                                >  2:831626/1‑228 (MQ=255)
                                                                                                                                            gACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTGGCTCTCCTCACTTCGCGGCGCCTTTTTTTAACATGGGGTGCGGGGGGGGACCCGTTTTTTAACTACTGAACCCAAAAACGGGATCAAAAAGATGCAGGGGGGGGGAGAAAAAAAAAACCTGTTCGCCGCCTTGTTTGTTGAGGGGGAGAAAgag                                                                                             >  2:1316219/1‑237 (MQ=255)
                                                                                                                                                      ttCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGc                                                                                   <  1:597906/239‑1 (MQ=255)
                                                                                                                                                      ttCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGc                                                                                   <  1:1111813/239‑1 (MQ=255)
                                                                                                                                                            ggAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGGCTTCCCCCCTTCGCCGCGCCTTTTTTAAAAATGGCGTGCGGGCGACGACCCCTTTTGTAACTACTGAACCTAAAAAAAGCAGCAAAGAGTTGTGGGAAGGGGGAGAAGAAAAAAACCACTTTGCCCCCTTGATTGTTGGGGGGGAGAAAACGACAAACGAGCCgaa                                                                               >  2:87448/1‑235 (MQ=255)
                                                                                                                                                                  gATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCtggtgg                                                                       <  1:173111/239‑1 (MQ=255)
                                                                                                                                                                   aTTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCtggtgg                                                                       <  2:691083/238‑1 (MQ=255)
                                                                                                                                                                       agagCGGATGCAAACCGGTTATCTAGCACCATGCCTCTCCGGCGTAAATGGGCGGGTCGCACGCCACATTGCGCGACGAATTTTTTAAAAGTTGGCCTGGGCGCGAGGACCTGTTTTGAAACTCCTGACCGTAACAACGCCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGaa                                                                     <  2:147565/237‑1 (MQ=255)
                                                                                                                                                                                      cttttctCGATTATCATGTCTCTTAGTCGTTGATCGCCGGGTCCCTCTACACATTGCGCGACGCATTGTTAAAAGTTGGCCGGCGCGAGATGACCTGATTTGAAGCTACCGACCTTAACACCGGCAACAAGGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTGAGCCGGTTTGATTCATGGGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACCTATCGGTGCTCAc                                                   <  2:533883/233‑1 (MQ=255)
                                                                                                                                                                                      cGGTTATCGATCCCCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCAc                                                   <  1:1103431/239‑1 (MQ=255)
                                                                                                                                                                                      cGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCAc                                                   <  1:351249/239‑1 (MQ=255)
                                                                                                                                                                                       ggTTATCCATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGCCAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCAcg                                                  <  2:841411/239‑1 (MQ=255)
                                                                                                                                                                                           aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGcc                                              <  2:640143/239‑1 (MQ=255)
                                                                                                                                                                                              gATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCCCTCCCCCTTTCCCCCCGCCTTTTTTAAAAATGGGCTGCGCGCGGTGACCTGTTTTTTAACTACCGACCGCAAAAAAAGCATCCAAAAGGTGGAGGAAGGGGTGGAAGAAAAAAACCGCTTCTCCCCCTTGATTGATTAGGGGGAGGAAAAGAACAACAACCCGGATCTGGTGGGAGGACATATCAGGGCTCCCGCGCCgct                                           >  2:960102/1‑237 (MQ=255)
                                                                                                                                                                                                  aCCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTAcc                                       <  1:770124/239‑1 (MQ=255)
                                                                                                                                                                                                    cATGCCTCTCGGGCGTATAGCGGCGGGTCGCTCTCCACATTGCGCGACGAATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGt                                     <  1:992369/239‑1 (MQ=255)
                                                                                                                                                                                                            tcCGGCGTAGATCGGCGGGTCGCTCTCCCCATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATg                             <  2:406966/239‑1 (MQ=255)
                                                                                                                                                                                                            gccGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTGTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATg                             <  2:1268008/238‑1 (MQ=255)
                                                                                                                                                                                                            cccGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATg                             <  1:434749/238‑1 (MQ=255)
                                                                                                                                                                                                                gCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGc                                                                                       <  1:1263726/177‑1 (MQ=255)
                                                                                                                                                                                                                  gTATATCGGCGGTTCGCGCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTc                       <  1:953527/239‑1 (MQ=255)
                                                                                                                                                                                                                      atCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGga                   <  1:596948/239‑1 (MQ=255)
                                                                                                                                                                                                                       tCGGCGGGTCGCTCTCCTCATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTATTTTGAAACTACTGACCTTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCTTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGgag                  <  1:570800/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ctcCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGaa      >  1:340475/1‑238 (MQ=255)
                                                                                                                                                                                                                                    ctcCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAg     <  1:1251250/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ctcCAAATTGCGCGTCGGTTTTTTAAAATGTGGCCTGCTCGCGATGACCTTGTTTGATACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTGCTCCCGTTTGATTGATGAGTCGAAGAAAGCGACAAGCGAGCCGTATCTTGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAg     <  2:499375/239‑1 (MQ=255)
                                                                                                                                                                                                                                        cacaTTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCg  >  1:1008459/1‑239 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
CTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACGCATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCG  >  NC_000913/3019340‑3019809

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