Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,225,050 A→G E1215G (GAG→GGG)  metH → B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,225,0500AG99.2% 404.9 / ‑4.6 127E1215G (GAG→GGG) metHB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (70/56);  total (71/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

Only 100 of 128 total aligned reads displayed.

CGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAG  >  NC_000913/4224960‑4225138
                                                                                          |                                                                                        
cGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGgc                                                                                        >  1:173764/1‑93 (MQ=255)
  gACAGCAATTACTACGCTGTAGCAAAAATCCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGgcgc                                                                                      <  2:1533238/93‑1 (MQ=255)
  gACAGAAAGTACTCCGTTGAAGCAAAAATCCAGCGCGATAAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGgcgc                                                                                      <  1:900124/93‑1 (MQ=255)
   aCAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgct                                                                                     <  1:1204883/93‑1 (MQ=255)
   aCAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgct                                                                                     >  2:4484166/1‑93 (MQ=255)
     aGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAGGTTGGGCgctgg                                                                                   >  1:1276065/1‑93 (MQ=255)
      gCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggt                                                                                  >  2:1279865/1‑92 (MQ=255)
      gCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggc                                                                                  <  2:1588489/93‑1 (MQ=255)
        aaGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggctg                                                                                <  2:2906276/93‑1 (MQ=255)
         aGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGctggctgg                                                                               >  1:726793/1‑93 (MQ=255)
         aGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGctggctgg                                                                               >  2:3427510/1‑93 (MQ=255)
          gTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCtggctggc                                                                              >  2:3642999/1‑93 (MQ=255)
          gTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCtggctggc                                                                              >  2:3702213/1‑93 (MQ=255)
          gTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCtggctggc                                                                              >  1:4182653/1‑93 (MQ=255)
            actacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAc                                                                            >  1:1093238/1‑93 (MQ=255)
             ctacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAcc                                                                           <  2:1998894/93‑1 (MQ=255)
              cacGCTGTAGCACAAATTCCGCGCGATCAGTTTGAAGTTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCg                                                                          <  1:1812085/92‑1 (MQ=255)
               acGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGa                                                                         >  1:2297482/1‑93 (MQ=255)
                cGCTGTAGCACAAATTCAGGGCGATGAGGGTGAAAGTTAGCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGaa                                                                        <  1:1893278/93‑1 (MQ=255)
                cGCTGTAGCAAAAATTCGGCGCGATCGGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGaa                                                                        <  1:1893279/93‑1 (MQ=255)
                 gCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAAt                                                                       >  2:2884569/1‑93 (MQ=255)
                  cTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATc                                                                      >  2:4055429/1‑93 (MQ=255)
                   tGTAGCACAATTTCAGCGCGATCAGGTGAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTTGCTGGCACCGAATCt                                                                     <  2:1884753/93‑1 (MQ=255)
                   tGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCt                                                                     >  1:3851146/1‑93 (MQ=255)
                     tAGCACAATTTAAGCGAGATCAGTTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgg                                                                   <  1:2668311/93‑1 (MQ=255)
                     tAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgg                                                                   <  1:1485805/93‑1 (MQ=255)
                       gCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                 >  2:1036797/1‑93 (MQ=255)
                       ccacaAATTCGGCGCGATCAGTTTGAAGTTTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                 <  1:2235769/92‑1 (MQ=255)
                        cacaAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGGATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGt                                                                >  2:2233497/1‑93 (MQ=255)
                          caAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTAt                                                              >  2:1118941/1‑93 (MQ=255)
                          acaaTGCAGCGCTATGAGGTTGAACATTAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTAt                                                              <  1:2580900/91‑1 (MQ=255)
                           aaaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATg                                                             >  2:3659385/1‑93 (MQ=255)
                           aaaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATg                                                             >  2:43858/1‑93 (MQ=255)
                            aaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                            >  1:2398501/1‑93 (MQ=255)
                            aaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                            >  2:628003/1‑93 (MQ=255)
                             aTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGag                                                           >  2:4491045/1‑92 (MQ=255)
                               tCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGGATCTGGGGTATGAcgg                                                         >  2:2704560/1‑92 (MQ=255)
                               tCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcgc                                                         <  1:3668157/93‑1 (MQ=255)
                               tCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcgc                                                         >  1:2075594/1‑93 (MQ=255)
                                cAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcgcg                                                        >  1:2887748/1‑93 (MQ=255)
                                  gcgcGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGa                                                      >  2:4412945/1‑93 (MQ=255)
                                  gcgcGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGa                                                      >  1:1635177/1‑93 (MQ=255)
                                   ggcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGAc                                                     <  1:4408131/92‑1 (MQ=255)
                                    gcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACt                                                    >  1:3165791/1‑93 (MQ=255)
                                    gcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACt                                                    >  1:3165790/1‑93 (MQ=255)
                                     cgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTg                                                   >  1:4519443/1‑93 (MQ=255)
                                      gATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGa                                                  >  1:1675493/1‑93 (MQ=255)
                                      gATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGa                                                  >  1:2224376/1‑93 (MQ=255)
                                      gATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGa                                                  >  1:535364/1‑93 (MQ=255)
                                       gccAGGTTGAAGATGAGCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                 <  2:1874184/91‑1 (MQ=255)
                                       aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                 <  1:4237229/93‑1 (MQ=255)
                                       aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                 >  1:4579741/1‑93 (MQ=255)
                                       aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                 >  1:481441/1‑93 (MQ=255)
                                         atggTTGAAAATTGTCGCCGCCGTAAAGGTATGTGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTc                                               <  1:4571703/91‑1 (MQ=255)
                                          aGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTca                                              >  2:2162695/1‑93 (MQ=255)
                                          aGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTca                                              >  1:1594584/1‑93 (MQ=255)
                                              tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACaaa                                          >  1:3691687/1‑93 (MQ=255)
                                                aaGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGAACCGAATCTGGGGTATGACGCGGACTGATTCACAAATc                                        >  2:4228758/1‑93 (MQ=255)
                                                 aGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCt                                       >  2:1310786/1‑93 (MQ=255)
                                                    ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTc                                    >  1:4319014/1‑93 (MQ=255)
                                                    ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTc                                    >  2:1512010/1‑93 (MQ=255)
                                                      tgccccGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCAc                                  <  1:1167720/90‑1 (MQ=255)
                                                       tGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACt                                 >  1:4707766/1‑93 (MQ=255)
                                                        gCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGGATGACGCGGACTGATTCACAAATCTGTCACtt                                >  1:392527/1‑93 (MQ=255)
                                                          ccgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtttt                              <  2:2013930/93‑1 (MQ=255)
                                                          ccgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtttt                              >  1:1902397/1‑93 (MQ=255)
                                                          ccgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtttt                              >  2:2662231/1‑93 (MQ=255)
                                                          ccgccgTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtttt                              >  1:4703306/1‑93 (MQ=255)
                                                           cgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTc                             <  2:3340603/93‑1 (MQ=255)
                                                             ccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCt                           <  1:4009806/93‑1 (MQ=255)
                                                             ccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCt                           <  1:3667571/93‑1 (MQ=255)
                                                              cgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCtt                          <  2:829322/93‑1 (MQ=255)
                                                              cgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCtt                          >  1:1757700/1‑93 (MQ=255)
                                                               gTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTa                         <  1:602848/93‑1 (MQ=255)
                                                                tAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTac                        >  1:3083387/1‑93 (MQ=255)
                                                                 aaaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaca                       <  1:2446862/93‑1 (MQ=255)
                                                                   aGGTATGAGCGTTGCCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaac                     <  1:2968687/93‑1 (MQ=255)
                                                                   aGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCATAAATCTGTCACTTTTCCTTacaac                     <  2:4593545/93‑1 (MQ=255)
                                                                   aGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaac                     <  2:4287562/93‑1 (MQ=255)
                                                                    ggTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaaca                    >  1:2314410/1‑93 (MQ=255)
                                                                     gTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTAcaacaa                   <  1:3942729/93‑1 (MQ=255)
                                                                       aTGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaac                 <  1:461062/93‑1 (MQ=255)
                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                <  1:254960/93‑1 (MQ=255)
                                                                         gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg               <  1:2971321/93‑1 (MQ=255)
                                                                         gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg               <  1:990640/93‑1 (MQ=255)
                                                                         gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg               <  2:1459624/93‑1 (MQ=255)
                                                                          aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg              <  2:2384560/93‑1 (MQ=255)
                                                                          aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg              <  2:4372561/93‑1 (MQ=255)
                                                                            cGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgc            >  2:1290287/1‑93 (MQ=255)
                                                                             gTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcg           >  1:2931242/1‑93 (MQ=255)
                                                                             gTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcg           >  2:4633291/1‑93 (MQ=255)
                                                                             gTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcg           >  1:1602137/1‑93 (MQ=255)
                                                                              ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc          <  2:616394/93‑1 (MQ=255)
                                                                               tACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCt         <  2:3761086/93‑1 (MQ=255)
                                                                                gccGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTc        <  2:1110149/92‑1 (MQ=255)
                                                                                aCCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTc        <  2:3859691/93‑1 (MQ=255)
                                                                                  cGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa      <  1:342606/93‑1 (MQ=255)
                                                                                   gAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt     >  1:4431909/1‑93 (MQ=255)
                                                                                     aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa   >  2:3895762/1‑93 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg  <  1:4246845/93‑1 (MQ=255)
                                                                                          |                                                                                        
CGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAG  >  NC_000913/4224960‑4225138

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGC  >  NC_000913/4224957‑4225141
                                                                                             |                                                                                           
ACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGG                                                                                             >  SRR4302544.4718137/1‑94 (MQ=60)
 ccgggacagaaggtactaccttgaAGCACAATTTCAGCGCGATTAGGTTGAAGATTAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGG                                                                                            <  SRR4302544.3504661/70‑1 (MQ=60)
 CCCGGACAGAAAGTACTACGCTTTAGCACAAATTCAGCGCGATCAGGTGAAAGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGG                                                                                            <  SRR4302544.1225287/94‑1 (MQ=60)
  CCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGC                                                                                           >  SRR4302544.175478/1‑94 (MQ=60)
    GGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGC                                                                                         >  SRR4302544.2265498/1‑94 (MQ=60)
     GACAGAAAGTACTCCGTTGAAGCAAAAATCCAGCGCGATAAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCT                                                                                        <  SRR4302544.909029/94‑1 (MQ=60)
      ACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTG                                                                                       <  SRR4302544.1216820/94‑1 (MQ=60)
       CAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAGGTTGGGCGCTGG                                                                                      >  SRR4302544.1288629/1‑94 (MQ=60)
          CAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTG                                                                                   >  SRR4302544.3738868/1‑94 (MQ=60)
           AAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGG                                                                                  >  SRR4302544.733925/1‑94 (MQ=60)
            AGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGC                                                                                 >  SRR4302544.4224075/1‑94 (MQ=60)
              TACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAC                                                                               >  SRR4302544.1104115/1‑94 (MQ=60)
                 TACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGA                                                                            >  SRR4302544.2319158/1‑94 (MQ=60)
                 CACGCTGTAGCACAAATTCCGCGCGATCAGTTTGAAGTTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGA                                                                            <  SRR4302544.1829779/94‑1 (MQ=60)
                   CGCTGTAGCACAAATTCAGGGCGATGAGGGTGAAAGTTAGCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAAT                                                                          <  SRR4302544.1911571/94‑1 (MQ=60)
                   CGCTGTAGCAAAAATTCGGCGCGATCGGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAAT                                                                          <  SRR4302544.1911572/94‑1 (MQ=60)
                     CTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCT                                                                        >  SRR4302544.3889373/1‑94 (MQ=60)
                        TAGCACAATTTAAGCGAGATCAGTTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGG                                                                     <  SRR4302544.2693889/94‑1 (MQ=60)
                        TAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGG                                                                     <  SRR4302544.1500485/94‑1 (MQ=60)
                          CCACAAATTCGGCGCGATCAGTTTGAAGTTTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGT                                                                   <  SRR4302544.2256794/94‑1 (MQ=60)
                             acaatgcagcgctatgAGGTTGAACATTAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATG                                                                <  SRR4302544.2605704/78‑1 (MQ=60)
                             CAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATG                                                                >  SRR4302544.4323748/1‑94 (MQ=60)
                              AAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGA                                                               >  SRR4302544.2421249/1‑94 (MQ=60)
                               aatccggcgcaaccagtgtgaagatttcGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAC                                                              <  SRR4302544.4187534/66‑1 (MQ=60)
                                 TTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGC                                                            >  SRR4302544.2095540/1‑94 (MQ=60)
                                  TCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCG                                                           >  SRR4302544.2915591/1‑94 (MQ=60)
                                  TCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCG                                                           <  SRR4302544.3704736/94‑1 (MQ=60)
                                    AGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGA                                                         >  SRR4302544.1651273/1‑94 (MQ=60)
                                      GGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACT                                                       <  SRR4302544.4452650/94‑1 (MQ=60)
                                      GGCGATCAGGTGGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTTGCACCGAATCTGGGGTATGACGCGGACT                                                       <  SRR4302544.1373836/94‑1 (MQ=60)
                                      CGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACT                                                       >  SRR4302544.3195914/1‑94 (MQ=60)
                                      CGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACT                                                       >  SRR4302544.3195913/1‑94 (MQ=60)
                                       gcaaagagaTTGAAGTTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTG                                                      <  SRR4302544.4502543/85‑1 (MQ=60)
                                       GCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTG                                                      >  SRR4302544.4564697/1‑94 (MQ=60)
                                        CGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGA                                                     >  SRR4302544.2245303/1‑94 (MQ=60)
                                        CGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGA                                                     >  SRR4302544.1691994/1‑94 (MQ=60)
                                        CGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGA                                                     >  SRR4302544.540660/1‑94 (MQ=60)
                                         GATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAT                                                    >  SRR4302544.4625355/1‑94 (MQ=60)
                                         GATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAT                                                    >  SRR4302544.486152/1‑94 (MQ=60)
                                          ATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATT                                                   <  SRR4302544.4280416/94‑1 (MQ=60)
                                            atggttgaaaattgtcgCCGCCGTAAAGGTATGTGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCA                                                 <  SRR4302544.4617267/77‑1 (MQ=60)
                                            CAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCA                                                 >  SRR4302544.1610287/1‑94 (MQ=60)
                                             AGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCAC                                                >  SRR4302544.3102109/1‑94 (MQ=60)
                                                TTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAA                                             >  SRR4302544.3728472/1‑94 (MQ=60)
                                                 TGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAAT                                            >  SRR4302544.4206234/1‑94 (MQ=60)
                                                      ATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTC                                       >  SRR4302544.4363012/1‑94 (MQ=60)
                                                      ATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTC                                       >  SRR4302544.3674203/1‑94 (MQ=60)
                                                         tgcCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACT                                    <  SRR4302544.1179364/91‑1 (MQ=60)
                                                         ATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACT                                    >  SRR4302544.4754282/1‑94 (MQ=60)
                                                         ATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACT                                    >  SRR4302544.1777946/1‑94 (MQ=60)
                                                          TGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGGATGACGCGGACTGATTCACAAATCTGTCACTT                                   >  SRR4302544.396507/1‑94 (MQ=60)
                                                            CCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTT                                 >  SRR4302544.1920779/1‑94 (MQ=60)
                                                            CCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTT                                 >  SRR4302544.4749792/1‑94 (MQ=60)
                                                                CCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT                             >  SRR4302544.1774902/1‑94 (MQ=60)
                                                                CCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT                             <  SRR4302544.4049695/94‑1 (MQ=60)
                                                                CCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT                             <  SRR4302544.3704145/94‑1 (MQ=60)
                                                                  GTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTAC                           <  SRR4302544.608761/94‑1 (MQ=60)
                                                                  GTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTAC                           >  SRR4302544.3112760/1‑94 (MQ=60)
                                                                   TAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCCTACA                          >  SRR4302544.2046251/1‑94 (MQ=60)
                                                                    AAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAA                         <  SRR4302544.2470190/94‑1 (MQ=60)
                                                                      AGGTATGAGCGTTGCCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACA                       <  SRR4302544.2997158/94‑1 (MQ=60)
                                                                      AGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACA                       >  SRR4302544.2336310/1‑94 (MQ=60)
                                                                      AGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACA                       <  SRR4302544.2738841/94‑1 (MQ=60)
                                                                        GTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAA                     <  SRR4302544.3981902/94‑1 (MQ=60)
                                                                          ATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACA                   <  SRR4302544.465594/94‑1 (MQ=60)
                                                                           TGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAG                  <  SRR4302544.257664/94‑1 (MQ=60)
                                                                           TGAGCGTTACCGAAGTTGGCCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAG                  <  SRR4302544.4662253/94‑1 (MQ=60)
                                                                            GAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGG                 <  SRR4302544.1688616/94‑1 (MQ=60)
                                                                            GAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGG                 >  SRR4302544.4586456/1‑94 (MQ=60)
                                                                            GAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGG                 <  SRR4302544.1000394/94‑1 (MQ=60)
                                                                            GAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGG                 <  SRR4302544.2999807/94‑1 (MQ=60)
                                                                               CGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCG              >  SRR4302544.2959410/1‑94 (MQ=60)
                                                                               CGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCG              >  SRR4302544.1617934/1‑94 (MQ=60)
                                                                                     CGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAT        <  SRR4302544.346146/94‑1 (MQ=60)
                                                                                     CGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAT        >  SRR4302544.4476621/1‑94 (MQ=60)
                                                                                      GAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATG       >  SRR4302544.546672/1‑94 (MQ=60)
                                                                                         GTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGC    <  SRR4302544.4290140/94‑1 (MQ=60)
                                                                                         GTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGC    <  SRR4302544.2472985/94‑1 (MQ=60)
                                                                                           TGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGC  >  SRR4302544.1635338/1‑94 (MQ=60)
                                                                                           TGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGC  >  SRR4302544.1046904/1‑94 (MQ=60)
                                                                                             |                                                                                           
ACCCGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGC  >  NC_000913/4224957‑4225141

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