Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,546 C→T L218L (CTG→TTG)  serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,5460CT96.3% 160.1 / ‑2.1 55L218L (CTG→TTG) serBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (28/24);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

CGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑TGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAA  >  NC_000913/4625301‑4625776
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
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)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      
CGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCG‑TGTGGCGACGCTGAAAGGCGCTGACGC‑CAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGC‑TCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAA  >  NC_000913/4625301‑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: ‑