Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,560 T→C 16.9% intergenic (+49/‑171) 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,5600TC16.9% 268.2 / 48.3 124intergenic (+49/‑171)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base T (61/42);  new base C (8/13);  total (69/55)
Fisher's exact test for biased strand distribution p-value = 9.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

Only 100 of 126 total aligned reads displayed.

ACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCA  >  NC_000913/4227461‑4227658
                                                                                                   |                                                                                                  
aCCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGt                                                                                                    >  2:111065/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                                                                                              >  2:404842/1‑100 (MQ=255)
      gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                              >  2:1377388/1‑100 (MQ=255)
      gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                              >  1:2807779/1‑100 (MQ=255)
      gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                              <  1:659400/100‑1 (MQ=255)
      gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                              <  1:1653966/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:1569080/100‑3 (MQ=38)
         gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGttt                                                                                           <  1:622194/100‑3 (MQ=38)
          aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtctc                                                                                          >  1:2811022/1‑100 (MQ=255)
           gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTctct                                                                                         <  2:2567753/100‑1 (MQ=255)
           gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTctct                                                                                         >  1:118882/1‑100 (MQ=255)
            cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgtt                                                                                        <  1:2173696/100‑6 (MQ=38)
             gctggctgGCACCGAGTGTCGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTGCCTTACAACAAACAGGGCGCTGAATGAGTGCCCTGGCTCgtt                                                                                       <  1:1203242/100‑1 (MQ=255)
             gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                                       <  2:2664184/100‑7 (MQ=37)
             gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                                       <  2:1189494/100‑7 (MQ=37)
              ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                      <  1:1058346/100‑1 (MQ=255)
              ctggctggCACCGAATCTGGCGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgt                                                                                      <  1:2298535/100‑8 (MQ=37)
                 gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgtaag                                                                                   <  2:694713/100‑11 (MQ=25)
                  ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                                  >  2:677990/1‑100 (MQ=255)
                   tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAat                                                                                 <  1:2340789/100‑1 (MQ=255)
                    ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                                <  2:1986996/100‑1 (MQ=255)
                     gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAatat                                                                               >  1:1167103/1‑100 (MQ=255)
                      cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                              <  1:487768/100‑1 (MQ=255)
                        ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATgn                                                                             <  1:2612033/99‑2 (MQ=255)
                          gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                          >  1:262277/1‑100 (MQ=255)
                                ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                                     >  2:1811542/1‑99 (MQ=255)
                                  ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                                  >  2:2253018/1‑100 (MQ=255)
                                   gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                                 >  1:439181/1‑100 (MQ=255)
                                     aTGACGCGGACTGATTGACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                               <  1:1696863/100‑1 (MQ=255)
                                      tGACGCGCACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                              >  2:1532650/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)
                                                tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                                    <  2:2366642/100‑1 (MQ=255)
                                                tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                                    <  2:383336/100‑1 (MQ=255)
                                                 gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                                   >  2:2496250/1‑100 (MQ=255)
                                                   ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                                 >  2:1281675/1‑100 (MQ=255)
                                                     cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                               >  1:2353898/1‑100 (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)
                                                         aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                           <  1:389219/100‑1 (MQ=255)
                                                          aTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                          <  1:1499351/100‑1 (MQ=255)
                                                          aTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                          <  2:20452/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:2153568/1‑100 (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)
                                                                 aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                                   >  1:1050320/1‑100 (MQ=255)
                                                                  cTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                                  >  2:111482/1‑100 (MQ=255)
                                                                     ttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                               <  2:2743132/100‑1 (MQ=255)
                                                                      tCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                              <  2:1167105/100‑1 (MQ=255)
                                                                        cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                            >  2:1595897/1‑100 (MQ=255)
                                                                        cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                            >  2:2767506/1‑100 (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:1890462/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)
                                                                                   caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAAAACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                 <  1:2452640/100‑1 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                 >  2:891480/1‑100 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                 >  2:1992789/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)
                                                                                           tCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaaa         <  2:2500742/100‑1 (MQ=255)
                                                                                            cAATGAGTGCCCTGTCTCTTAATTAATACGAAACACATATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACCGAATGTCAAAAGACACAGTAAAg        >  1:2728441/1‑100 (MQ=255)
                                                                                             ccgatctGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGa       <  1:1244281/94‑1 (MQ=255)
                                                                                             aaTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGa       <  1:522240/100‑1 (MQ=255)
                                                                                                 tgTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc   <  1:721143/99‑1 (MQ=255)
                                                                                                 aGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc   >  2:2556409/1‑100 (MQ=255)
                                                                                                 aGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc   <  1:2207786/100‑1 (MQ=255)
                                                                                                  gTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCa  >  2:556744/1‑100 (MQ=255)
                                                                                                  gTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCa  >  2:721145/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
ACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCA  >  NC_000913/4227461‑4227658

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