Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,130,398 A→G D68G (GAC→GGC)  metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,130,3980AG85.5% 671.6 / 103.8 290D68G (GAC→GGC) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (15/27);  new base G (119/129);  total (134/156)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

Only 100 of 291 total aligned reads displayed.

GAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGC  >  NC_000913/4130307‑4130490
                                                                                           |                                                                                            
gAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGAc                                                                                             <  2:1338046/93‑1 (MQ=255)
gAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGAc                                                                                             <  1:4160112/93‑1 (MQ=255)
   cAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGt                                                                                          <  1:1831288/93‑1 (MQ=255)
     gACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTAc                                                                                        <  1:4701122/93‑1 (MQ=255)
      aCCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACg                                                                                       >  2:4473546/1‑93 (MQ=255)
       cccTGTGGAACTCCATCGATCGCCTTACCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGc                                                                                      >  1:855861/1‑93 (MQ=255)
        ccTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGca                                                                                     >  2:1488317/1‑93 (MQ=255)
          tgtgGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGcaca                                                                                   <  2:1562847/93‑1 (MQ=255)
          tgtgGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGcaca                                                                                   >  1:5451329/1‑93 (MQ=255)
           gtgGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAg                                                                                  >  1:5280743/1‑93 (MQ=255)
           gtgGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAg                                                                                  >  2:6630368/1‑93 (MQ=255)
             ggAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCa                                                                                <  2:5375804/93‑1 (MQ=255)
             ggAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCa                                                                                <  1:2474157/93‑1 (MQ=255)
                aCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCatta                                                                             <  1:2080516/93‑1 (MQ=255)
                aCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCatta                                                                             >  1:3424965/1‑93 (MQ=255)
                 cTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCattat                                                                            <  2:4058741/93‑1 (MQ=255)
                 cTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCattat                                                                            <  1:1736725/93‑1 (MQ=255)
                   ccATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCAttatta                                                                          <  2:1028350/93‑1 (MQ=255)
                     atcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTaaa                                                                        <  2:4974213/93‑1 (MQ=255)
                      tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAg                                                                       >  1:6852381/1‑93 (MQ=255)
                      tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAg                                                                       >  2:3182372/1‑93 (MQ=255)
                      tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAg                                                                       >  2:5098738/1‑93 (MQ=255)
                       cgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAgg                                                                      >  1:5596124/1‑93 (MQ=255)
                        gatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGc                                                                     >  1:1228013/1‑93 (MQ=255)
                         atcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCa                                                                    >  2:7112662/1‑93 (MQ=255)
                         atcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCa                                                                    <  2:7404801/93‑1 (MQ=255)
                          tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCAt                                                                   >  1:6604130/1‑93 (MQ=255)
                          tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCAt                                                                   >  1:5385572/1‑93 (MQ=255)
                          tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCAt                                                                   >  2:396858/1‑93 (MQ=255)
                           cgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCAtt                                                                  <  1:7170757/93‑1 (MQ=255)
                            gCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTa                                                                 <  2:7500037/93‑1 (MQ=255)
                            gCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTa                                                                 >  2:650799/1‑93 (MQ=255)
                            gCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTa                                                                 >  1:4868145/1‑93 (MQ=255)
                            gCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTa                                                                 >  1:3300362/1‑93 (MQ=255)
                             ccTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaa                                                                <  2:2759538/93‑1 (MQ=255)
                             ccTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaa                                                                >  1:6233607/1‑93 (MQ=255)
                             ccTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaa                                                                <  2:1606229/93‑1 (MQ=255)
                              cTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaaa                                                               >  1:3541202/1‑93 (MQ=255)
                              cTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaaa                                                               <  1:955822/93‑1 (MQ=255)
                              cTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTaaa                                                               <  2:5242416/93‑1 (MQ=255)
                               ttAGCAGCCTGAAACCGAAGTTTTTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAg                                                              <  2:6266077/93‑1 (MQ=255)
                               ttAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAg                                                              >  1:4030692/1‑93 (MQ=255)
                                 agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGAt                                                            <  2:4120818/93‑1 (MQ=255)
                                  gcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATc                                                           <  1:459693/93‑1 (MQ=255)
                                  gcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATc                                                           >  1:3132235/1‑93 (MQ=255)
                                  gcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATc                                                           <  2:1187206/93‑1 (MQ=255)
                                   cagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCg                                                          <  2:3432175/93‑1 (MQ=255)
                                    agcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGc                                                         <  2:7333219/93‑1 (MQ=255)
                                    agcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGc                                                         >  1:2203239/1‑93 (MQ=255)
                                     gcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCa                                                        <  1:50406/93‑1 (MQ=255)
                                      ccTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCAc                                                       <  2:5273730/93‑1 (MQ=255)
                                       cTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACt                                                      <  1:6082652/93‑1 (MQ=255)
                                       cTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACt                                                      >  2:7293799/1‑93 (MQ=255)
                                       cTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACt                                                      <  1:1894667/93‑1 (MQ=255)
                                        tGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTg                                                     <  2:5669163/93‑1 (MQ=255)
                                         gAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTgg                                                    >  2:324170/1‑93 (MQ=255)
                                         gAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTgg                                                    <  1:1007639/93‑1 (MQ=255)
                                          aaaCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGt                                                   >  2:5773043/1‑93 (MQ=255)
                                           aaCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTc                                                  >  1:3710643/1‑93 (MQ=255)
                                            aCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGctct                                                 <  1:6969969/93‑1 (MQ=255)
                                              cGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTgg                                               >  2:2641585/1‑93 (MQ=255)
                                               gAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGa                                              >  2:1751396/1‑93 (MQ=255)
                                               gAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGa                                              <  1:3924868/93‑1 (MQ=255)
                                                aaGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGaa                                             <  1:6683529/93‑1 (MQ=255)
                                                 aGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAg                                            <  1:3215882/93‑1 (MQ=255)
                                                  gTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGc                                           <  2:6249976/93‑1 (MQ=255)
                                                  gTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGc                                           <  1:1851936/93‑1 (MQ=255)
                                                    ttGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCgg                                         <  1:7356078/93‑1 (MQ=255)
                                                     tGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGc                                        <  1:1145189/93‑1 (MQ=255)
                                                      gTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCa                                       <  2:550659/93‑1 (MQ=255)
                                                       tATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAc                                      <  1:2148637/93‑1 (MQ=255)
                                                       tATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAc                                      <  2:1386991/93‑1 (MQ=255)
                                                         tCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCg                                    <  1:520870/93‑1 (MQ=255)
                                                         tCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCg                                    >  1:567746/1‑93 (MQ=255)
                                                           ggTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCa                                  <  1:1556916/93‑1 (MQ=255)
                                                             tGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATc                                <  2:5574772/93‑1 (MQ=255)
                                                              gACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCt                               <  2:6362061/93‑1 (MQ=255)
                                                                 cTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTAc                            <  1:4356896/93‑1 (MQ=255)
                                                                 cTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTAc                            <  2:3008541/93‑1 (MQ=255)
                                                                 cTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTAc                            >  1:5955389/1‑93 (MQ=255)
                                                                  tATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACt                           <  2:2690000/93‑1 (MQ=255)
                                                                  tATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACt                           >  1:1298362/1‑93 (MQ=255)
                                                                        gcgAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTTAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCAtt                     >  1:7502479/1‑93 (MQ=255)
                                                                        gcgAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCAtt                     <  1:7382519/93‑1 (MQ=255)
                                                                          gAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGa                   <  2:3635225/93‑1 (MQ=255)
                                                                          gAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGa                   <  2:2454706/93‑1 (MQ=255)
                                                                          gAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGa                   >  1:393897/1‑93 (MQ=255)
                                                                            aCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATg                 >  1:2890404/1‑93 (MQ=255)
                                                                            aCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATg                 <  2:6063360/93‑1 (MQ=255)
                                                                             cTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGc                <  1:1441337/93‑1 (MQ=255)
                                                                             cTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGc                >  2:69791/1‑93 (MQ=255)
                                                                              tCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCg               <  2:4806256/93‑1 (MQ=255)
                                                                               ccGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGa              >  1:6524574/1‑93 (MQ=255)
                                                                                cGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGAc             >  2:6897712/1‑93 (MQ=255)
                                                                                 ggCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACg            >  1:5000365/1‑93 (MQ=255)
                                                                                    gAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGccc         >  2:2915871/1‑93 (MQ=255)
                                                                                     aGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcg        <  2:4611730/93‑1 (MQ=255)
                                                                                     aGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcg        >  2:5317674/1‑93 (MQ=255)
                                                                                      gcgcGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcga       <  2:2129559/93‑1 (MQ=255)
                                                                                          gACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAg   <  1:6399887/93‑1 (MQ=255)
                                                                                           |                                                                                            
GAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGC  >  NC_000913/4130307‑4130490

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