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 521 238.2 12299412 96.0% 11807435 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,059,022 T→C G51G (GGT→GGC)  sbm → methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,059,0220TC99.5% 745.5 / NA 214G51G (GGT→GGC) sbmmethylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (104/108);  total (105/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

TGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTT  >  NC_000913/3058932‑3059111
                                                                                          |                                                                                         
  tCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCTATCTCGATAATCTGGAGGTGACAGGCAccc                                                                                       <  1:1881877/93‑1 (MQ=255)
  tCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCAccc                                                                                       <  2:540130/93‑1 (MQ=255)
  gcGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCAccc                                                                                       <  2:662180/92‑1 (MQ=255)
    gACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCtt                                                                                     >  1:5142902/1‑93 (MQ=255)
     aCTCGCTGGTTCATCAAACCGCGGAAGGTATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTc                                                                                    >  2:4235117/1‑93 (MQ=255)
     aCTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTc                                                                                    >  1:276592/1‑93 (MQ=255)
     aCTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTc                                                                                    <  2:4426101/93‑1 (MQ=255)
       tCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCt                                                                                  >  2:1856714/1‑93 (MQ=255)
       tCGCTGGTTCACCAAACCGCTGAAGGGATCCCCATCAAGCCGCTTTATACCGAAGCCCATCTCGATAATCTGGAGGTGACAGGCACCCTTCCt                                                                                  <  1:1735925/93‑1 (MQ=255)
         gCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTgg                                                                                >  2:5557448/1‑93 (MQ=255)
           tGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGtt                                                                              >  1:651778/1‑93 (MQ=255)
            ggTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGttt                                                                             >  2:636317/1‑93 (MQ=255)
             gTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTg                                                                            >  1:3184399/1‑93 (MQ=255)
              ttCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgc                                                                           <  2:1235591/93‑1 (MQ=255)
              taaatcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgc                                                                           <  1:5817117/90‑1 (MQ=255)
               tcatcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgcc                                                                          <  2:3844367/93‑1 (MQ=255)
                 atcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgccgc                                                                        >  2:4702082/1‑93 (MQ=255)
                 atcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgccgc                                                                        <  1:2563201/93‑1 (MQ=255)
                    aaaCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCt                                                                     <  1:3030252/93‑1 (MQ=255)
                    aaaCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCt                                                                     <  1:2271171/93‑1 (MQ=255)
                     aaCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTa                                                                    <  2:5800883/93‑1 (MQ=255)
                      aCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTAc                                                                   <  2:5162822/93‑1 (MQ=255)
                      aCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTAc                                                                   <  1:4332239/93‑1 (MQ=255)
                       ccGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACg                                                                  <  2:3730179/93‑1 (MQ=255)
                       ccGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACg                                                                  <  1:4109353/93‑1 (MQ=255)
                        cgcgGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGt                                                                 <  2:6126794/93‑1 (MQ=255)
                         gcgGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGtt                                                                <  1:963938/93‑1 (MQ=255)
                         gcgGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGtt                                                                <  1:1127031/93‑1 (MQ=255)
                           ggAGGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCg                                                              <  1:5398757/93‑1 (MQ=255)
                           ggAAGGGATCGCCATCAAGCCGCTGTATGCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCg                                                              <  2:4121555/93‑1 (MQ=255)
                           ggAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCg                                                              <  1:5544207/93‑1 (MQ=255)
                           ggAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCg                                                              <  2:3345213/93‑1 (MQ=255)
                            gAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGt                                                             >  2:2284246/1‑93 (MQ=255)
                            gAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGt                                                             <  1:1311204/93‑1 (MQ=255)
                            gAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGt                                                             >  2:3016560/1‑93 (MQ=255)
                            gAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGt                                                             >  1:4166052/1‑93 (MQ=255)
                             aaGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTg                                                            <  2:2473227/93‑1 (MQ=255)
                             aaGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTg                                                            <  1:179507/93‑1 (MQ=255)
                              aGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTgg                                                           >  1:5707312/1‑93 (MQ=255)
                                ggATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTggtc                                                         >  1:228021/1‑91 (MQ=255)
                                ggATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGcc                                                         >  2:1100635/1‑93 (MQ=255)
                                 gATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGccc                                                        >  1:761513/1‑93 (MQ=255)
                                 gATCGCCATCAAGCCGCTGTATACCAAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGccc                                                        <  1:4826117/93‑1 (MQ=255)
                                   tCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCcgc                                                      <  1:3195207/93‑1 (MQ=255)
                                    cGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCcgcg                                                     >  1:3749872/1‑93 (MQ=255)
                                     gCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCcgcgg                                                    >  2:1677176/1‑92 (MQ=255)
                                      ccATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTg                                                   >  1:541407/1‑93 (MQ=255)
                                       cATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGc                                                  <  2:3625042/93‑1 (MQ=255)
                                        aTCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGcc                                                 >  2:3040391/1‑93 (MQ=255)
                                        aTCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGcc                                                 <  2:4550379/93‑1 (MQ=255)
                                         tCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCa                                                >  1:6082360/1‑93 (MQ=255)
                                          cAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCAc                                               <  2:162156/93‑1 (MQ=255)
                                           aaGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGAACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCCTGCCACt                                              >  1:2805039/1‑93 (MQ=255)
                                            aGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTa                                             <  1:5876174/93‑1 (MQ=255)
                                              ccGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATg                                           <  1:4703486/93‑1 (MQ=255)
                                               cGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                          <  2:3951987/93‑1 (MQ=255)
                                               cGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                          >  2:1047694/1‑93 (MQ=255)
                                                gCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGta                                         <  2:2941985/93‑1 (MQ=255)
                                                gCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGta                                         >  2:1546828/1‑93 (MQ=255)
                                                 cTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGtat                                        >  2:286483/1‑93 (MQ=255)
                                                  tGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACtatatata                                       >  1:5535195/1‑93 (MQ=255)
                                                  tGTATACCGAAGCCGATCTAGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGtata                                       >  1:5448517/1‑93 (MQ=255)
                                                    tataCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATAcc                                     >  1:3286686/1‑93 (MQ=255)
                                                     ataCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTCCCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCg                                    <  1:2191320/93‑1 (MQ=255)
                                                     ataCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCg                                    >  1:3839491/1‑93 (MQ=255)
                                                       aCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGcc                                  >  2:5617267/1‑93 (MQ=255)
                                                        ccGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGccc                                 >  2:254459/1‑93 (MQ=255)
                                                           aaGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAAc                              <  2:5852815/93‑1 (MQ=255)
                                                            aGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAAcc                             <  2:3988824/93‑1 (MQ=255)
                                                             gCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCg                            >  2:3483489/1‑93 (MQ=255)
                                                             accGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCg                            <  1:5179086/92‑1 (MQ=255)
                                                               cGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTg                          <  2:4893710/93‑1 (MQ=255)
                                                                gATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTgg                         >  1:5804201/1‑93 (MQ=255)
                                                                 aTCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGGTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGa                        >  1:716947/1‑93 (MQ=255)
                                                                    tcGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCa                     <  1:4030022/93‑1 (MQ=255)
                                                                     cGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACGGCGGACCAt                    >  2:680477/1‑93 (MQ=255)
                                                                      gATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATc                   <  2:2571337/93‑1 (MQ=255)
                                                                      cagcaTCTGGAGGAGACGGGCACCCTTCCTGGTTTTCCGCCCTGCGTACGTGGCCCGCGGGCCACTATGTATACCGCCCAACCGGGGACCATc                   <  1:5086363/89‑1 (MQ=255)
                                                                       aTAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATcc                  >  2:2172677/1‑93 (MQ=255)
                                                                       aTAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATcc                  <  2:5130924/93‑1 (MQ=255)
                                                                        tAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCg                 <  1:1235379/93‑1 (MQ=255)
                                                                        tAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCg                 >  2:2381080/1‑93 (MQ=255)
                                                                        tAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCg                 <  2:1467759/93‑1 (MQ=255)
                                                                        tAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCg                 <  1:5265541/93‑1 (MQ=255)
                                                                        tAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCg                 >  2:3187809/1‑93 (MQ=255)
                                                                        tAATCGGGAGGTGACAGGCCCCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCg                 <  1:1270686/93‑1 (MQ=255)
                                                                         aaTCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCgc                >  2:4303350/1‑92 (MQ=255)
                                                                         aaTCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGt                >  1:2692521/1‑93 (MQ=255)
                                                                          aTCTGGAGGTGACAGGCACCCTTCCTGGTTTGACGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATATCGTCCAACCGCGGAACAACCGTc               >  1:3890707/1‑93 (MQ=255)
                                                                           tCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCa              <  1:2087549/93‑1 (MQ=255)
                                                                           ccTGGAGGGGACAGGCACCCTTCCTGACTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCa              <  1:1372957/92‑1 (MQ=255)
                                                                             tGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGt            >  2:5024163/1‑93 (MQ=255)
                                                                             tGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGt            <  2:118284/93‑1 (MQ=255)
                                                                              ggAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTa           >  2:2728808/1‑93 (MQ=255)
                                                                               gAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTAt          >  1:4217970/1‑93 (MQ=255)
                                                                                 ggTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGc        <  1:2474943/93‑1 (MQ=255)
                                                                                 ggTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGc        <  2:1333227/93‑1 (MQ=255)
                                                                                  gTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATgat       >  1:2162848/1‑91 (MQ=255)
                                                                                  gTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCt       >  2:1605234/1‑93 (MQ=255)
                                                                                   tGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCACGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTg      <  1:4492287/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTT  >  NC_000913/3058932‑3059111

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCA  >  NC_000913/3058930‑3059115
                                                                                            |                                                                                             
ACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCA                                                                                              >  SRR4302532.489778/1‑94 (MQ=60)
    TCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCTATCTCGATAATCTGGAGGTGACAGGCACCCT                                                                                          <  SRR4302532.1944267/94‑1 (MQ=60)
     CGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTT                                                                                         >  SRR4302532.5319826/1‑94 (MQ=60)
      GACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTC                                                                                        >  SRR4302532.286174/1‑94 (MQ=60)
      GACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTC                                                                                        >  SRR4302532.4910430/1‑94 (MQ=60)
       cccccCTGGTTCATCAACCCGCGGAAGGGATCTCCAGCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCGCCCTTCC                                                                                       <  SRR4302532.3020373/89‑1 (MQ=60)
        CTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCT                                                                                      >  SRR4302532.5616902/1‑94 (MQ=60)
         TCGCTGGTTCATCAAACCGCTGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTG                                                                                     <  SRR4302532.1794356/94‑1 (MQ=60)
         TCGCTGGTTCACCAAACCGCTGAAGGGATCCCCATCAAGCCGCTTTATACCGAAGCCCATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTG                                                                                     <  SRR4302532.1794355/94‑1 (MQ=60)
         CCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTG                                                                                     <  SRR4302532.3897124/94‑1 (MQ=60)
            CTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTT                                                                                  >  SRR4302532.674370/1‑94 (MQ=60)
             TGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTT                                                                                 <  SRR4302532.3803892/94‑1 (MQ=60)
              GGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTG                                                                                >  SRR4302532.5131366/1‑94 (MQ=60)
              GGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTG                                                                                >  SRR4302532.3290137/1‑94 (MQ=60)
              GGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTG                                                                                <  SRR4302532.1934912/94‑1 (MQ=60)
              GGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTG                                                                                >  SRR4302532.5733440/1‑94 (MQ=60)
              GGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTACTGGTTTC                                                                                >  SRR4302532.2223625/1‑94 (MQ=60)
               GTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGC                                                                               <  SRR4302532.1596360/94‑1 (MQ=60)
                taaATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCC                                                                              <  SRR4302532.6019164/91‑1 (MQ=60)
                TTCATCAACCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCC                                                                              <  SRR4302532.6019163/94‑1 (MQ=60)
                TTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCC                                                                              >  SRR4302532.872405/1‑94 (MQ=60)
                 TCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCG                                                                             >  SRR4302532.2844765/1‑94 (MQ=60)
                   ATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCC                                                                           <  SRR4302532.2648739/94‑1 (MQ=60)
                      AAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTA                                                                        <  SRR4302532.2346542/94‑1 (MQ=60)
                      AAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTA                                                                        <  SRR4302532.3130262/94‑1 (MQ=60)
                        ACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACG                                                                      <  SRR4302532.4478593/94‑1 (MQ=60)
                         CCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGT                                                                     <  SRR4302532.4247290/94‑1 (MQ=60)
                           GCGGAGGGGCTCGCCATCAAGCCGCTGTATGCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTGC                                                                   <  SRR4302532.996829/94‑1 (MQ=60)
                           GCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTC                                                                   >  SRR4302532.1355726/1‑94 (MQ=60)
                           GCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTC                                                                   <  SRR4302532.996830/94‑1 (MQ=60)
                           GCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTC                                                                   <  SRR4302532.1164636/94‑1 (MQ=60)
                             GGAGGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGT                                                                 <  SRR4302532.5729684/94‑1 (MQ=60)
                             GGAGGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGT                                                                 <  SRR4302532.5585826/94‑1 (MQ=60)
                             GGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGT                                                                 >  SRR4302532.3119324/1‑94 (MQ=60)
                             GGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGT                                                                 >  SRR4302532.4305759/1‑94 (MQ=60)
                             GGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGT                                                                 <  SRR4302532.5736634/94‑1 (MQ=60)
                              GAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTG                                                                <  SRR4302532.1355075/94‑1 (MQ=60)
                              GAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTG                                                                <  SRR4302532.2581870/94‑1 (MQ=60)
                               AAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGG                                                               <  SRR4302532.185549/94‑1 (MQ=60)
                               AAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGG                                                               >  SRR4302532.5906289/1‑94 (MQ=60)
                                 GGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGTC                                                             >  SRR4302532.235811/1‑94 (MQ=60)
                                  GGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCC                                                            >  SRR4302532.787738/1‑94 (MQ=60)
                                   GATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCG                                                           >  SRR4302532.4151157/1‑94 (MQ=60)
                                   GATCGCCATCAAGCCGCTGTATACCAAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCG                                                           <  SRR4302532.4992260/94‑1 (MQ=60)
                                     TCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCG                                                         <  SRR4302532.3301300/94‑1 (MQ=60)
                                     TCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCG                                                         >  SRR4302532.3875938/1‑94 (MQ=60)
                                       GCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTG                                                       >  SRR4302532.560267/1‑94 (MQ=60)
                                       GCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTG                                                       >  SRR4302532.3296633/1‑94 (MQ=60)
                                       GCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTG                                                       >  SRR4302532.5562028/1‑94 (MQ=60)
                                       CCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTG                                                       <  SRR4302532.2742435/94‑1 (MQ=60)
                                          ATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCA                                                    >  SRR4302532.6292394/1‑94 (MQ=60)
                                          ATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGGGACAGGCACCCTTCCTGGTTTGCCGCCGTACGTTCttttcccgcgtgccc                                                    >  SRR4302532.2351656/1‑79 (MQ=60)
                                            CAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGAACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCCTGCCACT                                                  >  SRR4302532.2899006/1‑94 (MQ=60)
                                              AGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTAT                                                <  SRR4302532.6079750/94‑1 (MQ=60)
                                              AGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTAT                                                <  SRR4302532.6079752/94‑1 (MQ=60)
                                                ccgctttttcccgcAGCCGATCACGATAATCTGGAGGTGACAGGCATCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGT                                              <  SRR4302532.4863944/80‑1 (MQ=60)
                                                CCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGT                                              <  SRR4302532.4863945/94‑1 (MQ=60)
                                                  cgTGTATAACGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTAT                                            <  SRR4302532.4862170/92‑1 (MQ=60)
                                                  GCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTAT                                            >  SRR4302532.741118/1‑94 (MQ=60)
                                                  GCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGACCGCGTGCCCCTATGTAT                                            >  SRR4302532.4319062/1‑94 (MQ=60)
                                                   CTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATA                                           <  SRR4302532.1256467/94‑1 (MQ=60)
                                                   CTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATA                                           >  SRR4302532.5716362/1‑94 (MQ=60)
                                                   CTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATA                                           >  SRR4302532.5482317/1‑94 (MQ=60)
                                                   CTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATATATA                                           >  SRR4302532.5727294/1‑94 (MQ=60)
                                                   CTGTATACCGAAGCCGATCTAGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATA                                           >  SRR4302532.5637357/1‑94 (MQ=60)
                                                     GTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACC                                         >  SRR4302532.3396013/1‑94 (MQ=60)
                                                     GTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACC                                         <  SRR4302532.447946/94‑1 (MQ=60)
                                                      TATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCG                                        >  SRR4302532.3969127/1‑94 (MQ=60)
                                                       ATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTCCCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGC                                       <  SRR4302532.2264356/94‑1 (MQ=60)
                                                          CCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCA                                    >  SRR4302532.226091/1‑94 (MQ=60)
                                                          CCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCA                                    >  SRR4302532.226092/1‑94 (MQ=60)
                                                             AAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACC                                 <  SRR4302532.5847091/94‑1 (MQ=60)
                                                               GCCGATCTCGAGAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGG                               >  SRR4302532.4211652/1‑94 (MQ=60)
                                                               ACCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGT                               <  SRR4302532.5357450/94‑1 (MQ=60)
                                                                 CGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGG                             >  SRR4302532.5416500/1‑94 (MQ=60)
                                                                 CGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGG                             >  SRR4302532.6005898/1‑94 (MQ=60)
                                                                  GATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGA                            >  SRR4302532.4006151/1‑94 (MQ=60)
                                                                  GATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGA                            >  SRR4302532.665315/1‑94 (MQ=60)
                                                                  GATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGA                            >  SRR4302532.1332901/1‑94 (MQ=60)
                                                                  GATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGGTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGA                            >  SRR4302532.741433/1‑94 (MQ=60)
                                                                   ATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGAC                           <  SRR4302532.559806/94‑1 (MQ=60)
                                                                   ATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGAC                           <  SRR4302532.1278272/94‑1 (MQ=60)
                                                                   ATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGAC                           <  SRR4302532.4596863/94‑1 (MQ=60)
                                                                   ATCTCGAGAATCTGGAGGTCACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGAC                           <  SRR4302532.559808/94‑1 (MQ=60)
                                                                    TCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACC                          <  SRR4302532.3183817/94‑1 (MQ=60)
                                                                    TCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACC                          <  SRR4302532.775050/94‑1 (MQ=60)
                                                                    TCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACC                          <  SRR4302532.930470/94‑1 (MQ=60)
                                                                     CTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCA                         <  SRR4302532.2864418/94‑1 (MQ=60)
                                                                     CTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCA                         >  SRR4302532.6272036/1‑94 (MQ=60)
                                                                      TCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCAT                        <  SRR4302532.4165808/94‑1 (MQ=60)
                                                                       CGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGCGCCACTATGTATACCGCCCAACCGTGGACCATC                       >  SRR4302532.5277267/1‑94 (MQ=60)
                                                                        cagcATCTGGAGGAGACGGGCACCCTTCCTGGTTTTCCGCCCTGCGTACGTGGCCCGCGGGCCACTATGTATACCGCCCAACCGGGGACCATCC                      <  SRR4302532.5261243/90‑1 (MQ=57)
                                                                        GATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCC                      <  SRR4302532.3966452/94‑1 (MQ=60)
                                                                        GATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCC                      <  SRR4302532.2585708/94‑1 (MQ=60)
                                                                        GATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCC                      <  SRR4302532.5346563/94‑1 (MQ=60)
                                                                          TAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGT                    >  SRR4302532.2782956/1‑94 (MQ=60)
                                                                          TAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGT                    <  SRR4302532.1276700/94‑1 (MQ=60)
                                                                          TAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGT                    <  SRR4302532.5447122/94‑1 (MQ=60)
                                                                          TAATCGGGAGGTGACAGGCCCCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGT                    <  SRR4302532.1313192/94‑1 (MQ=60)
                                                                           AATCTGGAGGTGACAGGCACCCTTCCTGGTTTGACGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATATCGTCCAACCGCGGAACAACCGTC                   >  SRR4302532.4021974/1‑94 (MQ=60)
                                                                            ATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCA                  <  SRR4302532.3274330/94‑1 (MQ=60)
                                                                            ATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGGTCGTGGCCCGCGTGCCACTATGTATACCGCCCACCCGTGGACCATCCGTCA                  >  SRR4302532.6239658/1‑94 (MQ=60)
                                                                             TCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAG                 <  SRR4302532.2480689/94‑1 (MQ=60)
                                                                             TCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAG                 <  SRR4302532.4347906/94‑1 (MQ=60)
                                                                             TCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAG                 <  SRR4302532.2157155/94‑1 (MQ=60)
                                                                             CCTGGAGGGGACAGGCACCCTTCCTGACTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAG                 <  SRR4302532.1418954/94‑1 (MQ=60)
                                                                              CTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGT                >  SRR4302532.3361307/1‑94 (MQ=60)
                                                                               TGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTA               <  SRR4302532.2296246/94‑1 (MQ=60)
                                                                                GGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTAT              <  SRR4302532.1755996/94‑1 (MQ=60)
                                                                                GGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTAT              >  SRR4302532.4359543/1‑94 (MQ=60)
                                                                                 GAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATG             <  SRR4302532.4282840/94‑1 (MQ=60)
                                                                                 GAGGTGACAGGCGCCCTTCCTGGTTTGCAGCCCTACGTTCGTGGCCCGCGTGCCACTATAAATACCGCCAAACCGTGGACCATCCGTCAGTATG             <  SRR4302532.4920081/94‑1 (MQ=60)
                                                                                   GGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCT           <  SRR4302532.2557366/94‑1 (MQ=60)
                                                                                   GGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGAT           >  SRR4302532.2235177/1‑94 (MQ=60)
                                                                                    GTGACAGGCACCCTGCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCCCTATTTATTCCGCCCAACCGTGGACCATCCGTCAGTATGCTG          <  SRR4302532.766111/94‑1 (MQ=60)
                                                                                     TGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCACGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGG         <  SRR4302532.4645651/94‑1 (MQ=60)
                                                                                      GACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGT        <  SRR4302532.5997361/94‑1 (MQ=60)
                                                                                      GACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGT        <  SRR4302532.2026448/94‑1 (MQ=60)
                                                                                         AGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTT     >  SRR4302532.5301305/1‑94 (MQ=60)
                                                                                         AGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTT     >  SRR4302532.1176476/1‑94 (MQ=60)
                                                                                          GGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTT    <  SRR4302532.5080339/94‑1 (MQ=60)
                                                                                            GACCCTTCCTGGTTTGCCGCCCTACGTGCGTGGCCCGCGTTCCACTATGTATCCCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCA  <  SRR4302532.67064/94‑1 (MQ=60)
                                                                                            CACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCA  <  SRR4302532.2019592/94‑1 (MQ=60)
                                                                                            |                                                                                             
ACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCA  >  NC_000913/3058930‑3059115

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 2 ≤ ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑