Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,560 T→C 16.4% 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.4% 317.4 / 50.4 145intergenic (+49/‑171)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base T (56/65);  new base C (11/13);  total (67/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 147 total aligned reads displayed.

ACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAG  >  NC_000913/4227461‑4227659
                                                                                                   |                                                                                                   
aCCGAAGTTGAGGGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGt                                                                                                     >  2:1916713/1‑100 (MQ=255)
 ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                                    <  1:2602060/100‑1 (MQ=255)
   gAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGcc                                                                                                  >  2:496501/1‑100 (MQ=255)
      gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                               <  2:378988/100‑1 (MQ=255)
      gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTg                                                                                               <  2:1836812/100‑1 (MQ=255)
        tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt                                                                                             <  1:1575410/100‑2 (MQ=255)
         gaccgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt                                                                                             <  1:1328837/96‑2 (MQ=38)
         gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGttt                                                                                            <  2:698542/100‑3 (MQ=38)
          aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttg                                                                                           <  2:2670244/100‑4 (MQ=38)
           gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACTAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                          <  2:2566401/100‑5 (MQ=38)
           gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                          <  1:1399482/100‑5 (MQ=38)
            cgcTGGCTCGCACCGAATCTGGGGTATGACGCGGACTGATGCACAACTCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                          <  1:1076348/99‑5 (MQ=37)
                ctctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                                     <  1:1729037/98‑1 (MQ=255)
                  ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttattag                                                                                   >  2:1729045/1‑99 (MQ=255)
                  ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                                   >  2:1029169/1‑100 (MQ=255)
                    ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTTACTTTTCCTTACAACAAACAGGGCGCGCAATGAGTGCCCTGTCTCTTTATTaa                                                                                   >  1:207506/1‑98 (MQ=255)
                    ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAata                                                                                 >  2:2443637/1‑100 (MQ=255)
                        ccGAATCAGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaa                                                                             >  2:1532633/1‑100 (MQ=255)
                         cGAATCTGGTGTATGACGCGGACTTATTCACAAATCTGTCATTTTTCCTTACAACAACCAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGa                                                                              <  1:2451497/98‑1 (MQ=255)
                          gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                           >  2:2140554/1‑100 (MQ=255)
                          gAATCTGGGGCATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAACGAGTGCCCTGTCTCTTTATTAATATGAAac                                                                           <  1:1916704/100‑1 (MQ=255)
                           aaTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAaca                                                                          <  1:2618985/100‑1 (MQ=255)
                           aaTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAaca                                                                          >  2:69988/1‑100 (MQ=255)
                             tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCGTGTCTCTTTATTAATATGAAACACt                                                                        >  2:210321/1‑100 (MQ=255)
                             tCTGCGGTATCACGCGGACTGATTCACAAATCTGCCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTCTTAATATg                                                                               <  1:1256148/93‑1 (MQ=255)
                              cTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACtt                                                                       >  1:424080/1‑100 (MQ=255)
                               tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                                      <  1:2797034/100‑1 (MQ=255)
                                ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTtat                                                                     <  1:934012/100‑1 (MQ=255)
                                  ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                                   >  2:104137/1‑100 (MQ=255)
                                   gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                                  <  2:2596082/100‑1 (MQ=255)
                                   gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                                  <  2:2189277/100‑1 (MQ=255)
                                    tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                                 <  1:2281954/100‑1 (MQ=255)
                                      tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                               >  2:1116238/1‑100 (MQ=255)
                                      tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                               <  2:1517551/100‑1 (MQ=255)
                                       gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                              >  2:316065/1‑100 (MQ=255)
                                       gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                              >  2:1509775/1‑100 (MQ=255)
                                         ctcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                            >  1:2027998/3‑100 (MQ=255)
                                         cgcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                            >  1:2737177/1‑100 (MQ=255)
                                            ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                                         <  2:2092767/100‑1 (MQ=255)
                                             gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                           <  1:372795/97‑1 (MQ=255)
                                               cTGATTCACAAATCTGTCACTTTTCCTTACACCAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGcn                                                       <  2:2323127/99‑2 (MQ=255)
                                                tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                                     <  1:2151285/100‑1 (MQ=255)
                                                 gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                                    >  2:2406031/1‑100 (MQ=255)
                                                  aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                                   <  1:1933427/100‑1 (MQ=255)
                                                   ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                                  >  2:230916/1‑100 (MQ=255)
                                                   ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                                  >  2:1505310/1‑100 (MQ=255)
                                                   ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                                  >  2:2807400/1‑100 (MQ=255)
                                                    tCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                                 >  1:2541766/1‑100 (MQ=255)
                                                     cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                                >  1:247281/1‑100 (MQ=255)
                                                       caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                              >  2:616315/1‑100 (MQ=255)
                                                       caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                              <  2:2650674/100‑1 (MQ=255)
                                                          aTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                           <  1:1564512/100‑1 (MQ=255)
                                                           tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                                          <  1:1023068/100‑1 (MQ=255)
                                                           tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                                          <  1:1654412/100‑1 (MQ=255)
                                                            cTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                                         <  2:121795/100‑1 (MQ=255)
                                                             tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGcc                                        <  2:2184537/100‑1 (MQ=255)
                                                             tGTCACTTTTCCTTACAACAAACAGGACGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGcc                                        >  2:1775073/1‑100 (MQ=255)
                                                               tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                                      <  2:2500947/100‑1 (MQ=255)
                                                                cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                                     >  2:821672/1‑100 (MQ=255)
                                                                 aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                                    >  1:2488814/1‑100 (MQ=255)
                                                                   ttttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTTGCCCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                                     <  1:303704/97‑1 (MQ=255)
                                                                    tttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                                 >  1:634165/1‑100 (MQ=255)
                                                                      tCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                               >  1:2055207/1‑100 (MQ=255)
                                                                      tCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                               <  2:1909595/100‑1 (MQ=255)
                                                                        cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                             >  1:2363867/1‑100 (MQ=255)
                                                                        cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                             <  1:976930/100‑1 (MQ=255)
                                                                         ttACAACAACCACGGCGCCCAATGCGTCCCCTGCCCCACCATTAACATGAAACACTTATACTGGAAACAGGCTCGAATCATCTTAGCCCGAAAAACGAAt                            <  1:858275/100‑1 (MQ=255)
                                                                         ttACAACAAACAGGGCGCTCAATGAGTGGCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                            <  1:900167/100‑1 (MQ=255)
                                                                          tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                            <  1:1832063/99‑1 (MQ=255)
                                                                          tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                           >  1:1296103/1‑100 (MQ=255)
                                                                            caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATGCTGGAAACAGGCTGGAATAATCTTAGCAGGAAAAACGAATGTc                         >  1:1868233/1‑100 (MQ=255)
                                                                            caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTAGTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                            <  1:1709719/97‑1 (MQ=255)
                                                                            aaagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAGACAGGATGGAATAATCTTAGCCGGAAAAACGAATg                           >  1:1844787/7‑98 (MQ=37)
                                                                            aaagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                         >  2:1717239/7‑100 (MQ=37)
                                                                             aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa                        <  1:1735208/100‑1 (MQ=255)
                                                                              acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAGACGAATGTCaa                       >  2:532516/1‑100 (MQ=255)
                                                                               caaccaCGGCGCTCAACGAACGCCCTGTCTCTTTATAAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAACCGAATGTCaaa                      <  2:2475009/100‑1 (MQ=255)
                                                                               caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa                      <  2:2415416/100‑1 (MQ=255)
                                                                                 gacaGGGCACTCATTGCGCGCCCTGTCTCTTTATTTATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACAAATGTCAAAAg                    >  1:508241/2‑100 (MQ=255)
                                                                                 aacaGGGCGCTCAATGAGTGCCCTGTCTCTTTAATAATATGAAACACTGATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg                    <  2:1331091/100‑1 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  2:2732471/1‑100 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  2:946850/1‑100 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  2:1286218/1‑100 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  1:702/1‑100 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  2:1835634/1‑100 (MQ=255)
                                                                                     gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac                >  1:842903/1‑100 (MQ=255)
                                                                                     gggCGCTCAATGAGTGCCCTGTCTCTTTAATAATATGAAACACTTATAATGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac                <  1:1269251/100‑1 (MQ=255)
                                                                                      ggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca               <  2:2112090/100‑1 (MQ=255)
                                                                                        cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt             >  2:2056335/1‑100 (MQ=255)
                                                                                        cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt             >  2:1742700/1‑100 (MQ=255)
                                                                                          cTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaa           >  2:972329/1‑100 (MQ=255)
                                                                                            cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAg         >  2:2592477/1‑100 (MQ=255)
                                                                                            cAATGAGTGCCCTGTCTCTGTATTAATATGAAACAATTATACTGGCAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac                <  1:1302213/93‑1 (MQ=255)
                                                                                              aTGAGTGTCCTGTCGCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAATCGAATGTCAAAAGGCACAGTAAAg         <  1:985447/98‑1 (MQ=255)
                                                                                              aTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGAt       <  1:2031165/100‑1 (MQ=255)
                                                                                               tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg      <  1:1087633/100‑1 (MQ=255)
                                                                                               tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg      <  2:896499/100‑1 (MQ=255)
                                                                                                 aGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGcc    <  2:2737195/100‑1 (MQ=255)
                                                                                                  gTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGCCAAAAGGCACCGTAAAGATGCCa   >  2:2666011/1‑100 (MQ=255)
                                                                                                   tGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAg  >  2:1231942/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAG  >  NC_000913/4227461‑4227659

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: