Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,555 A→T 17.2% intergenic (+44/‑176) metH → / → yjbB homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,227,5550AT17.2% 298.2 / 54.6 122intergenic (+44/‑176)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base A (60/41);  new base T (8/13);  total (68/54)
Fisher's exact test for biased strand distribution p-value = 9.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

Only 100 of 124 total aligned reads displayed.

CGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGA  >  NC_000913/4227457‑4227653
                                                                                                  |                                                                                                  
cGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                                   >  2:742226/1‑100 (MQ=255)
  ttACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa                                                                                                 <  1:674205/100‑1 (MQ=255)
   tACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg                                                                                                >  2:2060961/1‑100 (MQ=255)
     ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                              >  1:1102257/1‑100 (MQ=255)
     ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                              >  2:2408917/1‑100 (MQ=255)
      cGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGc                                                                                             >  1:245029/1‑100 (MQ=255)
       gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGcc                                                                                            <  1:1037993/100‑1 (MQ=255)
        aaGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGTACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                           <  1:390730/100‑1 (MQ=255)
          gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                         >  1:2807779/1‑100 (MQ=255)
          gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                         <  1:165145/100‑1 (MQ=255)
          gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                         >  2:404842/1‑100 (MQ=255)
          gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                         >  2:1377388/1‑100 (MQ=255)
          gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                         <  1:1653966/100‑1 (MQ=255)
          gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                         <  1:659400/100‑1 (MQ=255)
           ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                        <  2:667031/100‑1 (MQ=255)
            tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt                                                                                       <  2:1679019/100‑2 (MQ=39)
            tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                        <  1:1042155/99‑1 (MQ=39)
             gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGttt                                                                                      <  1:622194/100‑3 (MQ=38)
             gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGttt                                                                                      <  1:1569080/100‑3 (MQ=38)
              aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtctc                                                                                     >  1:2811022/1‑100 (MQ=255)
               gggctggctgGCACCGAATCTGGGGGATGACGCTGACGGATTCACAAATCTGTCACTTTCCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                           <  1:2429879/91‑1 (MQ=255)
               gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTctct                                                                                    <  2:2567753/100‑1 (MQ=255)
               gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTctct                                                                                    >  1:118882/1‑100 (MQ=255)
                 gctggctgGCACCGAGTGTCGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTGCCTTACAACAAACAGGGCGCTGAATGAGTGCCCTGGCTCgtt                                                                                  <  1:1203242/100‑1 (MQ=255)
                 gctggctgGCACCGAATCTGGGGTATGACGCTGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                                  <  2:703254/100‑1 (MQ=255)
                 gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCttt                                                                                  >  2:1088682/1‑100 (MQ=255)
                 gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                                  <  2:1189494/100‑7 (MQ=37)
                  ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                 <  1:1058346/100‑1 (MQ=255)
                  ctggctggCACCGAATCTGGCGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgt                                                                                 <  1:2298535/100‑8 (MQ=37)
                   tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgta                                                                                <  1:2120012/100‑9 (MQ=37)
                      ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                             >  2:677990/1‑100 (MQ=255)
                       tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAat                                                                            <  1:2340789/100‑1 (MQ=255)
                        ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                           >  2:2460992/1‑100 (MQ=255)
                        ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                           <  2:1986996/100‑1 (MQ=255)
                          cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                         <  1:487768/100‑1 (MQ=255)
                           aCCGAATCTGGGGTATGACGCGGACTGATTAACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGa                                                                        >  2:2548375/1‑100 (MQ=255)
                            ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATgn                                                                        <  1:2612033/99‑2 (MQ=255)
                              gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                     >  1:262277/1‑100 (MQ=255)
                                      ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                             >  2:2253018/1‑100 (MQ=255)
                                       gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                            >  1:439181/1‑100 (MQ=255)
                                         aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                          >  2:2122024/1‑100 (MQ=255)
                                           gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                        <  2:2466930/100‑1 (MQ=255)
                                              gcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                     <  2:1868255/100‑1 (MQ=255)
                                               cgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACan                                                    <  2:1147359/100‑2 (MQ=255)
                                               cgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                                    >  1:2252629/1‑100 (MQ=255)
                                                   cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                                >  2:1388374/1‑100 (MQ=255)
                                                    tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                               <  2:383336/100‑1 (MQ=255)
                                                         cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                          >  2:1495697/1‑100 (MQ=255)
                                                         cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                          >  2:926456/1‑100 (MQ=255)
                                                          acaAATCTGTCACTTTTCCTTACAATAAACAGGGCGCTCAATGAGTGCCCTGTCTCATTATTAATATGATACACTTATACTGGAAACAGTCTGGgataa                                          >  1:2122492/1‑99 (MQ=255)
                                                           caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTTTCTCTTTATTAATATGAAACACTTATACTCGAAACAGGCTGGAACAATc                                        >  1:2201937/1‑100 (MQ=255)
                                                           caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                        <  1:88565/100‑1 (MQ=255)
                                                           caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                        >  2:2485118/1‑100 (MQ=255)
                                                           caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                        <  2:270789/100‑1 (MQ=255)
                                                             aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                      <  1:389219/100‑1 (MQ=255)
                                                              aTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                     <  2:20452/100‑1 (MQ=255)
                                                               tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                                    <  2:118883/100‑1 (MQ=255)
                                                               tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                                    <  1:718859/100‑1 (MQ=255)
                                                                 tGTCACTTTTCCTTACAACAAACAGGGTGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGcc                                  >  2:2024331/1‑100 (MQ=255)
                                                                  gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGCACAAGGATGGCCTAATCtt                                      >  2:2717710/1‑95 (MQ=255)
                                                                  gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                                 >  2:1626467/1‑100 (MQ=255)
                                                                   tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                                >  1:1461250/1‑100 (MQ=255)
                                                                   tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                                <  1:971262/100‑1 (MQ=255)
                                                                    cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                                 >  1:984854/1‑98 (MQ=255)
                                                                    cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                               >  2:1629205/1‑100 (MQ=255)
                                                                     aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                              >  1:1050320/1‑100 (MQ=255)
                                                                      cTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                             >  2:111482/1‑100 (MQ=255)
                                                                         ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                          <  2:354975/100‑1 (MQ=255)
                                                                         ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                          <  2:2743132/100‑1 (MQ=255)
                                                                          tCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                         <  2:1167105/100‑1 (MQ=255)
                                                                            cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                       >  2:1595897/1‑100 (MQ=255)
                                                                             ttACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                      <  2:558849/100‑1 (MQ=255)
                                                                             taacaacaAACAGGGCGGTCAATGAGTGCCCTGTCTCGTTAGTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTCGCCGGAAAAACGAAt                      <  1:2460981/98‑1 (MQ=255)
                                                                              tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                     >  1:2398407/1‑100 (MQ=255)
                                                                              tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                     <  1:327833/100‑1 (MQ=255)
                                                                               taaagagacaGAGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGACTAATCTTAGCCGGAAAAACGAATg                     >  2:364320/8‑99 (MQ=37)
                                                                               acaacaAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt                    >  2:182815/1‑100 (MQ=255)
                                                                                tctcaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                   <  2:2446860/97‑1 (MQ=255)
                                                                                gaagagagaGGGCACTGAATGAGTTACGTGTCTCTTTAATAATAAGAAACACGTAAAATTGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                       <  1:2502141/88‑1 (MQ=255)
                                                                                caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                   >  1:1325998/1‑100 (MQ=255)
                                                                                caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                   >  1:1890462/1‑100 (MQ=255)
                                                                                 aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa                  >  2:228104/1‑100 (MQ=255)
                                                                                  acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa                 >  1:483592/1‑100 (MQ=255)
                                                                                   gagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATTCTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa                >  1:2874883/4‑100 (MQ=38)
                                                                                   caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaga                >  1:2446849/1‑98 (MQ=255)
                                                                                    agacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa               >  1:1896294/3‑100 (MQ=39)
                                                                                     gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg              >  2:2566065/2‑100 (MQ=39)
                                                                                     gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg              >  2:1314504/2‑100 (MQ=39)
                                                                                       caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg              >  2:629127/1‑98 (MQ=255)
                                                                                       caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc            >  1:1624389/1‑100 (MQ=255)
                                                                                       caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc            >  2:891480/1‑100 (MQ=255)
                                                                                         gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca           <  1:522802/99‑1 (MQ=255)
                                                                                           gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACCGAATGTCAAAAGGCACAg        >  1:2844533/1‑100 (MQ=255)
                                                                                             gcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTa      >  2:2170317/1‑100 (MQ=255)
                                                                                             gcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTa      >  2:2165614/1‑100 (MQ=255)
                                                                                              cTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTAAAAAGGAACAGTaa     >  2:2717725/1‑100 (MQ=255)
                                                                                               taaaTGAGTGCCCGGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaaa    >  2:1963005/3‑100 (MQ=255)
                                                                                               tCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaaa    <  2:2500742/100‑1 (MQ=255)
                                                                                                cAATGAGTGCCCTGTCTCTTAATTAATACGAAACACATATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACCGAATGTCAAAAGACACAGTAAAg   >  1:2728441/1‑100 (MQ=255)
                                                                                                 aaTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGa  <  1:522240/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
CGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGA  >  NC_000913/4227457‑4227653

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