Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,550 G→A 18.5% intergenic (+39/‑181) 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,5500GA18.5% 261.8 / 46.7 118intergenic (+39/‑181)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (56/40);  new base A (8/14);  total (64/54)
Fisher's exact test for biased strand distribution p-value = 9.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 124 total aligned reads displayed.

TATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTA  >  NC_000913/4227451‑4227649
                                                                                                   |                                                                                                   
tATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcg                                                                                                     <  1:314249/100‑1 (MQ=255)
 aTGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGTGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc                                                                                                    <  1:2412702/100‑1 (MQ=255)
  tGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCt                                                                                                   <  2:457994/100‑1 (MQ=255)
   gAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTc                                                                                                  >  2:385820/1‑100 (MQ=255)
      cGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                               >  2:742226/1‑100 (MQ=255)
        ttACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa                                                                                             <  1:674205/100‑1 (MQ=255)
         tACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg                                                                                            >  2:2060961/1‑100 (MQ=255)
           ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                          >  2:2408917/1‑100 (MQ=255)
           ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                          >  1:1102257/1‑100 (MQ=255)
             gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGcc                                                                                        <  1:1037993/100‑1 (MQ=255)
              aaGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGTACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                       <  1:390730/100‑1 (MQ=255)
                gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                     >  2:1377388/1‑100 (MQ=255)
                gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                     <  1:165145/100‑1 (MQ=255)
                gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                     >  2:404842/1‑100 (MQ=255)
                gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                     >  1:2807779/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)
                       gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                              <  2:2664184/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)
                                          ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                            >  2:1811542/1‑99 (MQ=255)
                                            ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                         >  2:2253018/1‑100 (MQ=255)
                                               aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                      <  1:464833/100‑1 (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                                    >  2:2485118/1‑100 (MQ=255)
                                                                 caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                    <  1:88565/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:971262/100‑1 (MQ=255)
                                                                         tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                            >  1:1461250/1‑100 (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)
                                                                                  cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                   >  2:2767506/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)
                                                                                     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)
                                                                                      caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc               >  1:1325998/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:1314504/2‑100 (MQ=39)
                                                                                           gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg          >  2:2566065/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)
                                                                                               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)
                                                                                                   |                                                                                                   
TATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTA  >  NC_000913/4227451‑4227649

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: