Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,306,540 T→C 6.6% S20G (AGC→GGC)  yjcS ← metallo‑beta‑lactamase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,306,5400TC 6.6% 382.9 / 11.4 121S20G (AGC→GGC) yjcSmetallo‑beta‑lactamase superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (56/57);  new base C (4/4);  total (60/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 121 total aligned reads displayed.

GATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATT  >  NC_000913/4306397‑4306685
                                                                                                                                               |                                                                                                                                                 
gATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcgg                                                                                                                                              >  1:229793/1‑149 (MQ=255)
gATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTggcgg                                                                                                                                              >  4:537007/1‑149 (MQ=255)
  tAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCa                                                                                                                                            >  1:783375/1‑149 (MQ=255)
    aaTCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGt                                                                                                                                          >  4:816698/1‑149 (MQ=255)
     aTCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatga                                                                                                                                         <  1:836438/149‑4 (MQ=255)
      tCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTCGCGGTATTTACCATCATTCCGCTGGCGGCAGTTa                                                                                                                                        >  2:362365/1‑149 (MQ=255)
      tCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatgat                                                                                                                                        <  1:949807/149‑5 (MQ=255)
      tCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatgat                                                                                                                                        >  2:949807/1‑145 (MQ=255)
       ccACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTaa                                                                                                                                       <  2:858725/149‑1 (MQ=255)
       ccACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatgatg                                                                                                                                       <  2:814515/149‑6 (MQ=255)
        cACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAAc                                                                                                                                      <  3:498859/149‑1 (MQ=255)
        cACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAAc                                                                                                                                      <  4:713817/149‑1 (MQ=255)
        cACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAAc                                                                                                                                      <  3:816698/149‑1 (MQ=255)
                    cGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATc                                                                                                                          <  2:783375/149‑1 (MQ=255)
                    cGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATc                                                                                                                          <  2:298867/149‑1 (MQ=255)
                    cGTCAAAATCCCGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATc                                                                                                                          <  3:682929/149‑1 (MQ=255)
                     gTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATcc                                                                                                                         >  2:892443/1‑149 (MQ=255)
                     gTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATcc                                                                                                                         >  1:349089/1‑149 (MQ=255)
                      tCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGCTTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatga                                                                                                                                         <  2:466949/132‑4 (MQ=255)
                      tCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCCAGCGGCatga                                                                                                                                         >  1:466949/1‑129 (MQ=255)
                      tCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGGTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCt                                                                                                                        >  1:408919/1‑149 (MQ=255)
                       cAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTg                                                                                                                       <  4:257869/149‑1 (MQ=255)
                        aaaaTCCTGACGATCGCTAAATGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCCTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGc                                                                                                                      <  3:755611/149‑1 (MQ=255)
                             ccTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCaaa                                                                                                                 <  2:419734/149‑1 (MQ=255)
                              cTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAAc                                                                                                                <  3:89892/149‑1 (MQ=255)
                              cTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAAc                                                                                                                <  2:63450/149‑1 (MQ=255)
                               tGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGCGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCACGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACg                                                                                                               >  1:870774/1‑149 (MQ=255)
                               tGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACg                                                                                                               >  2:1017101/1‑149 (MQ=255)
                                  cGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGaa                                                                                                            <  4:506228/149‑1 (MQ=255)
                                    aTCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATa                                                                                                          >  2:583961/1‑149 (MQ=255)
                                              ggTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGtt                                                                                                <  4:614402/149‑1 (MQ=255)
                                               gTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTa                                                                                               >  3:771540/1‑149 (MQ=255)
                                                tAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTAt                                                                                              <  1:69615/149‑1 (MQ=255)
                                                 aaTGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTAtt                                                                                             >  2:1016284/1‑149 (MQ=255)
                                                  aTGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTc                                                                                            <  2:1029869/149‑1 (MQ=255)
                                                   tGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCa                                                                                           <  3:595104/149‑1 (MQ=255)
                                                    gATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCat                                                                                          <  4:600587/149‑1 (MQ=255)
                                                    gATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCat                                                                                          >  3:797652/1‑149 (MQ=255)
                                                         ggCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATaaa                                                                                     <  2:967522/149‑1 (MQ=255)
                                                          gCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAt                                                                                    >  4:32184/1‑149 (MQ=255)
                                                             tAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAtctc                                                                                 >  4:1011401/1‑149 (MQ=255)
                                                               aTTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCca                                                                               <  4:170649/149‑1 (MQ=255)
                                                                ttCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcat                                                                              >  1:838738/1‑149 (MQ=255)
                                                                ttCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcat                                                                              >  3:529536/1‑149 (MQ=255)
                                                                   tGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcattta                                                                           <  3:730770/149‑1 (MQ=255)
                                                                   tGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcattta                                                                           <  3:568100/149‑1 (MQ=255)
                                                                     aTTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatc                                                                         >  4:92935/1‑149 (MQ=255)
                                                                     aTTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatc                                                                         <  2:617534/149‑1 (MQ=255)
                                                                      ttAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatca                                                                        <  2:456813/149‑1 (MQ=255)
                                                                         aCCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcattt                                                                     <  3:865934/149‑1 (MQ=255)
                                                                          ccTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcattta                                                                    <  3:611228/149‑1 (MQ=255)
                                                                           cTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcatttat                                                                   <  2:388870/149‑1 (MQ=255)
                                                                              ttGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTc                                                                <  1:780916/149‑1 (MQ=255)
                                                                              ttGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTc                                                                <  3:218342/149‑1 (MQ=255)
                                                                              ttGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTc                                                                <  4:283717/149‑1 (MQ=255)
                                                                              ttGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTACCATTTATGTc                                                                <  1:7773/149‑1 (MQ=255)
                                                                                gggTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgaga                                                                                                   >  4:302236/1‑112 (MQ=255)
                                                                                gggTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCgaga                                                                                                   <  3:302236/112‑1 (MQ=255)
                                                                                 ggTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCaa                                                                                                                                  <  3:743718/80‑1 (MQ=255)
                                                                                 ggTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCaa                                                                                                                                  >  4:743718/1‑80 (MQ=255)
                                                                                    atatTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAat                                                          <  3:18195/149‑1 (MQ=255)
                                                                                     tatTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAata                                                         <  3:720127/149‑1 (MQ=255)
                                                                                      atTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAatat                                                        >  2:925558/1‑149 (MQ=255)
                                                                                       ttGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAAtctc                                                                                 >  2:825960/1‑123 (MQ=255)
                                                                                          aGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcatttat                                                                   <  2:469411/134‑1 (MQ=255)
                                                                                          aGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcatttat                                                                   >  1:469411/1‑134 (MQ=255)
                                                                                            tGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATc                                                  <  4:858800/149‑1 (MQ=255)
                                                                                            tGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATc                                                  >  4:459167/1‑149 (MQ=255)
                                                                                             ggcggcTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATcc                                                 >  2:1031211/1‑149 (MQ=255)
                                                                                               cggcTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCAt                                               <  1:1049311/149‑1 (MQ=255)
                                                                                                ggcTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGt                                                              <  3:61747/133‑1 (MQ=255)
                                                                                                ggcTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGt                                                              >  4:61747/1‑133 (MQ=255)
                                                                                                     tCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCaaaa                                         <  1:622503/149‑1 (MQ=255)
                                                                                                         ttCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACg                                     >  4:121343/1‑149 (MQ=255)
                                                                                                          tCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATATAATCCATGCAAAATACGt                                    <  1:925558/149‑1 (MQ=255)
                                                                                                           ccTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATcc                                                 >  2:750232/1‑135 (MQ=255)
                                                                                                           ccTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATcc                                                 <  1:750232/135‑1 (MQ=255)
                                                                                                           ccTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTc                                   >  4:452415/1‑149 (MQ=255)
                                                                                                            ctctTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCa                                  >  1:139998/1‑149 (MQ=255)
                                                                                                               tttCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCAttt                               >  3:875543/1‑149 (MQ=255)
                                                                                                               tttCGCGCTAGCGGTATTTTCCATCATGCCGCTGGCGGCAGTTAACGCAATATCAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATATCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCAttt                               >  4:661983/1‑149 (MQ=255)
                                                                                                                  cgcgTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTaaa                            <  2:392193/149‑1 (MQ=255)
                                                                                                                   gcgTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAg                           >  3:200476/1‑149 (MQ=255)
                                                                                                                      ttAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTg                        >  4:1035269/1‑149 (MQ=255)
                                                                                                                          cGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTaaa                            <  1:471064/141‑1 (MQ=255)
                                                                                                                          cGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTaaa                            >  2:471064/1‑141 (MQ=255)
                                                                                                                            gTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcatttat                                                                   <  3:350172/100‑1 (MQ=255)
                                                                                                                            gTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCcatttatcatttat                                                                   >  4:350172/1‑100 (MQ=255)
                                                                                                                            gTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTa                  <  3:1031375/149‑1 (MQ=255)
                                                                                                                            gTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTa                  >  3:662457/1‑149 (MQ=255)
                                                                                                                               tttACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCa               >  3:182667/1‑149 (MQ=255)
                                                                                                                                 tACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATa             >  1:33522/1‑149 (MQ=255)
                                                                                                                                   gttaactGCCGCCAGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAAt           >  4:346799/7‑149 (MQ=255)
                                                                                                                                    catcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAAtt          >  2:181133/1‑149 (MQ=255)
                                                                                                                                     atcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTg         >  4:473275/1‑149 (MQ=255)
                                                                                                                                      tcatGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAACCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTgg        <  2:650770/149‑1 (MQ=255)
                                                                                                                                            cGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAAtt  >  3:299907/1‑149 (MQ=255)
                                                                                                                                            cGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAAtt  >  3:835788/1‑149 (MQ=255)
                                                                                                                                             gCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATa                                       >  3:266418/1‑111 (MQ=255)
                                                                                                                                             gCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATa                                       <  4:266418/111‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GATAAATCCACGCTGGGCATCGTCAAAATCCTGACGATCGCTAAACGGTAATGATTTGGCGTAATTCTGATTAACCTGTTGGGTATATTGAGTGGCGGCTTTCGCTTCCTCTTTCGCGTTAGCGGTATTTACCATCATGCCGCTGGCGGCAGTTAACGCAATAACAATCCTGCTCAAACGGAATAACCGAGAGTTATTCATATAAATCTCCATTTATCATTTATGTCGTAAATATGTAATCCATGCAAAATACGTCATTTAAAGCTGGCGTTATCATAATTGGGAAATT  >  NC_000913/4306397‑4306685

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