Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,555 A→T 16.8% 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,5550AT16.8% 414.4 / 70.9 171intergenic (+44/‑176)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base A (70/72);  new base T (15/14);  total (85/86)
Fisher's exact test for biased strand distribution p-value = 8.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

Only 100 of 175 total aligned reads displayed.

GCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGAT  >  NC_000913/4227456‑4227654
                                                                                                   |                                                                                                   
 cGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                                    <  1:146517/100‑1 (MQ=255)
  gTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg                                                                                                   <  2:581448/100‑1 (MQ=255)
    cgccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg                                                                                                   <  1:1537056/96‑1 (MQ=255)
      ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgat                                                                                               >  2:1537065/1‑98 (MQ=255)
      ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGtn                                                                                               <  1:1971288/100‑2 (MQ=255)
         aaGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                            <  2:3476691/100‑1 (MQ=255)
         aaGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGccc                                                                                            <  2:322332/100‑1 (MQ=255)
           ggtGAGCGCTGGCTGGCACCGAATCTGTGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                          <  1:3149529/98‑1 (MQ=255)
           gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                          >  2:2029568/1‑100 (MQ=255)
           gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                          >  2:1489796/1‑100 (MQ=255)
           gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                          <  2:705569/100‑1 (MQ=255)
             tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt                                                                                        <  2:1807456/100‑2 (MQ=39)
              gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGttt                                                                                       <  2:359673/100‑3 (MQ=39)
               aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtctc                                                                                      >  1:757663/1‑100 (MQ=255)
                 cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCtt                                                                                    <  1:2604245/100‑1 (MQ=255)
                 cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCtt                                                                                    <  2:1287626/100‑1 (MQ=255)
                 cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCtt                                                                                    <  1:264266/100‑1 (MQ=255)
                    tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttat                                                                                 <  2:3748500/100‑1 (MQ=255)
                    tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttat                                                                                 <  2:3108814/100‑1 (MQ=255)
                     ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                                >  1:2939691/1‑100 (MQ=255)
                      gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttgtaag                                                                               <  1:819803/100‑11 (MQ=37)
                      gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCGTACATCAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                                <  1:90550/99‑1 (MQ=255)
                       ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                              >  2:2617451/1‑100 (MQ=255)
                         ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTGCTTACGACAAACATGGCGATAAAAGAGTGCCCTTTTGCTttattc                                                                               >  2:2547984/1‑96 (MQ=255)
                         ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                            <  2:3403828/100‑1 (MQ=255)
                         ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                            <  1:676346/100‑1 (MQ=255)
                           cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                          >  1:1416037/1‑100 (MQ=255)
                           cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                          >  2:2466729/1‑100 (MQ=255)
                            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttagatcggaag                                                                         >  1:1538246/1‑91 (MQ=255)
                               gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                      >  1:2212661/1‑100 (MQ=255)
                               gAATCTGGGGTATGACGCCGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                      <  1:1591670/100‑1 (MQ=255)
                                 aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCGCTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacac                                                                    >  1:1039172/1‑100 (MQ=255)
                                  tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                                   <  2:1466086/100‑1 (MQ=255)
                                    tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                                 >  2:1001967/1‑100 (MQ=255)
                                     ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTtat                                                                <  1:3185732/100‑1 (MQ=255)
                                     ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTtat                                                                >  2:2675585/1‑100 (MQ=255)
                                      gggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTtat                                                                <  1:3542371/99‑1 (MQ=255)
                                       ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                              <  2:3607277/100‑1 (MQ=255)
                                       ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                              >  2:1366542/1‑100 (MQ=255)
                                        gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                             >  2:1036759/1‑100 (MQ=255)
                                         tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                            >  2:3657584/1‑100 (MQ=255)
                                         tATGACGCGGACTGACTCACAAATCTGTCACTTTTCCTTACAATAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                            <  1:949746/100‑1 (MQ=255)
                                            gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                         >  2:3095004/1‑100 (MQ=255)
                                             aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                        <  2:3413431/100‑1 (MQ=255)
                                             aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                        >  1:1017081/1‑100 (MQ=255)
                                              cgcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACCAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                       >  1:1622138/1‑100 (MQ=255)
                                               gcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                      >  2:217497/1‑100 (MQ=255)
                                                  gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                                   <  1:3756892/100‑1 (MQ=255)
                                                    cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                                 <  1:1195159/100‑1 (MQ=255)
                                                      gcttCACAACTCTGTCACTTTTCCTTACAACCAACAGGGCGGTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                                <  1:2677954/97‑1 (MQ=255)
                                                       aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                              <  2:977897/100‑1 (MQ=255)
                                                       aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                              >  1:3365787/1‑100 (MQ=255)
                                                        ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                             >  2:3703871/1‑100 (MQ=255)
                                                        ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                             >  2:2473767/1‑100 (MQ=255)
                                                         tCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                            >  2:2317388/1‑100 (MQ=255)
                                                         acaccAATCTCTCACTTTTCCATACAACACACAGGGCCCTCAATGAGTCACCTGTCTATTTATTAATATGAAACACTTATACTGAAAACAGGCt                                                  >  1:552365/2‑94 (MQ=255)
                                                          cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                           >  1:1240935/1‑100 (MQ=255)
                                                          cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                           >  2:2828337/1‑100 (MQ=255)
                                                            caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                         >  2:507427/1‑100 (MQ=255)
                                                            caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                         <  2:1059821/100‑1 (MQ=255)
                                                             aaaTCTGTAACGTTTCCCTACAACAAGCAGGGCGCGCAATGAGTGCACGGTCTCTTTATCAATATGAAACACGTATACTGGGAACAGGCTGGAATAATc                                         <  1:3222402/99‑1 (MQ=255)
                                                              aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                       >  1:2548163/1‑100 (MQ=255)
                                                                tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                                     <  2:1017088/100‑1 (MQ=255)
                                                                  tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGcc                                   <  2:2398580/100‑1 (MQ=255)
                                                                     cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                                <  1:1733129/100‑1 (MQ=255)
                                                                        ctttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                                >  1:2988138/2‑97 (MQ=255)
                                                                           tCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                          <  1:1809286/100‑1 (MQ=255)
                                                                             cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                        <  2:1508329/100‑1 (MQ=255)
                                                                             cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAAAACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                        >  1:882729/1‑100 (MQ=255)
                                                                              gatctagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                       <  2:160648/91‑1 (MQ=37)
                                                                              gatctagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                       <  2:161551/91‑1 (MQ=38)
                                                                               tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                      <  2:623580/100‑1 (MQ=255)
                                                                               tACAACAAACAGGGCGCTCAATGAGTGCCCTGGCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                      <  1:1301720/100‑1 (MQ=255)
                                                                                acaacaAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt                     >  2:3812001/1‑100 (MQ=255)
                                                                                   agagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAtagatc                  >  1:161548/5‑95 (MQ=37)
                                                                                    caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa                 >  1:962039/1‑100 (MQ=255)
                                                                                      gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg               >  2:2804283/2‑100 (MQ=39)
                                                                                      gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg               >  2:19473/2‑100 (MQ=39)
                                                                                      gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg               >  2:559769/2‑100 (MQ=39)
                                                                                       acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg              >  1:1363645/1‑100 (MQ=255)
                                                                                       acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg              >  1:1702481/1‑100 (MQ=255)
                                                                                       acaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg              >  1:1873932/1‑100 (MQ=255)
                                                                                        caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc             >  2:1139469/1‑100 (MQ=255)
                                                                                        caGGGCGCTCAATGAGGGCCCTGTCTCTTTATTAAGAGGAAACACTTATACTGGAAACAGGCAGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg              <  1:1585655/99‑1 (MQ=255)
                                                                                        caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc             >  1:2824178/1‑100 (MQ=255)
                                                                                        caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc             >  2:1252533/1‑100 (MQ=255)
                                                                                        caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc             >  1:2465041/1‑100 (MQ=255)
                                                                                        caGGGCACTCATTAAGCGCCCTGTCTCTTTATTAATATAAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc             >  1:2350915/1‑100 (MQ=255)
                                                                                         aGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca            >  1:3673363/1‑100 (MQ=255)
                                                                                          gggCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac           >  1:3237990/1‑100 (MQ=255)
                                                                                           ggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca          <  2:3463077/100‑1 (MQ=255)
                                                                                           ggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca          <  1:3657561/100‑1 (MQ=255)
                                                                                            gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg         >  1:1810822/1‑100 (MQ=255)
                                                                                             cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg         >  2:2731366/1‑99 (MQ=255)
                                                                                             cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt        >  2:1582875/1‑100 (MQ=255)
                                                                                                 cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAg    >  2:250934/1‑100 (MQ=255)
                                                                                                 cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAg    >  1:661253/1‑100 (MQ=255)
                                                                                                 cAAAGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAACGGCACAGTaaa     <  1:2275902/99‑1 (MQ=255)
                                                                                                   aTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGAt  >  1:3040602/1‑100 (MQ=255)
                                                                                                   aTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGAt  <  1:336730/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGAT  >  NC_000913/4227456‑4227654

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: