Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,560 T→C 16.1% 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.1% 426.8 / 72.5 189intergenic (+49/‑171)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base T (79/78);  new base C (12/18);  total (93/96)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-02

Only 100 of 192 total aligned reads displayed.

ACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAG  >  NC_000913/4227461‑4227659
                                                                                                   |                                                                                                   
aCCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAGCAAACAGGGCGCTCAATGAGt                                                                                                     <  2:3888359/100‑1 (MQ=255)
aCCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGt                                                                                                     >  1:314002/1‑100 (MQ=255)
 ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                                    >  2:2659110/1‑100 (MQ=255)
     aGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCATTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCt                                                                                                <  2:2586860/100‑1 (MQ=255)
     aGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCt                                                                                                <  1:2737495/100‑1 (MQ=255)
      gTTGAGCGCTGGCTGGCATCGAATCTGGATTATGACGAAGACTGATTCACAAATCTGTCACTTTTCGTTACAACAATCAGGCCGCTCAATGAGTGCCCt                                                                                                <  2:634977/99‑1 (MQ=255)
       ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                               <  2:3699183/99‑1 (MQ=255)
       ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                              <  1:2345112/100‑1 (MQ=255)
        tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                             >  1:86086/1‑100 (MQ=255)
        tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                             >  2:2559624/1‑100 (MQ=255)
        tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                             >  2:597179/1‑100 (MQ=255)
        tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt                                                                                             <  2:3248971/100‑2 (MQ=39)
         gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCAAAAATCTGTCACTTTTCCTTGCAACAAACAGGGCGCTCAATGAGTGCCCTGt                                                                                              >  2:868883/1‑98 (MQ=255)
          aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtct                                                                                            <  2:1204923/99‑1 (MQ=255)
           gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                          <  1:3509421/100‑5 (MQ=38)
            cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgtt                                                                                         <  2:2513425/100‑6 (MQ=38)
             gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                                        <  1:1334850/100‑7 (MQ=37)
              ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                                       >  1:3089427/1‑100 (MQ=255)
               tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttat                                                                                      <  1:138307/100‑1 (MQ=255)
                ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                                     <  1:7457/100‑1 (MQ=255)
                ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTctatt                                                                                     <  2:1657282/100‑1 (MQ=255)
                  ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                                   >  1:1797794/1‑100 (MQ=255)
                  ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                                   >  2:1692826/1‑100 (MQ=255)
                     gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAatat                                                                                >  1:2357530/1‑100 (MQ=255)
                     gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAatat                                                                                >  1:3049730/1‑100 (MQ=255)
                      cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                               >  2:4098570/1‑100 (MQ=255)
                        ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaa                                                                             <  2:1573933/100‑1 (MQ=255)
                        ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaa                                                                             >  1:3385900/1‑100 (MQ=255)
                         cGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaaa                                                                            >  1:158737/1‑100 (MQ=255)
                         cGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaaa                                                                            >  2:1067950/1‑100 (MQ=255)
                              cTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACtt                                                                       >  2:687606/1‑100 (MQ=255)
                               tGGTGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                                      <  2:1991903/100‑1 (MQ=255)
                                ggggTATTACGCGGACTGATTCACAAATCTGTCATTTTTCCTTACAATAAACAGGGGGCTCAATGAGTGCAGTGTCTCTTTATTAATATGAAACACTtag                                                                     <  2:174497/100‑2 (MQ=255)
                                  ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                                   <  2:2134227/100‑1 (MQ=255)
                                  ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                                   >  1:4008871/1‑100 (MQ=255)
                                  ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCATTATTAATATGAAACACTTATAc                                                                   >  1:1133997/1‑100 (MQ=255)
                                   gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                                  <  2:525864/100‑1 (MQ=255)
                                   gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                                  >  1:2382093/1‑100 (MQ=255)
                                    tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                                 >  1:1005822/1‑100 (MQ=255)
                                    tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                                 >  2:3041008/1‑100 (MQ=255)
                                    tATGACGCAGACTTATTCACAAATCTGTCACTTTTTCTTACAACAAACAGTGCGCTCAATGATTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                                  <  2:3485902/99‑1 (MQ=255)
                                      tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                               <  1:2190309/100‑1 (MQ=255)
                                       gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                              <  2:682685/100‑1 (MQ=255)
                                       gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                              >  2:4046287/1‑100 (MQ=255)
                                       gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                              <  1:3181943/100‑1 (MQ=255)
                                         cgcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                            <  1:1966253/100‑1 (MQ=255)
                                            ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                                         <  1:1606643/100‑1 (MQ=255)
                                            ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                                         >  1:2923457/1‑100 (MQ=255)
                                             gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                                        <  2:2402419/100‑1 (MQ=255)
                                                tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                                     <  1:529813/100‑1 (MQ=255)
                                                tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                                     <  2:592108/100‑1 (MQ=255)
                                                tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                                     <  2:2727648/100‑1 (MQ=255)
                                                tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                                     <  1:3196643/100‑1 (MQ=255)
                                                 gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                                    >  2:2287889/1‑100 (MQ=255)
                                                  aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                                   <  2:393664/100‑1 (MQ=255)
                                                  aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                                   <  1:838429/100‑1 (MQ=255)
                                                   ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATCTGAAACACTTATACTGGAAACAGGCTGGaat                                                  >  1:444165/1‑100 (MQ=255)
                                                   ttCACAAATCTGTAACTTTTCCTTACAACAAACAGGGCGATAAATGAGTGCACTGTCTCTTTATTAATCTGAAACACTTATACTGGAAACAAGCTAGa                                                    >  2:2933249/1‑98 (MQ=255)
                                                    tCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                                 >  2:747351/1‑100 (MQ=255)
                                                     cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                                >  2:1847164/1‑100 (MQ=255)
                                                     cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                                >  2:1001751/1‑100 (MQ=255)
                                                       caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                              >  1:1702016/1‑100 (MQ=255)
                                                       caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                              <  1:2990656/100‑1 (MQ=255)
                                                         aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                            >  1:3134775/1‑100 (MQ=255)
                                                             tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGcc                                        <  2:3004508/100‑1 (MQ=255)
                                                             tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGcc                                        >  2:283626/1‑100 (MQ=255)
                                                               tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGGGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                                      >  1:2962855/1‑100 (MQ=255)
                                                                cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                                     >  1:3279341/1‑100 (MQ=255)
                                                                       ccTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                              >  1:761026/1‑100 (MQ=255)
                                                                       ccTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                              >  2:1197782/1‑100 (MQ=255)
                                                                        cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                             >  1:3432343/1‑100 (MQ=255)
                                                                        cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                             <  2:2038070/100‑1 (MQ=255)
                                                                        cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                             >  2:3159385/1‑100 (MQ=255)
                                                                         ttACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                            <  1:217617/100‑1 (MQ=255)
                                                                          tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                           <  2:2991842/100‑1 (MQ=255)
                                                                          ataaagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                           >  1:2257633/9‑100 (MQ=37)
                                                                          ataaagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACACGCTGCAATAATCTAACCCAAAAAAACaaa                             >  2:3498713/9‑98 (MQ=16)
                                                                           acaacaAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt                          >  2:3926381/1‑100 (MQ=255)
                                                                            caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc                         >  2:1383535/1‑100 (MQ=255)
                                                                              acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGACGGAAAAACGAATGTCaa                       >  2:1868436/1‑100 (MQ=255)
                                                                               caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa                      >  1:3643596/1‑100 (MQ=255)
                                                                               caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa                      >  1:3964977/1‑100 (MQ=255)
                                                                                 gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg                    >  1:1370122/2‑100 (MQ=39)
                                                                                 gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGAAAAAACGAATGTCaaaa                     >  2:3391636/2‑99 (MQ=38)
                                                                                 aacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg                    <  1:937183/100‑1 (MQ=255)
                                                                                  acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg                   <  2:3134761/100‑1 (MQ=255)
                                                                                   caGGGCGCTCAATGCGTGCCCTGTCTCTTTATTAATATGCAACACTTATACTGGAACCAGGCTGGAATAATCTTATCCGGAAAAACGAATGTCTAAAGGc                  >  2:1921837/1‑100 (MQ=255)
                                                                                   caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  1:378476/1‑100 (MQ=255)
                                                                                   caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc                  >  2:1772776/1‑100 (MQ=255)
                                                                                      ggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca               >  2:2964539/1‑100 (MQ=255)
                                                                                      ggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac                >  1:3989769/1‑99 (MQ=255)
                                                                                        cgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt             >  1:1225944/1‑100 (MQ=255)
                                                                                         gcTCAATCAGTCCCCGGTCTCTTTATTAATATGAAACACTTAGACTGCAAACAGGCTGCAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTa            <  2:3708885/100‑1 (MQ=255)
                                                                                          cacAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCAGTAAAAACGAATGTCAAAAGGCACAGTaa           >  2:727067/3‑100 (MQ=255)
                                                                                           tCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaaa          <  1:598962/100‑1 (MQ=255)
                                                                                           tCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTaaa          <  1:2140058/100‑1 (MQ=255)
                                                                                            cAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAg         >  1:4022658/1‑100 (MQ=255)
                                                                                              aTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGATGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGAt       <  2:1342394/100‑1 (MQ=255)
                                                                                               tGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATg      <  1:507066/100‑1 (MQ=255)
                                                                                                gAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGc     <  1:402176/100‑1 (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: