Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,225,139 T→C intergenic (+49/‑171) metH → / → yjbB B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase/putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,225,1390TC98.5% 614.5 / NA 208intergenic (+49/‑171)metH/yjbBB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase/putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base C (105/97);  minor base G (3/0);  total (111/97)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

GCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGT  >  NC_000913/4225051‑4225227
                                                                                        |                                                                                        
  gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTg                                                                                    >  1:2248352/1‑93 (MQ=255)
    tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                  >  2:6346998/1‑93 (MQ=255)
    tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                  <  2:3842481/93‑1 (MQ=255)
    tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                  <  2:609541/93‑1 (MQ=255)
     ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtct                                                                                 >  2:6698116/1‑93 (MQ=255)
          gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTtta                                                                            <  1:5177214/93‑1 (MQ=255)
          gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTtta                                                                            <  1:2503491/93‑1 (MQ=255)
           cACCGAATCTGGGTTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttat                                                                           <  1:1287193/93‑1 (MQ=255)
           cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttat                                                                           >  2:6401926/1‑93 (MQ=255)
            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttatt                                                                          <  1:914555/93‑1 (MQ=255)
            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCCGTCTCTttatt                                                                          >  2:4768060/1‑93 (MQ=255)
            aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGCGCGCTCAATGAGCGCCCTGTCTCTttatt                                                                          >  1:6724886/1‑93 (MQ=255)
             ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttatta                                                                         >  2:6022098/1‑93 (MQ=255)
             ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttatta                                                                         <  1:6708746/93‑1 (MQ=255)
               gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCCCTTTTCCTTACAACAAACAGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAat                                                                       <  1:2467434/93‑1 (MQ=255)
               gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAat                                                                       >  2:3420907/1‑93 (MQ=255)
                 aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAGGGCCCTGTCTCTTTTTTTatat                                                                     >  1:5595396/1‑93 (MQ=255)
                  tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAGCGCCCCGTCTCTTTTTTTatatt                                                                    >  1:2888636/1‑92 (MQ=255)
                  tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATg                                                                    >  1:2048460/1‑93 (MQ=255)
                   cTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGa                                                                   >  1:4659684/1‑93 (MQ=255)
                    tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCATTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                  <  2:1235812/93‑1 (MQ=255)
                     ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAACCAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaaa                                                                 <  2:2510444/93‑1 (MQ=255)
                      gggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAac                                                                <  1:4121732/93‑1 (MQ=255)
                      gggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAac                                                                <  1:6131674/93‑1 (MQ=255)
                        gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAacac                                                              >  2:4458299/1‑93 (MQ=255)
                           tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta                                                           <  2:5221324/93‑1 (MQ=255)
                           tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta                                                           >  2:2635455/1‑93 (MQ=255)
                           tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGATCGATTATCCGCCTACTTCTTTTTTAATATTAAACACttt                                                           >  1:7193128/1‑92 (MQ=255)
                           tGACGCGGACTGATTCACAAATCTGTCACTCTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta                                                           <  2:7731163/93‑1 (MQ=255)
                             aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTtata                                                         <  1:2654569/93‑1 (MQ=255)
                             aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTtata                                                         <  2:3787524/93‑1 (MQ=255)
                             aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGATCAATAAGCGCCCTGTCTCTTTATTAATATGAAACACTtata                                                         <  1:4944504/93‑1 (MQ=255)
                               gcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                       >  2:3883726/1‑93 (MQ=255)
                                cgGACTGATTCACAAACCTGTCACTTTTCCTTACAACAAACAGGGCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                      <  2:6532750/93‑1 (MQ=255)
                                  gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                    >  1:7499019/1‑93 (MQ=255)
                                  gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                    >  1:6529968/1‑93 (MQ=255)
                                  gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                    >  1:6188471/1‑93 (MQ=255)
                                  gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                    >  1:542557/1‑93 (MQ=255)
                                  gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                    >  1:4596210/1‑93 (MQ=255)
                                   aCTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                   >  2:5488445/1‑93 (MQ=255)
                                      gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                <  1:7632878/93‑1 (MQ=255)
                                      gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                <  2:4065385/93‑1 (MQ=255)
                                      gATTCACAAATCTGTCACTTTTCCTAACAACAAACGGGGCGCTCAATGGGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                <  1:2140840/93‑1 (MQ=255)
                                       aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                               <  2:7746336/93‑1 (MQ=255)
                                       aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                               >  1:4908045/1‑93 (MQ=255)
                                       aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                               <  2:4984439/93‑1 (MQ=255)
                                        ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTTATATGAAACACTTATACTGGAAACAgg                                              >  1:7143599/1‑93 (MQ=255)
                                        ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                              <  1:3668784/93‑1 (MQ=255)
                                         tCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                             <  2:1040419/93‑1 (MQ=255)
                                            caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                          >  2:428676/1‑93 (MQ=255)
                                             aaaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                         >  2:5359030/1‑93 (MQ=255)
                                               aTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                       >  1:4054261/1‑93 (MQ=255)
                                               aTCTCTCACTTCTGCTAAAACCAAGCAGGGCCCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                       <  1:3479222/93‑1 (MQ=255)
                                                 cTGTCTTTTTTCCTAACAAAAAACGGGGCGCTCAAGGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                     <  2:7162620/93‑1 (MQ=255)
                                                 cTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                     <  2:7283908/93‑1 (MQ=255)
                                                   gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                   >  1:1101997/1‑93 (MQ=255)
                                                    tcttttttCCTAACAAAAAGCGGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                  <  2:1699379/89‑1 (MQ=255)
                                                    tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                  <  2:4809095/93‑1 (MQ=255)
                                                    tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                  <  2:6841361/93‑1 (MQ=255)
                                                    tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                  <  1:7731875/93‑1 (MQ=255)
                                                     ccttttcCCTTACAAAAACCGGGGCGCTGAAGGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                 <  2:5836245/90‑1 (MQ=255)
                                                     cACTTTCCCTTACAACAAACAGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                 <  2:2402065/93‑1 (MQ=255)
                                                      aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                >  1:3245309/1‑93 (MQ=255)
                                                      aCTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                <  2:2684442/93‑1 (MQ=255)
                                                       tttttCCACACAACAAGCAGGGCCCTAAATGGCCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                               <  1:6940003/92‑1 (MQ=255)
                                                       cTTTTCCTTACAACAAACAGGGCGCTAAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                               <  2:7424147/93‑1 (MQ=255)
                                                        ttttCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                              >  1:3494163/1‑93 (MQ=255)
                                                          ttCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                            >  2:894201/1‑93 (MQ=255)
                                                           tCCTTACAACAAACAGGGCGCTCAATGAGCGCCCCGTCTCTTTTTTTATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                           >  2:1787443/1‑93 (MQ=255)
                                                            ccTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACCGGCTGGAATAATCTTAGCCGGa                          >  2:2569314/1‑93 (MQ=255)
                                                            ccTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                          >  1:2069910/1‑93 (MQ=255)
                                                            ccTTACAAAAAACAGGCCCGCAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                          <  2:3402703/93‑1 (MQ=255)
                                                              ttACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                        >  1:659554/1‑93 (MQ=255)
                                                               taaaacaaacaGGGCGCTAAAGGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaa                       <  1:1477907/90‑1 (MQ=255)
                                                                acaacaAACAGGTCTTCCAATTAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGTCTTGAATAATCTTAGCCGGaaaaa                      >  2:7237979/1‑93 (MQ=255)
                                                                acaacaAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                      <  1:6847367/93‑1 (MQ=255)
                                                                acaacaAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                      <  1:1148214/93‑1 (MQ=255)
                                                                acaaaaAACGGGCCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                      <  2:2863999/93‑1 (MQ=255)
                                                                  aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                    >  1:1006470/1‑93 (MQ=255)
                                                                   acaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCTTAGCCGGAAAAAAGa                   >  2:1958464/1‑93 (MQ=255)
                                                                    caaacaGGGCTCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                  <  2:2334789/93‑1 (MQ=255)
                                                                    caaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                  <  1:4819770/93‑1 (MQ=255)
                                                                    caaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                  <  1:3741992/93‑1 (MQ=255)
                                                                      aacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                <  1:1266650/93‑1 (MQ=255)
                                                                       acaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt               >  2:663417/1‑93 (MQ=255)
                                                                        caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc              >  2:1196270/1‑93 (MQ=255)
                                                                        caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc              <  2:7899116/93‑1 (MQ=255)
                                                                          gggCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa            >  1:2730259/1‑93 (MQ=255)
                                                                           ggCGCTCAATGAGGGCCCTGTCTCTTTTTTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCCTTTCCGGGAAAACGAAAGTCaaa           >  2:1877978/1‑93 (MQ=255)
                                                                            ggtcTCAAGGGGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa          <  2:3245308/90‑1 (MQ=255)
                                                                            gcgcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa          <  2:5035075/93‑1 (MQ=255)
                                                                            gcgcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa          <  1:586327/93‑1 (MQ=255)
                                                                            gcgcTCAATGACCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa          <  2:4548713/93‑1 (MQ=255)
                                                                            gcgcTCAACGACCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa          <  2:397698/93‑1 (MQ=255)
                                                                             ctcTCAACGAGCCCGCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg         <  1:4819545/91‑1 (MQ=255)
                                                                               cTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc       >  2:5470397/1‑93 (MQ=255)
                                                                                 cAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac     >  2:1299946/1‑93 (MQ=255)
                                                                                 cAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac     >  2:5785297/1‑93 (MQ=255)
                                                                                    tGAGCGCCCTGTCTCTTTATTAATATGAAACACTtata                                                         <  1:7793094/38‑1 (MQ=255)
                                                                                    tGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAga  >  2:6432125/1‑92 (MQ=255)
                                                                                        |                                                                                        
GCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGT  >  NC_000913/4225051‑4225227

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