Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,828,510 G→T 10.7% M258I (ATG→ATT)  srlD → sorbitol‑6‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,828,5100GT10.7% 233.7 / 2.6 81M258I (ATG→ATT) srlDsorbitol‑6‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (32/39);  new base T (2/8);  total (34/47)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

GCATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGTGCCTTCGACACACTATGCCAGCAAGGGCGGG  >  NC_000913/2828273‑2828738
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     
gCATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATg                                                                                                                                                                                                                                       <  2:939494/238‑1 (MQ=255)
 cATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGtt                                                                                                                                                                                                                                     >  1:13332/1‑239 (MQ=255)
      aCTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGtt                                                                                                                                                                                                                                     >  1:1190676/1‑234 (MQ=255)
                ggTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCCAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCACGGGATGTTCTTATCAACaaccgg                                                                                                                                                                                                                      >  1:821684/1‑234 (MQ=255)
                ggTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGGTGGTCTGATCAACaaccgg                                                                                                                                                                                                                      >  2:265475/1‑234 (MQ=255)
                  tAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGgag                                                                                                                                                                                                                     <  1:541835/238‑1 (MQ=255)
                  tAACCTGCTGAAATCGCCGATGTGCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTg                                                                                                                                                                                                                                  <  1:475534/225‑1 (MQ=255)
                   aaCCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGgaga                                                                                                                                                                                                                    <  1:18982/238‑1 (MQ=255)
                    aCCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATc                                                                                                                                                                                                                  <  1:485939/239‑1 (MQ=255)
                            aaaTCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGGTCTGGTTAACAACCgaggacccctaagg                                                                                                                                                                                                          >  2:252641/1‑228 (MQ=255)
                            aaaTCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GCACAGTCGTTCAATGTCACCGGCTGGTAGGAGATGATTTGCTCTACAGCAgagaatcactctag                                                                                                                                                                                                          >  2:1176345/1‑229 (MQ=255)
                              aTCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGAt                                                                                                                                                                                                       <  2:1180945/239‑1 (MQ=255)
                                   cGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCg                                                                                                                                                                                                  <  2:13332/239‑1 (MQ=255)
                                   cGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCg                                                                                                                                                                                                  <  1:1156363/239‑1 (MQ=255)
                                   cGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCg                                                                                                                                                                                                  <  1:1176345/239‑1 (MQ=255)
                                       gTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTgag                                                                                                                                                                                              <  1:183732/239‑1 (MQ=255)
                                         tCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGAc                                                                                                                                                                                            <  1:252641/239‑1 (MQ=255)
                                                              aCGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTa                                                                                                                                                                       <  2:550425/239‑1 (MQ=255)
                                                                  gACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGgag                                                                                                                                                                                                                     <  1:132755/190‑1 (MQ=255)
                                                                  aaCCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACg                                                                                                                                                                     <  2:624392/236‑1 (MQ=255)
                                                                     cAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCtcat                                                                                                                                                                <  1:479749/239‑1 (MQ=255)
                                                                       aGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCACGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGTATTTTCTGGTTAACAACGGGGAACCCTT‑AAGGAGTTCCCTGGGACCATAGCATGCCTGATGGGCTACGCccatca                                                                                                                                                              >  1:988790/1‑239 (MQ=255)
                                                                       aGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCCGGGGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGGTGGACCATAGAATGCCTGATGCGCTaaccccatcc                                                                                                                                                              >  1:761247/1‑230 (MQ=255)
                                                                       aGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTTTGATTCAAAGCGGGGAGCCCTT‑AAGGGGATCCGTGAGGCTATAGTATGCTCGAGACGCTACCCtcatca                                                                                                                                                              >  2:850772/1‑239 (MQ=255)
                                                                          tGGGTAACAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTGAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAg                                                                                                                                                             <  2:1064237/237‑1 (MQ=255)
                                                                            ggTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCCGGTGGTGTTCTTATCAACCGCCGGGATCCCTT‑aaa                                                                                                                                                                                                          >  2:132755/1‑189 (MQ=255)
                                                                              tATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGAt                                                                                                                                                                                                       >  1:420627/1‑191 (MQ=255)
                                                                              tATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGAt                                                                                                                                                                                                       <  2:420627/191‑1 (MQ=255)
                                                                                 cAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAgg                                                                                                                                                    <  2:341550/239‑1 (MQ=255)
                                                                                            aaGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCtcatc                                                                                                                                                               <  1:551200/217‑1 (MQ=255)
                                                                                                cGAGCAGGATCACAGCGACAAAGCACCGCACAAC‑CGCGGCTGCGATTAACAAGATGTGCAGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATAAATGTCACCGGCGGTCAGGTGATTTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACa                                                                                                                                     <  2:761247/239‑1 (MQ=255)
                                                                                                 gagaAGTATTACATCAACAAGGT‑CCCCTCAAACCGCGGCTGCGATAATCAAGTTGCGATGAATCTGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCCCC‑GGACAGTCGATAAATGTCCCCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACAt                                                                                                                                    <  2:1190676/239‑1 (MQ=255)
                                                                                                      agagggaacagACAATGTCCCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTAAATTGGCTGCTGTTCTCCCCCAGTCCTAAGGCGTCGTACTGCCCC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATTTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAgg                                                                                                                               <  2:1161384/229‑1 (MQ=255)
                                                                                                       tATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGgag                                                                                                                                                                                                                     <  1:515293/153‑1 (MQ=255)
                                                                                                       tATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTa                                                                                                                                                                       <  1:1092826/198‑1 (MQ=255)
                                                                                                        aTTACTGCGACAAAGTACCGCTCAAC‑CGGGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaa                                                                                                                             <  1:274212/239‑1 (MQ=255)
                                                                                                        aTTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCCCCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaa                                                                                                                             <  1:265475/239‑1 (MQ=255)
                                                                                                          tACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCCCCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGaaaa                                                                                                                           <  1:628179/239‑1 (MQ=255)
                                                                                                            cATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCTGATCAACAGCGGGGATCCATT‑AAGGGTCTCCGTGAGCATATAGAATGCCTGATGGGCTCCGCCCCTTAGGCACACAAGACTTCCGCCCCTACATTAAAGAAAAgg                                                                                                                         >  2:316427/1‑238 (MQ=255)
                                                                                                               cGAAAAGGTCCCGCTCAAC‑CGCGCATGCAGTTATAAAGTGGTCTTGAATGTGCTACTTTCCCACGCAATCCATAGGGCGTCGT‑CCGCACCGGGAAGGTCGATAAATGTCAGCGGGGTTAAGGTGATTTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCGGGGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTAt                                                                                                                      <  2:815675/239‑1 (MQ=255)
                                                                                                                gACAAAGTCCCGCTCAAA‑CGCGGCTGCGATTATCAAGATTTGCTGAATAGGCTTCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATg                                                                                                                     <  1:1048584/239‑1 (MQ=255)
                                                                                                                 aCAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGAt                                                                                                                                                                                                                                <  1:455318/132‑1 (MQ=255)
                                                                                                                 aCAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGAt                                                                                                                                                                                                                                >  2:455318/1‑132 (MQ=255)
                                                                                                                 aCAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCg                                                                                                                                                                                                  >  1:783778/1‑161 (MQ=255)
                                                                                                                 aCAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCg                                                                                                                                                                                                  <  2:783778/161‑1 (MQ=255)
                                                                                                                 aCAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGAc                                                                                                                                                                                            <  2:311945/167‑1 (MQ=255)
                                                                                                                 aCAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGAc                                                                                                                                                                                            >  1:311945/1‑167 (MQ=255)
                                                                                                                 aCAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCCGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATgg                                                                                                                    >  1:354725/1‑239 (MQ=255)
                                                                                                                 aCAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATgg                                                                                                                    >  1:707608/1‑239 (MQ=255)
                                                                                                                 aCAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGGTTTGATCAACAGCGGCGATCCATT‑AAGGATCTCCGTGGGACTATAGAATGCCCGATGCGCTACGCTCATCAg                                                                                                                                                             >  2:1092826/1‑198 (MQ=255)
                                                                                                                     aGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATc                                                                                                                                                                                                                  <  1:725769/142‑1 (MQ=255)
                                                                                                                     aGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCTGATCAACAGCGGAGATc                                                                                                                                                                                                                  >  2:725769/1‑142 (MQ=255)
                                                                                                                       tACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAg                                                                                                                                                             <  1:65887/192‑1 (MQ=255)
                                                                                                                           gCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCCCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCAc                                                                                                           >  1:1060421/1‑238 (MQ=255)
                                                                                                                            cTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCCCATCAGGCATACAGGACTTCCGCCCCTACATTAAGGAAAAGTTATGGGATCCGCACtc                                                                                                         >  1:1150773/1‑239 (MQ=255)
                                                                                                                            cTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATTTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTAtt                                                                                                                     >  1:1161384/1‑226 (MQ=255)
                                                                                                                            cTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCTGATCAACAGCGGGGGTCCATT‑AAGGATCTCCGGGAGACTTTAGAATGCCTGATTGGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACtc                                                                                                         >  1:331128/1‑239 (MQ=255)
                                                                                                                                 a‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGAt                                                                                                                                                                                                       >  1:954213/1‑140 (MQ=255)
                                                                                                                                 a‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGAt                                                                                                                                                                                                       <  2:954213/140‑1 (MQ=255)
                                                                                                                                 a‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGAACCATT‑AAGGATCTCCGTGGGACTATAGAATGCCCGATGCGCTAAGCTCATCAGGCATACAggg                                                                                                                                                   >  2:65887/1‑191 (MQ=255)
                                                                                                                                                    agttgtgATGAATATGCTGCTTTTCTACGCCAGCCCTAAGGCGTCGTACTGCACG‑GGACGGTCGATAAATGTCACCGGGGGTAGGGTGATTTTCTGATCAACGGCGGAGATCCATT‑AAGGATCTCGGTGAGCCTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTg                                                                                  <  2:1051985/236‑1 (MQ=255)
                                                                                                                                                                        tgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGGTTGCGCACCTGGCCTTAg                                                              >  1:327897/1‑239 (MQ=255)
                                                                                                                                                                          ttCTACGCCAGTCCTAAGGGGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGGGGAGACCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATcc                                                                                                               <  1:956704/188‑1 (MQ=255)
                                                                                                                                                                            cTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCTGATCCACAGCGGAGATCCATT‑AAGGGTCTCCGTGGGACTATAGAATGCCTGATTCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGGTATGGTATCCGCACTCATCACCGTCGCCCGTATCGCCTGGTGTGCGCCACAGGCCtttagcg                                                           >  2:942814/1‑233 (MQ=255)
                                                                                                                                                                                    gTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCCGGTGGTGTTCTGATCAACAGCGGGGATCCCTT‑AAGGATCTCCGTGGGACTATAGAAAGCCTGATGCGCTACGCTCATCAGGGATACATGACTTCCGCCCCCACATTAAGGAAAAATTATTGTAACCGccaccatta                                                                                                     >  2:956704/1‑180 (MQ=255)
                                                                                                                                                                                     tCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTAt                                                                                                                      <  1:1019261/170‑1 (MQ=255)
                                                                                                                                                                                               tcgtACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGGTATCGCCTGGTGTGGGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGtttt                                       >  1:756397/1‑239 (MQ=255)
                                                                                                                                                                                                     tccccc‑GGACATTCGATCAATGTCACCGGCGGTCAGGTTATGTTCTGATCAACGGCGGAGATCCTTT‑AAGGATCTCCGTGAGCCTATAAAATGCCTGATGCCCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCg                                                                                            <  2:966410/176‑1 (MQ=255)
                                                                                                                                                                                                     tGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATTTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGt                                 <  1:862189/239‑1 (MQ=255)
                                                                                                                                                                                                     tGCACA‑GGACGTTCGATAACTGTCACGGGCGGTCAGGTGATGTTCTGATCAAGAGCGGAGATCAATT‑AAGGATCTCCGTGAGACTAAAAAATGCCTGATGCGCTACCCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGt                                 <  2:756397/239‑1 (MQ=255)
                                                                                                                                                                                                      gcaacagaaccgTCGATCAATGTCCCCGGCGTTCAGGTGATTTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCGGTGAGCCTATAGAATGCCTGATGCGCTACGCTCATCGGGCATACAGGCCTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGt                                 <  2:354725/228‑1 (MQ=255)
                                                                                                                                                                                                      gCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGGGATGTTCTGATCAACAGCGGAGATCCATT‑TAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCCCATCAGGCATACAGGGCTTTCCCCCCTACATTAAGGGAAAGTTATGGTATCCGCACTCCTCACCCTCCCCGTTTTCCCCTTGGGTGGGCAACTGGGCCTAGGCGGCTGGGCAAATTCTCGGTTTTATCGt                                 >  2:862189/1‑238 (MQ=255)
                                                                                                                                                                                                       cccc‑GGACAGTCGTTCAATGTCCCCGGCGTTCAGGTATTTTTCTAACCAACGGCGGACACCCATT‑AAGGATCCCCGTGAGACTATAGAATCCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAAttt                                               <  2:734121/221‑1 (MQ=255)
                                                                                                                                                                                                        ccc‑GGAAAGTCGATCAATGCCCCCGGCGGTAAGGTGATTTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGTGc                               <  1:316427/237‑1 (MQ=255)
                                                                                                                                                                                                          c‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGGTTTTACCGTGCCt                             >  1:973157/1‑238 (MQ=255)
                                                                                                                                                                                                                aGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTg                                                                                  >  1:966410/1‑180 (MQ=255)
                                                                                                                                                                                                                  tCGATCAATGTCACCGGCGGTCAGGGGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGGGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTaa                                     >  1:734121/1‑223 (MQ=255)
                                                                                                                                                                                                                   cGATCAATGTCACCGGGGGTCAGGTGATGTTCTGATCAACAGCGGAGTTCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGTGCCTTCGa                         <  1:742969/234‑1 (MQ=255)
                                                                                                                                                                                                                     aTCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGCCTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGtt                                                                                          <  2:605261/167‑1 (MQ=255)
                                                                                                                                                                                                                     aTCAATGTCACCGGCGGTCAGGGGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGtt                                                                                          >  1:605261/1‑167 (MQ=255)
                                                                                                                                                                                                                                cggcgtTCAGGATATGTTCATAACAAC‑GGGGAACTCCAATAAAGGCTCTCGGGGAAACTAAAAAATCCCTGATGCCCTACGCTCATCGGGCAAACAGGCCTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGTGCCTTCGACACACTATGCCAGCAAgg       <  2:331128/239‑1 (MQ=255)
                                                                                                                                                                                                                                      taagttGATGTTCGGATCACCAGCGGAGATCCATT‑AAGGATCTCCGGGAGACTATAGAATGCCTGATGCGCAACGCTCATCAGGCATACAGGCCTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGTGCCTTCGACACACTATGCCAGCAAGGGCggg  <  1:850772/233‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     
GCATTCACTGATGCTCGGTAACCTGCTGAAATCGCCGATGTTCCAGTCACTGTTGCCACAATACGCGACCAAGCTGGGTATCAAACCGGATCAAGTCGAGCAGTATTACATCGACAAAGTACCGCTCAAA‑CGCGGCTGCGATTATCAAGATGTGCTGAATATGCTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACC‑GGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATT‑AAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGTGCCTTCGACACACTATGCCAGCAAGGGCGGG  >  NC_000913/2828273‑2828738

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑