Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,533 C→T Y213Y (TAC→TAT)  serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,5330CT96.4% 162.1 / ‑2.3 57Y213Y (TAC→TAT) serBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (29/25);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

CCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑TGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAA  >  NC_000913/4625286‑4625776
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
ccAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTg                                                                                                                                                                                                                                                  <  1:185881/251‑1 (MQ=255)
 aaaaCTGGCCGGAACGGGGGAGAGGGGGGGGGAAGTACCCGAACGGGGGAGGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTgc                                                                                                                                                                                                                                                 <  2:60887/250‑1 (MQ=255)
               cGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTgcgc                                                                                                                                                                                                                                   >  1:245572/1‑251 (MQ=255)
                gggCGAGAGGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGcc                                                                                                                                                                                                                                  <  1:449862/251‑1 (MQ=255)
                    gagaTGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGAAAAGTTGCGCCTGAc                                                                                                                                                                                                                              >  2:189852/1‑251 (MQ=255)
                                gaagaaccctaaCGGGCGATGCGCGGCGACTTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATCCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCGGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGc                                                                                                                                                                                                                  <  1:485913/241‑1 (MQ=255)
                                             gggggattcgcgggcaacccGATTTTACCGCAAGCCTGCGCAGCGG‑TGTGGCGACGCTGAAAGGCGCTGCCGC‑CAATATTCTGCAACAGGTGCGAGAAAATCGGCCGCTGATGCCATGCTTAACGCATCTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTTACCGCCGTGGTAGCCAATGAACTGgag                                                                                                                                                                                                     <  2:312086/232‑1 (MQ=255)
                                             gggCGAGGCGCGGAAACCTCGTTTTTCCCGCCAGCCGGCGCAGCCG‑CGTGGCGACGCGGAAATGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTACCGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCTGTGGCCTCCGGCGGCTTTACTTTCTTGGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGgag                                                                                                                                                                                                     <  2:33711/251‑1 (MQ=255)
                                                cGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTgg                                                                                                                                                                                                       >  1:92325/1‑246 (MQ=255)
                                                     cgcgGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTATGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATgg                                                                                                                                                                                              >  1:22900/1‑250 (MQ=255)
                                                        ggcgAACTCGTTTTTCCCGCCAGCCTGCGCAGCGGCTG‑GGCGACGCTGAAAGGCGTTGA‑GCACAATATTCTGAAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGTTTAACGCAACTGGTGC‑CTAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGTTTTACTTTCTTTGCTGAATATCTGCGCGAAAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACgg                                                                                                                                                                                          <  2:436564/251‑1 (MQ=255)
                                                        ggcgAACTCGATTTTCCAGCCAGCCGGCGCAGCCG‑CGTGGCGACGCGGAAAGGCGTTGACGC‑CAATATTCTGCAAAAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCGGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCCGAAGATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGCCgg                                                                                                                                                                                          <  2:153471/251‑1 (MQ=255)
                                                           gAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACggaa                                                                                                                                                                                        >  1:372213/1‑248 (MQ=255)
                                                                cGATTTTACCGCCAGCCTGCGCAGCCG‑CGTGGGGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAAtt                                                                                                                                                                                    <  1:321889/249‑1 (MQ=255)
                                                                 gtttttACCGCCAGCCTGCGAAGCCG‑CGTGGCGACGCTGAAAGTGGCGGACGC‑CAATTTTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGGGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAc                                                                                                                                                                                 <  2:157111/249‑1 (MQ=255)
                                                                  aTTTTACCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTAcc                                                                                                                                                                                >  2:346183/1‑251 (MQ=255)
                                                                      aaCCGCCAGCCTGCGCAGCCG‑CGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGc                                                                                                                                                                             <  1:410689/249‑1 (MQ=255)
                                                                        ccGCCAGTCTGCGCAGCCT‑CGGGGCGACGCGGAAAGGCGCTGACGC‑CAATATCCTGCAACAGGTGCGTGAAACTCTGCCGATTATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCTCGACAAGTTGCGCCTTACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAAt                                                                                                                                                                          <  1:312712/251‑1 (MQ=255)
                                                                                    gcAGCCG‑CGTGGAGACGATGAAAGGAGCTGACGC‑CCATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGATGGAAAGTGGCGATTGCCTCCGGCGGATTTACTTTCTTTGCTGAATATCTGCGCGAACGGTTGAGCCTGAACGCAGTGGTAGACAATGAACTGGGGATCATGGACGGTAAATTTCCCGGGAATGTGAGCGGGGAc                                                                                                                                                              >  2:276662/1‑251 (MQ=255)
                                                                                                 cGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGAcgcg                                                                                                                                                  >  2:257626/1‑251 (MQ=255)
                                                                                                   aCGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACgcgc                                                                                                                                                 >  1:407433/1‑250 (MQ=255)
                                                                                                     gCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATAAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGt                                                                                                                                              >  2:53541/1‑251 (MQ=255)
                                                                                                               gcTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGa                                                                                                                                     >  1:305748/1‑250 (MQ=255)
                                                                                                                  gACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACAGCGTAGACGCGCAGTACAAAGCGaaa                                                                                                                                   <  1:412162/249‑1 (MQ=255)
                                                                                                                   ttgaccatttttCTGCAACGGGTGCGTTAAATACGGCCGCGGATGCCAGGTTTAACGCAACTG‑TGCCTCAAGCTGGAAACGCTGTGTTGGAAAGTGGCTATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCGGCGCGACAAGTTGCCCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAGGCGaaaa                                                                                                                                  <  2:347024/240‑1 (MQ=255)
                                                                                                                      c‑cAATATTCTGCAGCAGGTGCGTGAAAATCGGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGTCTGGAAAGTGGCGCTTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGATGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGCCGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGGACAAAGCGAAAACTCTg                                                                                                                             >  2:263859/1‑251 (MQ=255)
                                                                                                                         aaTATTCTGCACCAGGCGCGTGAAAATCTGCCGCTAATGCCAGGCTTAACGCAACTGGTGC‑TCACGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTAATTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGAGCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAACACTCTGAc                                                                                                                           >  2:248169/1‑251 (MQ=255)
                                                                                                                         aaTATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGCCGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGAAGCGCAGTACAAAGCGAAAACTCTGAc                                                                                                                           >  1:239874/1‑251 (MQ=255)
                                                                                                                         aaTATTCCTCAACAGGTTCGGGAAAATCTGCCGCTGAGGCCAGGCTTAACGCAGCTGGGGC‑TCAAGCGTGAAGCGCTGGGCTGGAAAGTGGCGATTGCCGCCGGCGGTTTTACTTTGTTGGCTGAATGTCTGCGCGGTAGGTTGCGCCTGGCCGCCGTGGTAGCCAATGAACTGGATATCATGGACGGTAAATTTGCCTGCAATGTGATCGGCGACATCTTAGACGCGCAGTACAAAGCGAAAACTCTGAc                                                                                                                           <  1:210669/251‑1 (MQ=255)
                                                                                                                            atTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACaaa                                                                                                                                         >  1:104976/1‑234 (MQ=255)
                                                                                                                            atTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGActcg                                                                                                                        <  1:272965/251‑1 (MQ=255)
                                                                                                                            atTCGGCAACAGGTGCGGGAAAATCGGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACaaa                                                                                                                                         <  2:104976/234‑1 (MQ=255)
                                                                                                                                      aggtcggtgcaaaTCTGCCGCTGATGCCAGGCTTAACGCACCTTGTGC‑TCAAGCTGGAAACGCTTGGCTTGAAAGTGGCTATTGCCTCCGGCGGCTTTTCTTTCTTTGCTGAATATCTGCGCTACAAGTTGCGCCTGTCCGCCGTGGTAGCCAATGAACTGTAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAg                                                                                                              <  1:266960/241‑1 (MQ=255)
                                                                                                                                      aGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCgcgc                                                                                                                >  2:478142/1‑249 (MQ=255)
                                                                                                                                            gTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTAt                                                                                                        >  1:72649/1‑251 (MQ=255)
                                                                                                                                                               gCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGa                                                                                     <  1:132375/251‑1 (MQ=255)
                                                                                                                                                                 aaGGCTTAACGCAACTGGGGT‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGAcc                                                                                   <  1:12268/250‑1 (MQ=255)
                                                                                                                                                                   ggCTTAACGCAACTGGTGC‑TCAAGCTGGAACCGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGCACTGGAGATCATGGACGGTAAATTTACCGGCACTGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCATCGCGCAGGGGTATGAAATCCCGCTGgcgacgaccgg                                                                                 >  2:491361/1‑243 (MQ=255)
                                                                                                                                                                            cAACTGGGGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGGAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgattgg                                                                        >  2:9054/1‑251 (MQ=255)
                                                                                                                                                                                 ggTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGGAAATTTACCGGCAATGAGATCGGCGCCATCGTAGACGCGCAGTACAAATCGAAAAATCTGACACGCCTCGCGCCGGAGTATGCAACCCCGCCGgcgccgcgcgaggcgattagcgagg                                                                   >  2:408774/1‑230 (MQ=255)
                                                                                                                                                                                  ggggcgcaggCTGGAACCGCGGGGCTGGAAGGGGGCGATGGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATg                                                                   <  2:108666/242‑1 (MQ=255)
                                                                                                                                                                                     c‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGc                                                               >  1:261114/1‑251 (MQ=255)
                                                                                                                                                                                         aaGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCa                                                             >  1:384471/1‑250 (MQ=255)
                                                                                                                                                                                         aaGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgattggc                                                                       >  1:136438/1‑240 (MQ=255)
                                                                                                                                                                                         aaGCGGGAAACGCTGGGCGGGAAAGGGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgattggc                                                                       <  2:136438/240‑1 (MQ=255)
                                                                                                                                                                                          agctggaaacgttgggcgagaggtgggCGATGGCCTCCGGCGGCTTTCCTTGCTTTGCTAAATCTCTGCGCGGCAAGTTGCGCCTGACCGCCGGGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAATCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAAGCCCGCTGGCGCAGACCGTGGCGATTGGCTATGGAGCCAAt                                                           <  2:284756/226‑1 (MQ=255)
                                                                                                                                                                                           gCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATg                                                          <  1:362486/251‑1 (MQ=255)
                                                                                                                                                                                              ggAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGa                                                         <  1:478142/249‑1 (MQ=255)
                                                                                                                                                                                               gAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGAGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGAACGTGGCGATTGGCGATGGAGCCAATAAcc                                                       >  1:220057/1‑250 (MQ=255)
                                                                                                                                                                                                      tgggctggAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGAACTGCCgatg                                               >  1:498824/1‑251 (MQ=255)
                                                                                                                                                                                                                              ccGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAGAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGACATCCCGCAGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCGGGTCTGGGGAtt                       >  2:8860/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tttACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGACAATAACCTGCAGATGATCAAAGCGGAAGGGCGgggggatgccaaccc               >  2:447705/1‑240 (MQ=255)
                                                                                                                                                                                                                                        tACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATg             >  1:176259/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ctttctttGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCCTgcc           >  2:231920/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ttctttGCTGAATAGCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgc            <  1:151433/248‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgccaagccaa    >  1:103094/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     gAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCaaaa  <  1:346183/249‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
CCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑TGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAA  >  NC_000913/4625286‑4625776

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